Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,415	A→G	K271K (AAA→AAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,415	0	A	G	98.6%	232.8 / ‑4.9	73	K271K (AAA→AAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/34);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAC  >  NC_000913/4620181‑4620660                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      gtgCGAAATCGCCCGTGAAGAGCTGTCCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTg                                                                                                                                                                                                                                       <  2:363562/251‑1 (MQ=255) gtgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACt                                                                                                                                                                                                                                        <  1:231169/250‑1 (MQ=255)          cGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAc                                                                                                                                                                                                                               <  1:168044/250‑1 (MQ=255)           gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAcgc                                                                                                                                                                                                                             >  1:156971/1‑251 (MQ=255)                  aaGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCTGTGTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGTGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGttt                                                                                                                                                                                                                      <  2:388344/251‑1 (MQ=255)                        tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                                                                                >  1:386476/1‑251 (MQ=255)                             aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                                                                                           <  1:294656/251‑1 (MQ=255)                                  cggcTACAACATCGGCCGTGTTATCGCTCGCCCGTTTATCGGCGACAAAGCCGGTAAATTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACAGCGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGCGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTACCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGAAGCCACCAGCaa                                                                                                                                                                                                        >  2:293571/1‑249 (MQ=255)                                    gcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                                                                                                                                                                                      >  1:315824/1‑249 (MQ=255)                                     cTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                                                                                                                                                                                      >  1:34968/1‑248 (MQ=255)                                         aacatCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTAACGGTAACCGTCACGAACTGGCTGTTGAGCCGCCAGCACCGCCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGATTCTGTCGGCAAACTTGCAGACATCTACGCCACCTGCGGCATCACCACGAAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACACCACCATCAAAGAGCTGa                                                                                                                                                                                                >  1:286307/4‑250 (MQ=255)                                              cgggcgttttATCGCGCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGACACGTCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAaagaag                                                                                                                                                                                          <  2:46657/243‑1 (MQ=255)                                                      ttATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGAt                                                                                                                                                                                  >  2:167342/1‑251 (MQ=255)                                                        agcGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCAGCTAGCACCGACCGTGCTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTTGTTTCTTTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAATAAAGTGAAAGCAACTGGCCTGGACGCGTTTTTTGACGCCTCCATCAAAGAGATGAAAGAAGCGGGTGATaa                                                                                                                                                                                <  1:605667/249‑1 (MQ=255)                                                          cgcctgtcctttaaacggtgACAAAGCCGGTAACTTCCAGCTTACTGGTAACCATCACGCCCTGGCTGTTGATCCGCCAGCACCGATCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCATGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACCTCGGTATCACCCAGAAAGTGAAAGTAACTGGCCTGGACGCGCTGTTTGACGCCACCTTCAAAGAGATGAAATAAGCCGGTGATatca                                                                                                                                                                              <  2:78078/232‑4 (MQ=255)                                                          cGCTCGTCCGTTTATCGGTGACAAAGCCGTTAACTTCCAGCGTACCTGTAACCTTCACGACCTGGCTTGTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCAAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGTGATGAAAGAAGCGGGTGATAac                                                                                                                                                                               <  1:400470/250‑1 (MQ=255)                                                          cGCTCGTCCGTTTACCGGTATCAAAGCCGGTAGCTTCCAGCGTGCCGGAAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGGAACTGGTTGATGAAAAACACGGTCAGGTGTTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTTGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAaca                                                                                                                                                                              <  1:285304/251‑1 (MQ=255)                                                              cGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCAt                                                                                                                                                                          >  1:562635/1‑251 (MQ=255)                                                                tCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCg                                                                                                                                                                        >  1:591204/1‑251 (MQ=255)                                                                    tttATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCATACATCTACGCCAACTGCGGCATCACCAAGAAAGAGAAGGCAACTGGCCTGGACGCGCTGTTTGACGCCCCCATCAAAGAGATGAAAGAATCggtgaataaccccaacgcctt                                                                                                                                                                    >  2:256456/1‑232 (MQ=255)                                                                      tATCGGTGAAAAAGCCGGAAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTc                                                                                                                                                                   <  2:156971/250‑1 (MQ=255)                                                                         cggtgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCAcc                                                                                                                                                                >  1:387987/5‑250 (MQ=255)                                                                         cggtgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCAcc                                                                                                                                                                >  1:573123/5‑250 (MQ=255)                                                                                    ccGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACACCACCCTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGAc                                                                                                                                                    >  2:261468/1‑251 (MQ=255)                                                                                     cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTAATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACt                                                                                                                                                   <  1:169618/251‑1 (MQ=255)                                                                                      ggTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACtt                                                                                                                                                  >  1:267888/1‑251 (MQ=255)                                                                                       gTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTc                                                                                                                                                 <  1:48938/251‑1 (MQ=255)                                                                                             cccctcgtcccGGTAACCGTCACGACCTGGCTGTGGAGCCGCCAGCACCGTCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGGCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGActct                                                                                                                                           <  2:387987/242‑1 (MQ=255)                                                                                               catcGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCACCAGCACCGACCGTGCTGCAGAAACTGGTTGCTGAAAAACACGATCAGGAGGCTTCTGTCGGTAACAATGCAGACATCTACGCCAACTGCGGTATCACCACGAAAGGGAAAGCAACTGGCCTGTACGCGCTGTTTGACGCCACCATCAAAGAGATACAAGAAGCTCGAAATAACACCCACGACTTCCCCCAATTCGTTGACTTCGActctac                                                                                                                                         >  2:60062/4‑249 (MQ=255)                                                                                                  cGTACCGGTAACCGTCACGACCGGGCTGTTGAGCCGCCAGCACCGACAGTGCAGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGGCGGAAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCTGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                                                                                      <  2:265870/251‑1 (MQ=255)                                                                                                    tACCGATAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGCAGCGGGTGATAACACCATCGTCTTCACCAACATCGTTGACTTCGCCTCTTCCTggg                                                                                                                                    >  2:311708/1‑251 (MQ=255)                                                                                                                    gACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGt                                                                                                                    >  1:546951/1‑251 (MQ=255)                                                                                                                      ccTGGCTGTTGAGCCGCCAGCACCGAACGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCCGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACGCTTCCTGGGGCCACCGTCGCGACATCg                                                                                                                  >  2:457009/1‑251 (MQ=255)                                                                                                                               tGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTg                                                                                                                                                                                                                                  >  2:640468/1‑130 (MQ=255)                                                                                                                               tGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTg                                                                                                                                                                                                                                  <  1:640468/130‑1 (MQ=255)                                                                                                                                         gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                                                                               <  1:399864/251‑1 (MQ=255)                                                                                                                                           accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                                                                                  >  1:189881/1‑246 (MQ=255)                                                                                                                                           accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                                                                             <  1:251500/251‑1 (MQ=255)                                                                                                                                           accgaccgTGCTGCAGAAACTGGGTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                                                                                  <  2:189881/246‑1 (MQ=255)                                                                                                                                                 cgTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGAGATACCACCATCGTCTTCACCAACTTCGTTGACTACGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGagcttt                                                                                       >  2:505988/1‑246 (MQ=255)                                                                                                                                                       gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                  <  1:440735/250‑1 (MQ=255)                                                                                                                                                       gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                  <  1:256456/250‑1 (MQ=255)                                                                                                                                                                ggTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTg                                                                         <  1:277375/250‑1 (MQ=255)                                                                                                                                                                 gTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGTCCGCCGTCTGcc                                                                       <  2:90776/251‑1 (MQ=255)                                                                                                                                                                    gatgaAAAACACGGACAGGTGGATTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACACTTCCTGGGGCCACCGTCGCTCCGTCGCCGGTTATGCTGCGGGCCTGGAGTTGTACCACCGCCGCCTGCCggg                                                                    >  2:276264/1‑250 (MQ=255)                                                                                                                                                                        aaaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTg                                                                >  1:408635/1‑251 (MQ=255)                                                                                                                                                                         aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCTGGTGATAACACCATCGTCTTCACCAACGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGa                                                               <  1:98925/251‑1 (MQ=255)                                                                                                                                                                          aaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGACAGTGAAAGCAACTGGCCTGGCCGAGCTGTTTGACGCCACCAACAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCCCCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCAGCCGCCGTCTGCCGGAGCTGAt                                                              >  2:601797/1‑251 (MQ=255)                                                                                                                                                                                    aGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGGGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTctgctg                                                    >  1:373386/1‑251 (MQ=255)                                                                                                                                                                                     ggaggtTACTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCGGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTTCGGGTCTGGAGCTGGTCGACCGAAGTCTGCCCGTGGTGATGtctttacttc                                                   >  2:571722/4‑244 (MQ=255)                                                                                                                                                                                      gtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcg                                                  <  1:675156/251‑1 (MQ=255)                                                                                                                                                                                      gtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCGGCCGGAGCTGATGTCTCTGCTgcg                                                  >  1:329055/1‑251 (MQ=255)                                                                                                                                                                                       tggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcgc                                                 >  1:519052/1‑251 (MQ=255)                                                                                                                                                                                           ttCTGTAGGTAAAATTGCAGACATATACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGCCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGCTCTGGAGCTGTTCGACCGCAGTCTGACGGAGCTGATGGCTCTGCTGCGCGATg                                             >  2:292990/1‑251 (MQ=255)                                                                                                                                                                                            tCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATg                                             <  1:311708/250‑1 (MQ=255)                                                                                                                                                                                              tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCTGCTGCGCGATgacg                                          >  2:240837/1‑251 (MQ=255)                                                                                                                                                                                                  ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                      <  1:9619/251‑1 (MQ=255)                                                                                                                                                                                                     aaaaTTGCAGACATCTACGCCAACTTCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCAGCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                                                                                      <  2:404623/152‑1 (MQ=255)                                                                                                                                                                                                     aaaaTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                                                                                      >  1:404623/1‑152 (MQ=255)                                                                                                                                                                                                      aatttGCGGACATCTACGTCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACATCATCGTCTTCACCAACTCCGTTGACTTCGTCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACGGTTCGACCGCCGTCTGCCGGAGCTGATGGCTCGGCCGCGCGATGACAGCATCCTg                                  >  2:290173/4‑251 (MQ=255)                                                                                                                                                                                                        aTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAt                                >  1:342506/1‑251 (MQ=255)                                                                                                                                                                                                                  tCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCt                      >  1:548329/1‑251 (MQ=255)                                                                                                                                                                                                                  tCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGGTTGACGCCACCATAAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTAGTTGACTTCGACTCTTCCTGGGGCCACTGTCGCGACGTCGCCGGTTATGCGGCGTGTCTGGAGCTGTTCGACCGCCTTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGTTCATCACTGCt                      >  2:535314/1‑251 (MQ=255)                                                                                                                                                                                                                     aCGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAACGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGctgc                                                                                                     >  2:254972/1‑166 (MQ=255)                                                                                                                                                                                                                     aCGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGctgc                                                                                                     <  1:254972/169‑4 (MQ=255)                                                                                                                                                                                                                     aCGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCActgc                       <  1:678706/247‑4 (MQ=255)                                                                                                                                                                                                                     aCGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCt                      <  1:398506/248‑1 (MQ=255)                                                                                                                                                                                                                     aCGCCAACTGCGGTAGCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGAc                   <  1:635946/251‑1 (MQ=255)                                                                                                                                                                                                                     aCGCCAACGGCGGTATCGACAAGAAATAGAAAGAAACGGGCCTGGACGCGCTGTGTGACGCCACCAACAAAGAGATGAAAGAAGCGGGTGATAACCCCATCGTCTTCACCAACTTCGTTGACTTCGACGCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGTTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAGCCTGATCCTctctgctg                     <  2:267888/249‑8 (MQ=255)                                                                                                                                                                                                                             cgCGGTATCACCAAGAAAGTGAAAGCCACTGGCCTGGACGCGCTGTTTGACTCCACCAGAAAAGAGATGAAAGAAGCGGTTGATACCACCATCGTCTTCACCATCGTCGTTGACTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCCCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTg           <  2:557871/250‑1 (MQ=255)                                                                                                                                                                                                                               cGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATAGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCAGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTgtg         >  1:402279/1‑249 (MQ=255)                                                                                                                                                                                                                                    tCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTc                                                                           >  1:482204/1‑180 (MQ=255)                                                                                                                                                                                                                                      aCCAAGAAATTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGAc  <  1:160310/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      GTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAC  >  NC_000913/4620181‑4620660

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