Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,064	T→A	12.0%	intergenic (+38/‑148)	rpsB → / → tsf	30S ribosomal protein S2/elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,064	0	T	A	12.0%	181.9 / 13.8	65	intergenic (+38/‑148)	rpsB/tsf	30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base A (4/4);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785922‑1786213                                                                                                                                                |                                                                                                                                                      cTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGcccc                                                                                                                                                >  2:1904239/1‑150 (MQ=255)   gTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCTAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCt                                                                                                                                              <  1:961236/150‑1 (MQ=255)      cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCtttt                                                                                                                                           <  1:784881/150‑1 (MQ=255)       aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                          <  2:2085604/150‑1 (MQ=255)         gACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGc                                                                                                                                                   >  1:948515/1‑139 (MQ=255)         gACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGc                                                                                                                                                   <  2:948515/139‑1 (MQ=255)                aTGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTg                                                                                                                                 >  2:1122310/1‑150 (MQ=255)                   gcggACGCAGTAATCGGCGGTGGCGAGAACGTTGTAGGCGGGAGCGTAGTTTACGCTGAGGAGGCGGCTGCAGGTGCTGAGGGGGAATAGG‑CGTTCGGGTTTACGTGGCACGCAAAAAGGGGGCTTTGCCCTCTTTTTGCCAGCTGGaaa                                                                                                                              <  2:607474/147‑1 (MQ=255)                   gCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGaaa                                                                                                                              <  1:1545453/150‑1 (MQ=255)                   gCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGAGGCTTTGCCCCCTTTTTGCCCCCTTGaaa                                                                                                                              >  1:2346678/1‑150 (MQ=255)                    cTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAAt                                                                                                                             >  1:1563834/1‑150 (MQ=255)                     tGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATc                                                                                                                            >  2:1888522/1‑150 (MQ=255)                       aCGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCct                                                                                                                          <  1:1122310/150‑1 (MQ=255)                        cGCACTCATGCGCGCCCGCCACCCCGTCGCCGGCCATGCCCAGCTTCATGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                         <  2:1841659/150‑1 (MQ=255)                           aGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCGCCTTGAAATCctgctg                                                                                                                      >  1:949464/1‑150 (MQ=255)                                tCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGc                                                                                                                 >  1:742287/1‑150 (MQ=255)                                        cGCAACAACGTTGCCGGCGACACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                         <  2:1447811/150‑1 (MQ=255)                                         gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGcccc                                                                                                        >  1:1786915/1‑150 (MQ=255)                                         gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGAGCACTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGcccc                                                                                                        >  2:1489191/1‑150 (MQ=255)                                           aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCAt                                                                                                               <  1:708648/141‑1 (MQ=255)                                           aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCAt                                                                                                               >  2:708648/1‑141 (MQ=255)                                           aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCg                                                                                                             <  2:1443144/143‑1 (MQ=255)                                           aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCg                                                                                                             >  1:1443144/1‑143 (MQ=255)                                            acaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGca                                                                                                     <  2:742287/150‑1 (MQ=255)                                              aacGTTGCAGGCGGCACCGAAGTTTACGCTGAAGAGGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCCAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTGGAAATCCTGTTGTCAGCATCGGCCCCGcatc                                                                                                   <  1:791707/150‑1 (MQ=255)                                                  ttGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCGGCTGTCAGCATCGGCCCCGCATCATgg                                                                                               <  2:1168029/150‑1 (MQ=255)                                                         cggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        <  1:795501/150‑1 (MQ=255)                                                         cggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        <  1:795510/150‑1 (MQ=255)                                                             aCCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaa                                                                                    <  1:836781/150‑1 (MQ=255)                                                              ccGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaat                                                                                   <  1:909254/150‑1 (MQ=255)                                                                gAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatga                                                                                 <  2:2112801/150‑1 (MQ=255)                                                                   gTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTCCTTGGTCCATATCG‑CTGGTTGGCACGCAAAAAGGGGACTGGGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                              <  1:1741766/150‑1 (MQ=255)                                                                        cGCTGAAGAAGCGGCTGCACGTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCg                                                                         >  1:1207638/1‑150 (MQ=255)                                                                         gCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCgg                                                                        >  1:248383/1‑150 (MQ=255)                                                                           tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAt                                                                      >  2:2223945/1‑150 (MQ=255)                                                                           tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAACTTGACATCCTGCTGTCACCATCCTCCCCTCAGCAGGGGCCGGCTGAAATGAAAACAGCGGAt                                                                      >  2:2230609/1‑150 (MQ=255)                                                                            gaagaaGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCGACGCATCATGGGCCTGCTGAAATGAAAACAGCGGAtt                                                                     >  1:318100/1‑150 (MQ=255)                                                                              agaagCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTg                                                                   <  1:1098760/150‑1 (MQ=255)                                                                                  gCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGa                                                               >  1:91895/1‑150 (MQ=255)                                                                                   cGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGaa                                                              <  1:1888522/150‑1 (MQ=255)                                                                                        gCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGaa                                                         <  2:1563834/150‑1 (MQ=255)                                                                                           ccTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGc                                                      >  1:738335/1‑150 (MQ=255)                                                                                                tgctgAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                      >  2:684289/1‑63 (MQ=25)                                                                                                tgctgAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                      <  1:684289/65‑3 (MQ=21)                                                                                                   tgANTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  2:1498813/1‑150 (MQ=255)                                                                                                   tgAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  1:1237238/1‑150 (MQ=255)                                                                                                         aTAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACggg                                        >  1:835540/1‑150 (MQ=255)                                                                                                                 cTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACATTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATg                                 <  2:2388515/150‑1 (MQ=255)                                                                                                                  tAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATgg                                <  2:1786915/150‑1 (MQ=255)                                                                                                                     cGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACTGCGGATTTGTATAGTTTGACGCCCCTGATTGACTGGTGGAATGGTTg                             <  2:1207638/150‑1 (MQ=255)                                                                                                                          aCTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGTAATGGTTGaaaaa                        >  2:817617/1‑150 (MQ=255)                                                                                                                            ttGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACt                      >  2:741187/1‑150 (MQ=255)                                                                                                                              ggCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACttt                    <  2:248383/150‑1 (MQ=255)                                                                                                                               ggtggCAAAAAGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                                 <  1:799912/100‑1 (MQ=38)                                                                                                                               ggtggCAAAAAGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                                 >  2:799912/5‑104 (MQ=38)                                                                                                                               gCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc                   >  2:1265986/1‑150 (MQ=255)                                                                                                                                  tgtaaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                <  2:2346678/147‑1 (MQ=255)                                                                                                                                     aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAtt                                                            >  1:379646/1‑103 (MQ=255)                                                                                                                                     aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAtt                                                            <  2:379646/103‑1 (MQ=255)                                                                                                                                     aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  1:1247925/1‑150 (MQ=255)                                                                                                                                      aaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAgg            >  2:559065/1‑150 (MQ=255)                                                                                                                                       aaaGGGGGCTTTGCCCCCTTTTTGCCGCCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGa           <  2:949464/150‑1 (MQ=255)                                                                                                                                           ggggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt       >  1:2312866/1‑150 (MQ=255)                                                                                                                                           ggggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt       <  2:91895/150‑1 (MQ=255)                                                                                                                                           ggggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt       <  2:351057/150‑1 (MQ=255)                                                                                                                                                tttGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  >  2:545261/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      CTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGG‑CGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785922‑1786213

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