Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	709,556	A→T	5.4%	K100* (AAG→TAG)	ribF →	bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	709,556	0	A	T	5.4%	268.2 / 3.7	89	K100* (AAG→TAG)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base T (3/2);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAACGCGTACGGCTTTACC  >  NC_002947/709409‑709705                                                                                                                                                     |                                                                                                                                                      gtggtgATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAg                                                                                                                                                     <  1:1816986/150‑1 (MQ=255)   ggtgATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCt                                                                                                                                                   >  1:1473483/1‑150 (MQ=255)    gtgATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGACAGGGTAGCCGGCG‑CGGAACTCGACCAGCGGGGGGCTGTGCCTGGCGTTCAACCAGCGTGGGGGCAAGccc                                                                                                                                                  >  1:696962/1‑148 (MQ=255)     tgATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCa                                                                                                                                                 >  2:1084148/1‑150 (MQ=255)     tgATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCa                                                                                                                                                 >  1:83398/1‑150 (MQ=255)       aTCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGt                                                                                                                                               >  2:1490353/1‑150 (MQ=255)       aTCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGt                                                                                                                                               >  1:1460122/1‑150 (MQ=255)          ttCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGATGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGcc                                                                                                                                            >  2:1148811/1‑150 (MQ=255)           tCGAACCACAACCACGCGAGTACTTAGCACCCGATACCGCGCCAGCACGACTGGCCCGGCTGCGCAAAACGAGCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGGGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCg                                                                                                                                           >  2:1161782/1‑150 (MQ=255)                ccacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGcc                                                                                                                                      >  2:364556/1‑150 (MQ=255)                ccacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGcc                                                                                                                                      >  1:1336551/1‑150 (MQ=255)                ccacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGcc                                                                                                                                      <  2:613273/150‑1 (MQ=255)                 cacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCTTGCTCAgacgct                                                                                                                                           >  2:906516/1‑139 (MQ=255)                 cacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCTTGCTCAgacgct                                                                                                                                           >  2:904398/1‑139 (MQ=255)                 cacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCTTGCTCAgacgct                                                                                                                                           <  1:904398/144‑6 (MQ=255)                 cacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCTTGCTCAgacgct                                                                                                                                           <  1:906516/144‑6 (MQ=255)                 cacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCCGTGCCGATGCCt                                                                                                                                     >  2:52896/1‑150 (MQ=255)                 cacaaccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCt                                                                                                                                     >  1:291869/1‑150 (MQ=255)                     cccacGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGt                                                                                                                                 <  1:52896/149‑1 (MQ=255)                        acGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACACGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCaa                                                                                                                              <  2:606878/150‑1 (MQ=255)                        acGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCaa                                                                                                                              <  2:610120/150‑1 (MQ=255)                         cgcgAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAg                                                                                                                             <  1:1664992/150‑1 (MQ=255)                         cgcgAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAg                                                                                                                             <  1:1786833/150‑1 (MQ=255)                          gcgAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAgg                                                                                                                            <  1:1035879/150‑1 (MQ=255)                          gcgAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAgg                                                                                                                            <  1:1165321/150‑1 (MQ=255)                               tACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGAATCACATTTAGCTCAGTGCCGATGCCTTCGTCAAGGCCATc                                                                                                                       >  2:1590923/1‑150 (MQ=255)                                  ttCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTg                                                                                                                    >  2:1710222/1‑150 (MQ=255)                                   tCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTgg                                                                                                                   <  1:349768/150‑1 (MQ=255)                                    cGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGt                                                                                                                  <  1:364556/150‑1 (MQ=255)                                    cGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGt                                                                                                                  <  1:1710222/150‑1 (MQ=255)                                      cACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCg                                                                                                                >  1:1194717/1‑150 (MQ=255)                                       aCCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGa                                                                                                               >  2:236775/1‑150 (MQ=255)                                        cccGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGTGCAAGCTCAGTGCCGATGCCTTCGTCaa                                                                                                                              <  2:67935/134‑1 (MQ=255)                                        cccGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCaa                                                                                                                              >  1:67935/1‑134 (MQ=255)                                        cccGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGAc                                                                                                              >  1:1840667/1‑150 (MQ=255)                                        cccGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGAc                                                                                                              <  1:236775/150‑1 (MQ=255)                                         ccGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGAGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACg                                                                                                             >  1:1368967/1‑150 (MQ=255)                                                cgcgCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTggg                                                                                                      >  1:1749189/1‑150 (MQ=255)                                                 gcgcCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgc                                                                                                     >  2:493561/1‑150 (MQ=255)                                                 gcgcCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgc                                                                                                     <  1:1166093/150‑1 (MQ=255)                                                 gcgcCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgc                                                                                                     <  2:317977/150‑1 (MQ=255)                                                    ccAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgcgtg                                                                                                  <  2:1107511/150‑1 (MQ=255)                                                     cAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgcgtgc                                                                                                 >  1:588268/1‑150 (MQ=255)                                                     cAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgcgtgc                                                                                                 <  1:1490353/150‑1 (MQ=255)                                                      aGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGgcgtgcg                                                                                                >  2:1166474/1‑150 (MQ=255)                                                             cTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCAtctc                                                                                         >  2:877623/1‑150 (MQ=255)                                                                 cccGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAg                                                                                     <  2:83398/150‑1 (MQ=255)                                                                  ccGGCTGCGCGACAAGATCGAGCTGCTGGCGAGCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGt                                                                                    >  2:1072163/1‑150 (MQ=255)                                                                  ccGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGt                                                                                    <  2:1761420/150‑1 (MQ=255)                                                                  ccGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGATGTGCCTGGCGTCCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGt                                                                                    >  2:1493047/1‑150 (MQ=255)                                                                  ccGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATAGGTAATCTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGACGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGt                                                                                    <  1:1645493/150‑1 (MQ=255)                                                                    ggCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTgg                                                                                  >  1:732814/1‑150 (MQ=255)                                                                        gcgcGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGa                                                                              <  2:1336551/150‑1 (MQ=255)                                                                        gcgcGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGa                                                                              <  2:121820/150‑1 (MQ=255)                                                                         cgcgACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGAt                                                                             <  1:621291/150‑1 (MQ=255)                                                                            gACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGAc                                                                          >  1:1033624/1‑150 (MQ=255)                                                                              cAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACtt                                                                        <  1:1158876/150‑1 (MQ=255)                                                                                aGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTcc                                                                      >  2:133301/1‑150 (MQ=255)                                                                                  aTCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGc                                                                    <  1:1148811/150‑1 (MQ=255)                                                                                     gAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTc                                                                 <  2:1180414/150‑1 (MQ=255)                                                                                     gAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTc                                                                 <  2:714632/150‑1 (MQ=255)                                                                                         tgctgGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCt                                                             <  1:1357449/150‑1 (MQ=255)                                                                                                gACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACcg                                                      >  1:1321126/1‑150 (MQ=255)                                                                                                 aCCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACcgc                                                     <  2:1522658/150‑1 (MQ=255)                                                                                                   cGAAGGC‑ATCGACCTGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTCCGGCAAGGCCAGCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCgcgc                                                   <  1:595379/150‑1 (MQ=255)                                                                                                   cGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCgcgc                                                   >  2:1387725/1‑150 (MQ=255)                                                                                                         c‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGa                                             <  2:1473483/150‑1 (MQ=255)                                                                                                         c‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGa                                             >  2:1855418/1‑150 (MQ=255)                                                                                                           aTCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGAc                                            >  2:473162/1‑150 (MQ=255)                                                                                                           aTCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCATGCGTCATCTCGAAGTGGGCGATGACATCACCCTCGGCCGCGCCCGCGCTGGCGAc                                            <  2:1460122/150‑1 (MQ=255)                                                                                                              gACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTc                                         <  1:473162/150‑1 (MQ=255)                                                                                                              gACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTc                                         <  1:1161782/150‑1 (MQ=255)                                                                                                                 cGGGTGCTGTGCCTGGCGTTCAACCAGTGTCTGAGCAAGCTCAGTGCCGAGGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCt                                      <  2:166152/150‑1 (MQ=255)                                                                                                                  gggTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCAtc                                                                                           <  1:107299/96‑1 (MQ=255)                                                                                                                  gggTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCAtc                                                                                           >  2:107299/1‑96 (MQ=255)                                                                                                                    gTGCGGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTc                                   <  2:1663542/150‑1 (MQ=255)                                                                                                                     tGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTcc                                  <  1:1855418/150‑1 (MQ=255)                                                                                                                       cTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGCTAAGGCCATCCTGGGCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTg                                <  2:1038809/150‑1 (MQ=255)                                                                                                                        tgtgCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCNTgg                               <  2:1194717/150‑1 (MQ=255)                                                                                                                         gtgCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGt                              <  2:291869/150‑1 (MQ=255)                                                                                                                            ccTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGa                           <  2:84884/150‑1 (MQ=255)                                                                                                                              tGGCGTTCAACCAGCGTCTGAGCCCGGACACGTCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAgg                         <  1:479660/150‑1 (MQ=255)                                                                                                                                   ttCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGCTGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGc                    >  2:591017/1‑150 (MQ=255)                                                                                                                                      aaCCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATATCGAAGTGGGCGATGACTCCTGCTTCCGCTGCGACGCCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAAc                 >  2:1058445/1‑150 (MQ=255)                                                                                                                                       aCCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAAGCCATAACGGTCGTCGGCCTGGGCGTGCGGGTGCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAAcg                >  1:37288/1‑150 (MQ=255)                                                                                                                                        ccAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAAcgc               >  1:1712266/1‑150 (MQ=255)                                                                                                                                          aGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAACGCGt             >  1:1064899/1‑150 (MQ=255)                                                                                                                                                tGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAACGCGTACGGCt       >  2:442497/1‑150 (MQ=255)                                                                                                                                                tGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAACGCGTACGGCt       >  2:1377717/1‑150 (MQ=255)                                                                                                                                                     aaGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAACGCGTACGGCTTTAcc  <  2:1840667/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTGGTGATCTTCGAACCACAACCACGCGAGTACTTCGCACCCGATACCGCGCCAGCGCGCCTGGCCCGGCTGCGCGACAAGATCGAGCTGCTGGCGACCGAAGGC‑ATCGACCGGGTGCTGTGCCTGGCGTTCAACCAGCGTCTGAGCAAGCTCAGTGCCGATGCCTTCGTCAAGGCCATCCTGGTCGACGGCCTGGGCGTGCGTCATCTCGAAGTGGGCGATGACTTCCGCTTCGGCTGCGACCGCGCTGGCGACTTCGCTTTCCTGGTCGAGGCCGGCAACGCGTACGGCTTTACC  >  NC_002947/709409‑709705

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