Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	83,587	A→C	10.9%	intergenic (‑58/+35)	leuA ← / ← leuL	2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	83,587	0	A	C	10.9%	286.6 / 23.5	100	intergenic (‑58/+35)	leuA/leuL	2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (47/42);  new base C (8/3);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATAT  >  NC_000913/83441‑83735                                                                                                                                                   |                                                                                                                                                      tGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                                    >  1:484474/1‑150 (MQ=255) tGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                                    >  2:318787/1‑150 (MQ=255)       ttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCATTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCgg                                                                                                                                              <  2:1563753/150‑1 (MQ=255)       ttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCATTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCgg                                                                                                                                              <  2:1571164/150‑1 (MQ=255)       ttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCgg                                                                                                                                              <  1:643692/150‑1 (MQ=255)            tttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttt                                                                                                                                        <  2:1985232/150‑1 (MQ=255)            tttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tttt                                                                                                                                        <  1:701213/150‑1 (MQ=255)              tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tttttt                                                                                                                                      >  1:1045839/1‑150 (MQ=255)              tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tttttt                                                                                                                                      >  1:190671/1‑150 (MQ=255)               cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttg                                                                                                                                     >  1:1494605/1‑150 (MQ=255)               cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttg                                                                                                                                     <  1:79571/150‑1 (MQ=255)                 cacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCGTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtt                                                                                                                                   <  1:79901/150‑1 (MQ=255)                  acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgttt                                                                                                                                  <  2:169506/150‑1 (MQ=255)                  acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgttt                                                                                                                                  >  2:29752/1‑150 (MQ=255)                  acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgttt                                                                                                                                  <  2:1449866/150‑1 (MQ=255)                     cAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGAc                                                                                                                               >  2:1537018/1‑150 (MQ=255)                           tGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTg                                                                                                                         <  1:1892034/150‑1 (MQ=255)                             ccTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCt                                                                                                                       <  2:1555297/150‑1 (MQ=255)                             ccTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGCCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGAGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCt                                                                                                                       >  1:1907862/1‑150 (MQ=255)                                gTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGc                                                                                                                    >  1:1159576/1‑150 (MQ=255)                                 tAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCt                                                                                                                   <  1:380009/150‑1 (MQ=255)                                  aaCGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCtt                                                                                                                  >  1:1329177/1‑150 (MQ=255)                                  aaCGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCtt                                                                                                                  <  1:1537018/150‑1 (MQ=255)                                      ccgggtCACCGGCGCGCAATGGGGTATCGAAAATTATGCCTTGCCTGGTCATGGGTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                                                                              <  1:1562278/145‑1 (MQ=255)                                      ccTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                                                                              <  2:1569238/150‑1 (MQ=255)                                        tGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                                                                            >  2:1881480/1‑150 (MQ=255)                                         gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                                                                           >  2:320850/1‑150 (MQ=255)                                          ttCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTgg                                                                                                          <  1:1567749/150‑1 (MQ=255)                                               cGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGTTGGATGcc                                                                                                     <  2:572164/150‑1 (MQ=255)                                                  cgcgCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCt                                                                                                  <  2:654201/150‑1 (MQ=255)                                                   gcgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCta                                                                                                 >  1:247491/1‑149 (MQ=255)                                                    cgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCa                                                                                                >  1:281583/1‑150 (MQ=255)                                                      cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACt                                                                                              <  1:320850/150‑1 (MQ=255)                                                           tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCt                                                                                         <  2:1329177/150‑1 (MQ=255)                                                           tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCt                                                                                         <  1:1181308/150‑1 (MQ=255)                                                               aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcg                                                                                     <  2:960688/150‑1 (MQ=255)                                                                 cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                                                                   >  1:1572720/1‑150 (MQ=255)                                                                 cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                                                                   >  1:534196/1‑150 (MQ=255)                                                                 cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                                                                   >  1:551380/1‑150 (MQ=255)                                                                 cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                                                                   >  1:1707120/1‑150 (MQ=255)                                                                  gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCgcgc                                                                                  >  2:268560/1‑150 (MQ=255)                                                                    aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaa                                                                                >  2:480543/1‑150 (MQ=255)                                                                         atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaag                                                                           >  2:682695/1‑150 (MQ=255)                                                                         atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TATTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaac                                                                           >  1:900375/1‑149 (MQ=255)                                                                            aCTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATg                                                                        <  1:753290/150‑1 (MQ=255)                                                                             cTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                                       >  2:226136/1‑150 (MQ=255)                                                                             cTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGAGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATCCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                                       >  2:1542130/1‑150 (MQ=255)                                                                                 cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                   <  1:653431/150‑1 (MQ=255)                                                                                  tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                                                                  >  1:1238052/1‑150 (MQ=255)                                                                                   ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                                                                 <  2:247491/150‑1 (MQ=255)                                                                                     cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                               >  1:1413232/1‑150 (MQ=255)                                                                                       cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagt                                                             <  2:1960599/150‑1 (MQ=255)                                                                                       cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagt                                                             <  1:480543/150‑1 (MQ=255)                                                                                        aTGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                            <  1:1214170/150‑1 (MQ=255)                                                                                        aTGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTCGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                            <  2:1045839/150‑1 (MQ=255)                                                                                         tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtag                                                           <  2:1494605/150‑1 (MQ=255)                                                                                           gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagta                                                         <  2:801901/150‑1 (MQ=255)                                                                                            tttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagtag                                                        >  2:1429297/1‑150 (MQ=255)                                                                                               gggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAcc                                                     <  1:318787/150‑1 (MQ=255)                                                                                                ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGCCCg                                                    >  1:1533294/1‑150 (MQ=255)                                                                                                 gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCg                                                    >  2:178554/1‑150 (MQ=255)                                                                                                 gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                   >  2:193529/1‑150 (MQ=255)                                                                                                    cTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaa                                                >  1:1294989/1‑150 (MQ=255)                                                                                                      tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                              >  1:1691917/1‑150 (MQ=255)                                                                                                       gTCTCCTTAAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                             <  1:873655/150‑1 (MQ=255)                                                                                                        tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCg                                            >  1:775355/1‑150 (MQ=255)                                                                                                         ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                   >  1:1595002/1‑142 (MQ=255)                                                                                                         ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                   <  2:1595002/142‑1 (MQ=255)                                                                                                         ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGa                                           >  1:1285808/1‑150 (MQ=255)                                                                                                           ctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                         >  2:873450/1‑150 (MQ=255)                                                                                                               tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                     >  1:1331150/1‑150 (MQ=255)                                                                                                                  agCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCCCTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                  >  1:1518467/1‑150 (MQ=255)                                                                                                                  agCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TGTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACGCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGAGAAAGCGAACGATGTGa                                  >  1:1511421/1‑150 (MQ=255)                                                                                                                   gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAg                                 <  2:1238052/150‑1 (MQ=255)                                                                                                                   gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAg                                 >  1:985407/1‑150 (MQ=255)                                                                                                                      ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                              >  1:1485396/1‑150 (MQ=255)                                                                                                                         cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACg                                        <  1:103840/137‑1 (MQ=255)                                                                                                                         cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACg                                        >  2:103840/1‑137 (MQ=255)                                                                                                                             tcggttaTAAAAAAACCCGCGCAATGGGGGGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATc                       <  2:900375/144‑1 (MQ=255)                                                                                                                               gtcaaacaaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCCCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                     >  2:747780/8‑150 (MQ=39)                                                                                                                                ggCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGc                    >  1:1972611/1‑150 (MQ=255)                                                                                                                                ggCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGc                    <  2:1572720/150‑1 (MQ=255)                                                                                                                                    tAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa                <  2:1159576/150‑1 (MQ=255)                                                                                                                                    tAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa                >  2:1719123/1‑150 (MQ=255)                                                                                                                                     caaaaaaCCCGCGCCATTCCGCGCG‑TTTATTGATTGACTCCGTGCTGGCTTAATGCTGGATGCCGCTCAAACGTATACCGCGCAAAGAAGATGCGTATAGTAGTAATAGACAGATAAAGCGAACGATGTGAGTCAATAAATCAGCTCCAg               >  2:1425033/2‑150 (MQ=255)                                                                                                                                      aaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGa              >  2:228315/1‑150 (MQ=255)                                                                                                                                       aaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGAt             <  2:484474/150‑1 (MQ=255)                                                                                                                                        aaatCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg            <  1:1425033/150‑1 (MQ=255)                                                                                                                                        aaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                           <  2:251610/135‑1 (MQ=255)                                                                                                                                        aaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                           >  1:251610/1‑135 (MQ=255)                                                                                                                                        aaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg            >  1:199285/1‑150 (MQ=255)                                                                                                                                        aaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg            <  2:534196/150‑1 (MQ=255)                                                                                                                                         aaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGa           >  1:306986/1‑150 (MQ=255)                                                                                                                                          aaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa          >  1:452654/1‑150 (MQ=255)                                                                                                                                          aaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa          >  2:1667945/1‑150 (MQ=255)                                                                                                                                             ccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGc       <  1:682695/150‑1 (MQ=255)                                                                                                                                              cgcgCAATGGCGCGGGTTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGc       <  1:178554/150‑1 (MQ=255)                                                                                                                                              cgcgCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCg      >  1:937164/1‑150 (MQ=255)                                                                                                                                               gcgcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGa     >  2:796576/1‑150 (MQ=255)                                                                                                                                                  cAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGatat  >  2:21529/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATAT  >  NC_000913/83441‑83735

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