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	90.5%	44.7 / ‑1.4	22	I105I (ATT→ATC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (13/6);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

GGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGC  >  NC_000913/4619682‑4620138                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              gggCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCCACGCTGAAGTTATCGGCGCGTAAGCATGGGCGCAC‑GAAATGTACTCCGGTAAAGATAACCCGGCTGGTGACTGGGACATAGCCGGTGGCCCgtgt                                                                                                                                                                                                                >  2:333415/1‑248 (MQ=255)                  cATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATTCCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGGTTGAGTGGGGat                                                                                                                                                                                                <  2:65361/250‑1 (MQ=255)                         cAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GGAATGTCATCCGGTAAAGCTACCCCGACTGGTAGCTGGGAAATCGCCGGTGTCCCGGTTCTGATTGAGTGGGGGTATTTCTcc                                                                                                                                                                                        >  1:263422/1‑251 (MQ=255)                            aaGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTCACGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGCATTGACGGCAACGCTGAAGTTATCGGCGCGTACCCATGGGCGGACTCAAATGCCA‑CCCGTATAGATAACCCCACTGGGCACTGGGAAATCACCCGGGTCCCGGTCCTGTTTGTGTGGGGaatttttcccgc                                                                                                                                                                                      >  2:161927/1‑239 (MQ=255)                                                gCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGc                                                                                                                                                                  >  1:142387/1‑250 (MQ=255)                                                                  ggCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACt                                                                                                                                               <  1:195662/251‑1 (MQ=255)                                                                           aaTCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGGAcg                                                                                                                                      >  1:100362/1‑251 (MQ=255)                                                                             acTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgcg                                                                                                                                    <  2:176541/250‑1 (MQ=255)                                                                                   aaaTCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGGGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATc                                                                                                                              >  1:243274/1‑251 (MQ=255)                                                                                       cTGACCCGACTGGGGCTGGCGACAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGCAGTTATCGGCGCGTACGCATAGGCGCAC‑GAAATATAATCCGGTAAAGAAACCCCGTCTGGTCACTGGGAAATCGCCGGGATCCCGGTTCTTTTTTAGTGGGGATATTTTTCCGCTCCCGCCAACCACTTCCCGCCAGTGCTGCTGGACAAACTGGTCGAACGCGCTactccacc                                                                                                                          >  2:360013/1‑244 (MQ=255)                                                                                                             aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                                                                                                          <  1:140633/245‑1 (MQ=255)                                                                                                             aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGGCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGGGTGGGGATGTTTCTCCGATCCCGAAAACAGCTTCCCGCAAGAGCTGCTGGGTAAACTGGTCGCACGGGCTAATCTGCCGGGTTCCCTcaggaac                                                                                                          >  2:140633/1‑239 (MQ=255)                                                                                                                 acacGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGGGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCtt                                                                                                >  1:207455/1‑251 (MQ=255)                                                                                                                                 ggTTTCATCCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCt                                                                                <  1:195258/251‑1 (MQ=255)                                                                                                                                   tttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTAAGCATGGGCGCAC‑GAAATGTCATCCGGTAAAAATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTgg                                                                              >  1:277491/1‑251 (MQ=255)                                                                                                                                            ggcggGAAGGTACGGAAACGCTGAGGTTATCGGGTCGTACGCTTGGGCGCAC‑GAAATGTCAGCCGGTAAAGAACCCCCGTCTGGTCCCTGGGAAACCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAACCGGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTgg                                                                     <  2:142387/251‑1 (MQ=255)                                                                                                                                                 gAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGa                                                                 >  1:293933/1‑250 (MQ=255)                                                                                                                                                      gACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGca                                                           <  1:239312/251‑1 (MQ=255)                                                                                                                                                                     gTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGTTCAGCTGGGCGAAGAGCACATGAAAACCGGCa                                             >  1:56004/1‑250 (MQ=255)                                                                                                                                                                                       tGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTc                          >  1:165880/1‑251 (MQ=255)                                                                                                                                                                                           cgcAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACGGGGAAATCGCAGGCGTACCGGTTCTGTTTGAGAGGGGATATTTCTCCGATCACGAAAACAGATTCCCGCAAGAGCTGATGGATAAACTGGTCGAACGCGCTAATCTTCCACGTTACCTAGGTCACTGCCTCACTTCCTGAACGGTCCTTCTGGATCAGCTGGGCGCACATCACATGAACACCGGCCAGCCGATTTTCtataaccccgcc                      >  2:374797/1‑243 (MQ=255)                                                                                                                                                                                                                aaaGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGCGTGGGGATATTTATCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAATGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgc  >  2:304194/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              GGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCAC‑GAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGC  >  NC_000913/4619682‑4620138

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