Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,160	T→C	G186G (GGT→GGC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,160	0	T	C	95.7%	218.2 / ‑0.4	70	G186G (GGT→GGC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (38/29);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGA  >  NC_000913/4619929‑4620387                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       ggTTCTGTTTGAGTGGGTATAGTTCTCCGATCACGAAAACAGCTTCCCGCAAGTGCTGCTGGATAAACTGTTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTCCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAAc                                                                                                                                                                                                                     <  1:61272/251‑1 (MQ=255)             ggggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCGGCCGGGTTACCTGGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCg                                                                                                                                                                                                         <  1:606614/249‑1 (MQ=255)                  gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGCTCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTAGATAAACTCTACGAACTGTGAGAAATCGCCCgg                                                                                                                                                                                                    >  2:169618/1‑250 (MQ=255)                  gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGt                                                                                                                                                                                                    >  1:363562/1‑251 (MQ=255)                      tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGCTTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGa                                                                                                                                                                                                  >  2:319516/1‑249 (MQ=255)                                cACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCa                                                                                                                                                                                       >  1:172842/1‑250 (MQ=255)                                 aCGAAACCAGTTTCCCCCAAGAGCTGCGGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                                                                                                     <  1:25639/251‑1 (MQ=255)                                 aCGAAAACAGCTTCCCGCAAGGGCTGGTGGATAAACTGTTCGATCGCGCTAATCTGCCGGTTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGTTCATTCTGGATCATCTGGGCGAAGATCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTTGATAAACTCTACGATCTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                                                                                                     <  1:102605/251‑1 (MQ=255)                                   gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                                                                                                   >  1:21994/1‑251 (MQ=255)                                   gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                                                                                                   >  1:224817/1‑251 (MQ=255)                                    aaaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAACCCGGCAAGCCGATTTTCTCTACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAAaggc                                                                                                                                                                                  >  1:50065/1‑251 (MQ=255)                                         aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                                                                                              >  1:388344/1‑250 (MQ=255)                                         aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                                                                                              >  1:46657/1‑250 (MQ=255)                                         aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                                                                                              >  1:345111/1‑250 (MQ=255)                                          gCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGCAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGCAGAGCTGACCAACGGCGGATAc                                                                                                                                                                            >  2:48938/1‑251 (MQ=255)                                           cTTCCCGCAAGAGCTGCTGGCTAAACTGGCCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTACGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCtcc                                                                                                                                                                            >  2:651832/1‑248 (MQ=255)                                           cTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAc                                                                                                                                                                            >  1:106666/1‑250 (MQ=255)                                              cccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACaacc                                                                                                                                                                        >  1:172989/1‑249 (MQ=255)                                                   aaggtttgctgGATAACCGGGCCGACCGCGCTAATGTGCCGG‑TGTCCCCCGTAACCTGCCACTCTTCCGTTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                                                                                                                   <  2:21994/245‑1 (MQ=255)                                                        ctaatagATAAACTGGTCGAACGAGCTAATCTGCCGGGTTACCTCGGTAAA‑TGCCACTCTTCCGGTACGGGCATTCTGGATCAGCTGGGCGCAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGGTCAAGATTGCCTGCCACGAAGGAACCTTGGGCCTGAATAAACTCTACGAACTGTGCGGAGTCGCCCGTGAAGAGCGGACCAACTGCTGCTaacccatcggtagggt                                                                                                                                                              >  2:160310/7‑236 (MQ=255)                                                           gtgGATAAACGGGCCGAACGCGTCAATCTGCCGGTTTGCCTCGGAAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TAt                                                                                                                                                          <  2:508139/250‑1 (MQ=255)                                                           ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcggc                                                                                                                                                                                  >  1:146440/1‑228 (MQ=255)                                                              gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCAAAACATCGGTCGGGT‑TATCGc                                                                                                                                                       >  2:294656/1‑251 (MQ=255)                                                                    cgggTCGAACGCGCTAATCGGCCGGGTTACCTCGGTAAC‑TGCAACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTcc                                                                                                                                                 <  1:508257/249‑1 (MQ=255)                                                                     tGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCg                                                                                                                                                >  1:564585/1‑251 (MQ=255)                                                                         caaccgcgcagttcttccGGGTTCCCCCTGTAAC‑GGCCCCTCTTCCGGGACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCGATACCTCCGCTGCCCCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTGTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCTTCCGTTTa                                                                                                                                            <  2:350011/235‑1 (MQ=255)                                                                                  tAATCTACCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAACCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCACACAATCGGTCGGGTGTA‑CGCCCATCCGTTTATCGGTGAc                                                                                                                                    >  2:587606/1‑250 (MQ=255)                                                                                   aaTCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATAGCTCGTCCGTTTATCGGTGACa                                                                                                                                   >  1:78078/1‑250 (MQ=255)                                                                                    atcggCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATGCGGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACaaa                                                                                                                                 <  1:607517/247‑1 (MQ=255)                                                                                            ggTTACCTCAGTAAC‑GGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTc                                                                                                                                                     <  1:27740/223‑1 (MQ=255)                                                                                                    cGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATa                                                                                                                                                                                                                                 <  1:127671/140‑1 (MQ=255)                                                                                                        aaccggCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTa                                                                                                              <  2:316872/246‑1 (MQ=255)                                                                                                         aC‑TGCCACTCTTCCGGTACGGTCATTCTGGAACAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTGCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCTGTAACTTCCAGCGTAcc                                                                                                            <  2:526310/251‑1 (MQ=255)                                                                                                             gCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCgg                                                                                                          <  1:319516/250‑1 (MQ=255)                                                                                                              cccgacgctcttccgaactaTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                        <  2:29512/232‑1 (MQ=255)                                                                                                               cACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                        >  1:301911/1‑250 (MQ=255)                                                                                                                         gTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAc                                                                                             <  1:26666/251‑1 (MQ=255)                                                                                                                           aCGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCt                                                                                           <  1:301163/251‑1 (MQ=255)                                                                                                                             ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAACTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACAACCGTCGTGC‑TATCGCTCGTCCGGTTATCGGTGACCAAGCCGGTAACATCCCGCGTACCGGTCACCGTCACGAACTgg                                                                                         >  2:251500/1‑251 (MQ=255)                                                                                                                                cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTg                                                                                      <  1:130413/251‑1 (MQ=255)                                                                                                                                 aTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                        >  1:29512/1‑232 (MQ=255)                                                                                                                                    cTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTgg                                                                                 >  1:94773/1‑250 (MQ=255)                                                                                                                                        aTCAGCTGGGCGAAGAGCACATGAAAAACGGCAAGCCGATTTTCTATACCTCCGCTGACTACGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCCGTCGGTT‑TATCGCTCGTCCCTTATTCGGTGACAAAACCGGTCACTTCCCACGTACCGGTAACCGtaccgaccgggaagtggcgacg                                                                              >  2:675156/1‑230 (MQ=255)                                                                                                                                               gggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAACCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCGG‑TTGAGCCGCCCGCac                                                                      >  1:465348/1‑251 (MQ=255)                                                                                                                                                    aaGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATAGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCaccgac‑c                                                                >  1:478383/1‑251 (MQ=255)                                                                                                                                                        gCACATGAAAACCGGCAAGACGATTTTCTATACCTCCGCTGACTCCGTGATCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGGGCGAAATCGCCCGAGAAGAGCTGACAAACGGCGGCTACAAACTCGGTCGGGG‑TAACGCACGAACGTTTAACGGGGACCAAGCAGGTAAATTCCAGCGTACCGGTAAACGTGACGACCTGGCTG‑GGGGGCCTCCAGCaccgac‑cgggc                                                            >  2:9619/1‑248 (MQ=255)                                                                                                                                                         cacaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGtgct                                                           >  1:539849/1‑251 (MQ=255)                                                                                                                                                           caTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTG‑TTGGGCCGCCAGCCACGAC‑CGtgctgg                                                         >  1:133902/1‑250 (MQ=255)                                                                                                                                                            aTGAAACCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCa                                                        <  1:57819/251‑1 (MQ=255)                                                                                                                                                               aaaaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGGGACAAAGACGGTAACTTCCAACGTACCGGGAACCGTCACGACCGGGCTG‑GTGGGCCCCCAGGacccgccctggtgcaga                                                       >  2:558897/1‑236 (MQ=255)                                                                                                                                                                  aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTTACAAAGCCGGTAACTTCCAGCGTACCGGGAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                  >  1:537280/1‑251 (MQ=255)                                                                                                                                                                  aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                  <  1:651832/251‑1 (MQ=255)                                                                                                                                                                  aCCGGAAAGCCTATTTTCGATCCCTCTGCTGATTCCGTGTTCCAGATTGCCTTCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTTCGAAATCGCCCTTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCGCGTCCGTTTACCGGTGACAAAGCCGGTAACTTCCAGCGTTCCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                  <  2:106666/251‑1 (MQ=255)                                                                                                                                                                   ccGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGGCACGCCCTGGCTG‑TTGGGCCGCCAGCACCGAA‑CGGGCTGCAGAAACTg                                                 >  1:557871/1‑251 (MQ=255)                                                                                                                                                                      gCAAGCCGATTTTCTATTCCTCCGCGGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGAAAAAGCCGTTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGtt                                              <  1:47343/251‑1 (MQ=255)                                                                                                                                                                            cGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTTTCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCCGCCCCCAC‑CGCGCTGCAGAAACTGGTtgatta                                         >  2:235200/1‑248 (MQ=255)                                                                                                                                                                               ttttCTATACCTCCGCTGACTCCGTGTTCCCGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGaaaa                                      >  1:520621/1‑250 (MQ=255)                                                                                                                                                                                 ttCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAAATGGTTGATGAAAAac                                    >  1:90776/1‑250 (MQ=255)                                                                                                                                                                                           ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTGCAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtg                         <  1:235200/251‑1 (MQ=255)                                                                                                                                                                                            cGCTGACTCCGTGTTCCAGTTTGCTGGACACGAAGGAACCGTTGGCCTGGATAAACACTACGAAATGTGCGAAATCGCCCGTGAAGAGCGGCCCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTGGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtgg                        <  2:172989/251‑1 (MQ=255)                                                                                                                                                                                              cgtactccgtttgtcAGATTGCCTGCCCAGAAGGAACCTTTGGCCTGGATAAACTCTACGTACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGtt                      <  2:564585/237‑1 (MQ=255)                                                                                                                                                                                               tGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGA‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGACACTTCCAGCGTACCGGGAAACGGCACGACCTGGCTGGTTG‑GCCGCCAGCCCCGACACGGCCTAC‑GAAACCTGTTGATGAAAAACACGGGCAggagggtt                     >  2:332042/1‑245 (MQ=255)                                                                                                                                                                                                  cTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAac                                    <  1:235364/233‑1 (MQ=255)                                                                                                                                                                                                  cTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGAGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAACTGGCTG‑TTGGGCCGCCAGCACCGAC‑CGGGCTGCAGAAACTGGGTGGTGTAAAacagaacggaagagaggcgg                  >  2:235364/1‑234 (MQ=255)                                                                                                                                                                                                    ccGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGTTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                <  1:587606/251‑1 (MQ=255)                                                                                                                                                                                                         ttCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGGGCCGCCAGCAACGAC‑CGTGCTGCAGAAACCGGTTGATGAAAAAAACGGTAAGGGGGGTTCTGTCGGTa            >  1:371230/1‑250 (MQ=255)                                                                                                                                                                                                           ccAGATGGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTTTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaaa          <  2:94773/250‑1 (MQ=255)                                                                                                                                                                                                             aGATTGCCTGCCACGAAGGAACCTTGGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACGTCCAGCGTACCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAtt       <  1:319198/251‑1 (MQ=255)                                                                                                                                                                                                                ttGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTtgatga                                         >  1:360305/1‑214 (MQ=255)                                                                                                                                                                                                                  gcctgccACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGcaga  <  1:626593/251‑4 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       GGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGA  >  NC_000913/4619929‑4620387

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