Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,454	A→G	N409S (AAC→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,454	0	A	G	94.7%	47.4 / ‑6.1	19	N409S (AAC→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (8/10);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGG  >  NC_000913/4619239‑4619684                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       gATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCCTTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAATGCGAAAACAGCCGGCAGGAAGGGGCGaaaaggggggaagcggggc                                                                                                                                                                                                     >  2:135155/1‑235 (MQ=255)         cGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGGCGAAAACAGCTGGCAGGAAGCGGCGAAAACGGTGAAAGCGGGAATTaaaatt                                                                                                                                                                                             >  2:70286/1‑248 (MQ=255)                         tGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAGAAAAAAGCTGGCAGGAAGCGGGGAAAGCGGTGGAAGCGGCAATTAAACTTGCCGATAAAGAACCgg                                                                                                                                                                             >  1:177147/1‑251 (MQ=255)                                               gAAATGGATTCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCg                                                                                                                                                       <  2:83814/251‑1 (MQ=255)                                                        ccccGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAg                                                                                                                                              <  1:1020/250‑1 (MQ=255)                                                             cgggcTGGGGATGGCAGTGGTTGCAATGGGCAGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGTCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGATGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCTGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGGCTAGCGCCGTATCAGCGAAt                                                                                                                                         <  2:51015/248‑1 (MQ=255)                                                              gcgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa                                                                                                                                        <  1:132595/251‑1 (MQ=255)                                                                  tGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATCTGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGAGAAAGACGAAAACAGCTGGCAGGGAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACgg                                                                                                                                    >  1:157100/1‑251 (MQ=255)                                                                   ggggATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACgg                                                                                                                                    <  1:182724/250‑1 (MQ=255)                                                                                   gCAATGGGCGGCTGTCGTCGTCAGGCATCTGACACCATTGAGTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCa                                                                                                                    <  2:146923/250‑1 (MQ=255)                                                                                                     cGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAACGGCAATTAAAATTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCCTTTAAATTTGAAGCACTgg                                                                                                 >  1:143201/1‑250 (MQ=255)                                                                                                       tCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGGAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAACACCCGAAAGCACACCAACTGTCTATCGCCGTATCAACCAATAAAGGTATACTGATCTGATCATTTAAATTTTAAGCACTGAgg                                                                                               >  1:24114/1‑250 (MQ=255)                                                                                                                  acacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGTTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACat                                                                                    <  1:16774/251‑1 (MQ=255)                                                                                                                                         gttACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGTGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACAAACTGTCTATCGCCGTATCAGCGAATAACGGTATACCGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAACCGTGCATTTATTATGGt                                                             <  1:60269/250‑1 (MQ=255)                                                                                                                                                      cgcgTATGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACCCGAAAGACGAAAACCACTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGGCAATTAAAATTGCCGATAAAGCACCGGAAAGCACACCAACCGTCTCTCGCCGTATCAGCGAAAAACGGTAAACTGATCTGATCCTTTAAATTTTAAGCACTGAGGACGGGGGACATATGAAACGGGCATTTCTTATGGGGCTGGGATCATTc                                                >  2:115239/1‑251 (MQ=255)                                                                                                                                                                 gACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGgcgc                                     >  1:110301/1‑251 (MQ=255)                                                                                                                                                                                  gCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAAc                    <  2:177147/251‑1 (MQ=255)                                                                                                                                                                                          tGGCGGTTATCCACGCGAAAGACGCAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACAACAACTGTCTATCGCCGTATCAGCGAATAACGGTATAATGATCTGATCATTTAAATTTGAAGCACCGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCCTTCGGCATCGGCGCTACAGCAAATGCAGAACGCTTTggg            >  2:534/1‑250 (MQ=255)                                                                                                                                                                                                    ccACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCgggg  <  2:69913/251‑1 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       GATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGG  >  NC_000913/4619239‑4619684

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