Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,917	T→C	I105I (ATT→ATC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,917	0	T	C	81.1%	95.5 / 14.2	53	I105I (ATT→ATC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (6/4);  new base C (22/21);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

TGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTG  >  NC_000913/4619674‑4620163                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       tGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGGACGCATGGGCGCACG‑AAATGGCATCCGGGAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGgtg                                                                                                                                                                                                                                                 >  1:95862/1‑251 (MQ=255)     gTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTGGTGCCAAAGGCGAAGCTGATAACGGGCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTGCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTccc                                                                                                                                                                                                                                             <  2:301046/251‑1 (MQ=255)        ggggCTGACTCCCTGTGTCATATCGCAGACGCTTGGGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAACCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGt                                                                                                                                                                                                                                          <  2:503623/251‑1 (MQ=255)              gACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGtt                                                                                                                                                                                                                                    >  1:136150/1‑251 (MQ=255)                  cccTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGCAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGCCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAAAGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATTC‑CCCGTCTGCTCCTTGGGAAATCACCGGTGTCCCGGTTCTGTTTg                                                                                                                                                                                                                                  >  1:31290/1‑249 (MQ=255)                     tGGGTCATATCGCAGAAGCTAGTGCAAGACGGGAAGCTGATAAAGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGATATCGGCGCGTACGCATGAGCGCACG‑GAACGTCATCCGGTACAGA‑ACACCAGTCAGGTCACTGGGGAAACGCCGGTGTCACGGGTCCGTTTGGGTgg                                                                                                                                                                                                                             >  2:246634/1‑251 (MQ=255)                           atatCGCAGAAGCTGGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGTTTCTACCTGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGatat                                                                                                                                                                                                                       <  1:223292/251‑1 (MQ=255)                                 cAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTcc                                                                                                                                                                                                                 >  1:258372/1‑251 (MQ=255)                                     aGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCCTCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGGGTCCCGGTTCTGTTTGCGGGGGGATATTTCTCCGATc                                                                                                                                                                                                             >  1:225887/1‑251 (MQ=255)                                      gCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATATGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCCTGGGCGCACG‑AAATTTCCTCCGGTCAAGATAC‑CCCGGCTTGGAACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATcc                                                                                                                                                                                                            >  2:448435/1‑250 (MQ=255)                                       cTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGCACGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACAACTGGGCGCACGGAAATG‑CATACGGTAAAAATGC‑CCAGGCTGGTACCTGGGAAATCGCCGGTGTCCCGGTTCAGTTTGAGTGGGGATATTTCCCCggccac                                                                                                                                                                                                           >  2:432244/1‑246 (MQ=255)                                            gCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAAGCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCTTGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaaa                                                                                                                                                                                                       <  1:358639/250‑1 (MQ=255)                                             ccAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAAc                                                                                                                                                                                                     <  1:14512/251‑1 (MQ=255)                                                aaGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAAG‑AAATGTCATCCGGGAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGATCTGTTTGAGTGGGGGTATTTCTCCGATCACGCAAACAGc                                                                                                                                                                                                  >  2:403273/1‑251 (MQ=255)                                                       aGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAcagatcgggagagcacc                                                                                                                                                                                           >  1:107776/1‑235 (MQ=255)                                                       aGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAc                                                                                                                                                                                                           <  2:107776/235‑1 (MQ=255)                                                          tGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCCTGGGCGCCCG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCCGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTTAGTGGGGATATTTCTCCGATCCCGAAAACAACTTcccccacg                                                                                                                                                                                        >  2:233809/1‑246 (MQ=255)                                                           gATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAag                                                                                                                                                                                        <  2:258372/250‑1 (MQ=255)                                                                 ggTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCTAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGCAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGGTCTGTTTGGGTGGGGATATTTTTCCGCACACGAAAACAGCTTCCCGCACAAgctgct                                                                                                                                                                                 >  2:219538/1‑251 (MQ=255)                                                                  gTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGctgctg                                                                                                                                                                                >  1:397852/1‑251 (MQ=255)                                                                      tAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAt                                                                                                                                                                             >  1:228601/1‑250 (MQ=255)                                                                          ggCCCGCTCAATCTTCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACt                                                                                                                                                                        <  1:489338/251‑1 (MQ=255)                                                                           gCCCGCTCAAACTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTg                                                                                                                                                                       <  2:82481/251‑1 (MQ=255)                                                                              cGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTc                                                                                                                                                                    <  1:103888/251‑1 (MQ=255)                                                                              cGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTc                                                                                                                                                                    <  1:46747/251‑1 (MQ=255)                                                                                cTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCg                                                                                                                                                                   >  1:30548/1‑250 (MQ=255)                                                                                   aaccTGCCAAATCTGACCCGGCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcg                                                                                                                                                               <  2:95862/248‑1 (MQ=255)                                                                                        gCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgc                                                                                                                                                              <  1:355067/247‑1 (MQ=255)                                                                                           aaaTCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATc                                                                                                                                                       <  1:442135/251‑1 (MQ=255)                                                                                                     cGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAGGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                                                                                                                                              >  1:320594/1‑250 (MQ=255)                                                                                                            ggCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                                                                                                                                                >  1:67418/1‑161 (MQ=255)                                                                                                            ggCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                                                                                                                                                <  2:67418/161‑1 (MQ=255)                                                                                                                        cacaCGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActct                                                                                                                          >  1:376775/1‑251 (MQ=255)                                                                                                                              aaGGTTCTACCGTTTTCATTCCGGCGGGAATGGCCGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGt                                                                                                                    <  1:47189/251‑1 (MQ=255)                                                                                                                                      aCCGGTTTCATTCCGGCGGGAAGGGACGGCAACGCTGAAGTAATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATCC‑CCCGGCGGGTAACTGGGAAATCGCCGGTGTCCCGTTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAt                                                                                                            <  2:162195/251‑1 (MQ=255)                                                                                                                                        cGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGGTACCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTc                                                                                                          <  1:109284/251‑1 (MQ=255)                                                                                                                                              cATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATc                                                                                                    >  1:370539/1‑251 (MQ=255)                                                                                                                                               aTTCAGGCGGGAATGGACCGCAAAGCTGAAGTTATCGGCGCGTACGCATGGGCGCCCG‑AAATTTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATTTTTCTCCGATCACGAAAACAGCTTCCCGCCAGAGCTGCTGGATAAACTTGTCGCACGCGCTAATCTGCCGGGTTACCTCGGTACCTTCCACTCTTCCGGTACGGTCATTCTGGATcc                                                                                                   >  1:163428/1‑250 (MQ=255)                                                                                                                                                 tCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCag                                                                                                  <  1:403273/250‑2 (MQ=255)                                                                                                                                                        ggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAAG‑AAAGGTCATCCGGTAAAGATAC‑CCCGGCTGGTCACTGGGAAAGCGCCGGGGTACAGGTTCTGTTTGAGTGGGGATATTTCTCCGATAACGAAAACAGATTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAAGCGATAATCTGCCGGGTTACCTCGGAAACTGCCACTCTTCCGGAACGGTAATTCTGGAACAGCTgggggg                                                                                          >  2:59256/1‑248 (MQ=255)                                                                                                                                                             ggACGGCAACGCTGAAGTTATCGGCGCGTACGCCTGGGCGCCCT‑GAATGGCATCCGATAAAATTAA‑CCCGTCTGGTCACTGGGAAATCGCCGGGGTACCGGCTCTGTTTGAGTGGGGATGTTTCTCCGATCACGAAAACAGCTTCCCGCAAGCGCTGCTGGATAAAATGGTCGAACGCGCTACTCTGCCGGGTTACCTCGGAAACTGACAATCTTACGGCACGCTCATTCAGGCTCAGCTGGGCGAAGAGc                                                                                     >  2:425158/1‑251 (MQ=255)                                                                                                                                                                           aaggcaTCGGCGCGGACGTATGGGCGCACG‑AAATGACCTCCGGAAAAGATAC‑CCCGTCTGGTAACTGGGAAATCGCCGGAGTCCCGGTTCAGTTTGAGTGGGGATATTTCTCCGATCACGAAAACCACTTCCAGCACGAGCTGCTGGATAAAATGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGACACTATTCAGGTACGGGAATTCTGGAACAGCGGGGCGTAGAACACATGAAAACCgg                                                                        >  2:3585/6‑250 (MQ=255)                                                                                                                                                                                     cgcgTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTTTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTG‑ATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAtttt                                                            <  1:219538/251‑1 (MQ=255)                                                                                                                                                                                      gcgaccGCATGGGCGCAAG‑AAATTGCACCCGGTAAAGATAC‑CCCGTCGGGTCCCTGGGAAACCGCCGGTGTCCCGGTTCTGTTTGAGGGGGAATATTTCCCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGGCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAtttt                                                            <  2:228601/246‑1 (MQ=255)                                                                                                                                                                                          ccGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCTGTTTTGTTTGAGTGGGGATATTTCTCCGATCCCGAAAACAGCTTCCCGCAAGAGCTTCGGGATAAACTGGTCGAACGCGCTAATCTGCCTGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCta                                                         <  1:492891/249‑1 (MQ=255)                                                                                                                                                                                                    gcACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCtata                                                       >  1:371693/1‑242 (MQ=255)                                                                                                                                                                                                    gcACG‑AAATGTCAATCGTAAAAGATAC‑CCCGTCTGTGCACGGGGAAATCGCCGGGGGCCCGGTTCTGTTGGAGTGGGGATATTTCCCCGTGAACAAAAACTGCTTCCCGAAAGAGCTGCTGGATAACCTGGGCGAACGGGCTAATCTGCGTGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACAGGAAAACCGGCAAGCCGATTTTCtata                                                       <  2:371693/242‑1 (MQ=255)                                                                                                                                                                                                                 aTCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTcc                                  >  2:99616/1‑251 (MQ=255)                                                                                                                                                                                                                     ggTAAAGATAC‑CCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTcc                                  >  1:7576/1‑247 (MQ=255)                                                                                                                                                                                                                              cc‑cccGTCTGTTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCCCCGATCACGAAAACAGCTCCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCa                     <  2:136150/250‑1 (MQ=255)                                                                                                                                                                                                                                                ggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATTAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTc    >  1:100512/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaaaTCGCCGGGGTCCCGGTTTTGTTTGAGTGGGGTTATTTCCCCGATCCCGAAAACAGCTTCCCGCAGGACCTGCTGGATAACCTGGTCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcct   <  2:320594/250‑10 (MQ=255)                                                                                                                                                                                                                                                  aaaTCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcctg  >  1:174773/1‑241 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACG‑AAATGTCATCCGGTAAAGATAC‑CCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTG  >  NC_000913/4619674‑4620163

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: ‑