Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	872,690	A→T	15.1%	Y158N (TAC→AAC)	PP_0756 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	872,690	0	A	T	15.1%	161.4 / 19.0	66	Y158N (TAC→AAC)	PP_0756	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/24);  new base T (5/5);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GGTAAGCAGCCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGATGGTCTGGTCGACGTAGTA  >  NC_002947/872546‑872834                                                                                                                                                  |                                                                                                                                                 ggTAAGCAGCCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACAcc                                                                                                                                              >  1:329594/1‑149 (MQ=255)  gTAAGCAGCCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCg                                                                                                                                             >  1:241397/1‑149 (MQ=255)      gcagcCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCa                                                                                                                                           >  7:43464/1‑147 (MQ=255)        agcCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCGCCCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAg                                                                                                                                        <  1:87391/149‑1 (MQ=255)         gcCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGtt                                                                                                                                      <  8:64524/149‑1 (MQ=255)           cAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGa                                                                                                                                    >  6:134288/1‑149 (MQ=255)           cAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGa                                                                                                                                    >  4:137376/1‑149 (MQ=255)            aGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGAc                                                                                                                                   >  1:148423/1‑149 (MQ=255)            aGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGAc                                                                                                                                   >  5:146678/1‑149 (MQ=255)                    gcgcCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGcc                                                                                                                            <  1:29133/148‑1 (MQ=255)                              cTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGAc                                                                                                                 <  2:207887/149‑1 (MQ=255)                                        ggTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACgg                                                                                                       >  7:295026/1‑149 (MQ=255)                                               gATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTcc                                                                                                                      >  4:297534/1‑127 (MQ=255)                                               gATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTcc                                                                                                                      <  3:297534/127‑1 (MQ=255)                                               gATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTgg                                                                                                >  5:87755/1‑149 (MQ=255)                                                 tGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGGGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGATGGCCGGGTCCATGACGATGTCACGGTGGTTGGcc                                                                                              <  8:295026/149‑1 (MQ=255)                                                 tGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGc                                                                                               >  6:64025/1‑148 (MQ=255)                                                   aaGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCa                                                                                             <  2:329594/148‑1 (MQ=255)                                                   aaGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAg                                                                                            <  4:174950/149‑1 (MQ=255)                                                   aaGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAg                                                                                            >  8:110788/1‑149 (MQ=255)                                                    aGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAgg                                                                                           <  1:204861/149‑1 (MQ=255)                                                      ctctTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGaa                                                                                         <  4:104334/149‑1 (MQ=255)                                                       tctTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCTGCGTGGTACACCGCATAGTTGACGAAGGCCCGGTCCATGACGATGTCACGGTGGTTGTCCAGGaa                                                                                         >  5:268775/1‑148 (MQ=255)                                                            ttCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAgg                                                                                     >  3:358879/1‑147 (MQ=255)                                                              cAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAgg                                                                                 >  3:239125/1‑149 (MQ=255)                                                              cAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAg                                                                                  <  8:43464/148‑1 (MQ=255)                                                               aGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGGGGTACCCCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGc                                                                                <  8:151042/149‑1 (MQ=255)                                                                ggCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCg                                                                               >  7:34488/1‑149 (MQ=255)                                                                   gcATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTg                                                                            >  4:23586/1‑149 (MQ=255)                                                                     atcatGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCa                                                                                                                                           >  1:146739/1‑84 (MQ=255)                                                                     atcatGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCa                                                                                                                                           <  2:146739/84‑1 (MQ=255)                                                                       tatttCGCTTAAAAATGGGTTTTGCAGCAAATTGTCG‑CCAATGGCGTTGCTGGGCGTGGGCAGGCCGGCTTGGTACACCGTATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCg                                                                         <  8:54113/144‑1 (MQ=255)                                                                        atGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCgg                                                                        <  3:23586/148‑1 (MQ=255)                                                                        atGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCgg                                                                        <  8:34488/148‑1 (MQ=255)                                                                          gTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCt                                                                     >  6:257654/1‑149 (MQ=255)                                                                           tCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCtt                                                                    <  3:98347/149‑1 (MQ=255)                                                                             gCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCTATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCgg                                                                        <  5:183393/143‑1 (MQ=255)                                                                             gCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCTATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCgg                                                                        >  6:183393/1‑143 (MQ=255)                                                                             gCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCtt                                                                    >  2:7052/1‑147 (MQ=255)                                                                                   aaaGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTTGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGtt                                                            >  7:345306/1‑149 (MQ=255)                                                                                   aaaGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGtt                                                            >  8:224498/1‑149 (MQ=255)                                                                                   aaaGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGtt                                                            >  3:307842/1‑149 (MQ=255)                                                                                    aaGGGCTTTTGCTGCAAATTGTCG‑GCAATGGCGATGCGCTGCGTGGGCAGTCCGGCTTGGTACACAGCATAGTTGTCGATGGCCGGGTCCATGAAGATGTCACTGTGGTTGGCCAGGAACAGGTAGGCGGTGCGGGCGTTTAGCTGTTc                                                           <  6:268775/149‑1 (MQ=255)                                                                                    aaGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTc                                                           <  7:110788/149‑1 (MQ=255)                                                                                        gCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCcaca                                                       >  2:337720/1‑149 (MQ=255)                                                                                        gCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGAGGTACACCGCATAGTTGACGAAGGCCGGCTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGCTAGGCGGTGCCGGCCTTGAGCTGTTCcaca                                                       >  8:197589/1‑149 (MQ=255)                                                                                         cTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCacac                                                      <  4:239125/149‑1 (MQ=255)                                                                                          ttttGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACAcc                                                     >  7:36343/1‑149 (MQ=255)                                                                                            ttGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACAccc                                                    <  5:257654/148‑1 (MQ=255)                                                                                              gcagcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAg                                                 >  8:64502/1‑149 (MQ=255)                                                                                               cagcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACAGCCGAg                                                 >  4:10768/1‑148 (MQ=255)                                                                                               cagcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAg                                                                                            >  6:127468/1‑105 (MQ=255)                                                                                               cagcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAg                                                                                            <  5:127468/105‑1 (MQ=255)                                                                                               cagcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAg                                                 >  5:358911/1‑148 (MQ=255)                                                                                                agcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAg                                                 >  4:53182/1‑147 (MQ=255)                                                                                                agcaAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAg                                                 >  2:299592/1‑147 (MQ=255)                                                                                                      ttGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGAc                                         <  3:137376/149‑1 (MQ=255)                                                                                                             ccAATGGCGATGCGCGGCGTTGGCAGGCCGGCGTGGGACCCCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGACCAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTc                                   <  4:293200/149‑1 (MQ=255)                                                                                                              cAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCg                                  >  4:240297/1‑149 (MQ=255)                                                                                                                     gATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGgcgc                           >  8:325905/1‑149 (MQ=255)                                                                                                                          gcggcgTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGcgatc                       >  2:76381/1‑148 (MQ=255)                                                                                                                             gcgTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGTGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGAtg                   <  6:146678/149‑1 (MQ=255)                                                                                                                                  ggCAGGCCGGCGTGGTTCACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGAtggtctg              <  8:36343/149‑1 (MQ=255)                                                                                                                                       gCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGATGGTCTGGTCGa         >  7:145689/1‑149 (MQ=255)                                                                                                                                       gCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGGGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGATGGTCTGGTCGa         <  2:241397/149‑1 (MQ=255)                                                                                                                                              tGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGATGGTCTGGTCGACgtagta  <  7:64502/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GGTAAGCAGCCAGTTTTTCGCGCCGCCCGCTGATCGAACGGTGCACGATGAAGCTCTTGTTCAGGCGCATCATGTCGCTGACAAAGGGCTTTTGCAGCAAATTGTCG‑CCAATGGCGATGCGCGGCGTGGGCAGGCCGGCGTGGTACACCGCATAGTTGACGAAGGCCGGGTCCATGACGATGTCACGGTGGTTGGCCAGGAACAGGTAGGCGGTGCCGGCCTTGAGCTGTTCCACACCCGAGTAGGTGACGCCGTCGGTGGCGCGATCGATGGTCTGGTCGACGTAGTA  >  NC_002947/872546‑872834

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑