Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,804	T→G	L192L (CTT→CTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,804	0	T	G	92.6%	73.4 / ‑1.5	27	L192L (CTT→CTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (14/11);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

CTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATG  >  NC_000913/4618566‑4619051                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        ctGGGCGCGGCCTCGGTCATACTGGCGGTACGCTCGAAAAACTGGAATCCACCCCGGGCTTCGACATTTCCCCGGATGACACCCGTTTCCGCGAAATTATTAAAGCCGTCGGCGTGGCGATTATCGGTCAGCCCAGCTCACTGGCTCCGGCGGATAAACTTTTCTACGCGCCCCGTGATATTCCCGCAACCGTGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAg                                                                                                                                                                                                                                                  <  1:148799/246‑1 (MQ=255)                gTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGAGATATTAACGCCACCGTGGAC‑TCCATCCCGCTGCTCCCCGCCTCTATTCTGGCTGAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATg                                                                                                                                                                                                                              >  1:276909/1‑251 (MQ=255)                         gctgTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGTAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaa                                                                                                                                                                                                                      >  1:197078/4‑250 (MQ=255)                                    aaCAAACTGGAACCCATCCCTTGCTGCTATATTTTCCCGGATCACATCCTTTTCCGCAAAATTATTAAAGGCGTCTGCGTGGCTATTATCGTTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACTCGCCCCGTGATTTTCCCCCAACCG‑GGACCTCCACCCCGCTGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcgg                                                                                                                                                                                                          <  1:125547/250‑1 (MQ=255)                                           gggAATCCATCCCTGCTTTCAACATTTTCCCGGATGCCAACCGTTCCCGAGAAATTATTAAAGGCGCCGGCGGGGCGTTTATCGGTCAGACCAGCTCACTGGCGCCGGCGGATAACCTTTTCTACGCGACCTGTTATATTGCCGCACCCGTGGCC‑GCCTTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGttt                                                                                                                                                                                                   <  1:48242/250‑1 (MQ=255)                                                 ccatcaCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGGGAAAGTGGGTAGCGGCGCGTTTATGCCg                                                                                                                                                                                             >  1:304051/1‑251 (MQ=255)                                                     cccTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGACATTATTAAAGACCTCGGCCTGGCGATTATCGGTCACACCAGCCCCCTGGCTCCGGCGGATAAACGTTTCTCCGCGACCAGAGATATTACCGCAACCGTGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAACCTGGCGGAAGGTCTGGTAGCGCTGGGGATGGACGTGAACGAGGGTAGCGGGGCGTTTATGGCGAccc                                                                                                                                                                                         >  1:165713/1‑250 (MQ=255)                                                     cccTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGATCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATggcggcct                                                                                                                                                                                         >  1:271611/1‑244 (MQ=255)                                                            ttCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACAGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGGAct                                                                                                                                                                                  >  1:149201/1‑251 (MQ=255)                                                                                aacccGTTTCGGCGAAATTATAAAATACGTCGGCGTGGCAATTATCGGTAAGACCAGCTCACTGCCTCCGGCGGATAAACTTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAg                                                                                                                                                              <  1:249160/248‑1 (MQ=255)                                                                                                           cGTCGGCGTGGCTATTATCGTTCAGCCCACCTCACTGCCTCCGGCGGAAAAACGTTTCTACGCGCCCCGTGATATCACCGCACCCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTg                                                                                                                                   <  1:286872/251‑1 (MQ=255)                                                                                                                        aTTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGcaccac                                                                                                                      >  1:23199/1‑251 (MQ=255)                                                                                                                                  aGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTc                                                                                                            >  1:116673/1‑251 (MQ=255)                                                                                                                                      cagctCACTGGCTCCGGAGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTACCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCg                                                                                                        >  1:312900/5‑251 (MQ=255)                                                                                                                                                ccaccGGCGGATAAACGTTTCTACGGGACCCGTGATATTCCCGCACCCGTGGAC‑TCCACCCCGCTTATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCa                                                                                              <  1:88629/248‑1 (MQ=255)                                                                                                                                                      gcggATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGCAGCAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACt                                                                                        >  1:195598/4‑251 (MQ=255)                                                                                                                                                           tAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTggact                                                                                   >  1:337555/1‑248 (MQ=255)                                                                                                                                                                             ccGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAGCAGGGAATGGCCTCCAGTGCAGGtgaccggg                                                                 >  2:223886/1‑244 (MQ=255)                                                                                                                                                                                ggttatTACCGCACCCGTGGCC‑CCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTg                                                              <  1:299718/248‑1 (MQ=255)                                                                                                                                                                                           cAACCGTGGCC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGAcc                                                                                                                                                                                          <  1:337128/116‑1 (MQ=255)                                                                                                                                                                                           cAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGTGAACGTGGGTAGCGGCGCGTTTATGCCGAcc                                                                                                                                                                                          >  2:337128/1‑116 (MQ=255)                                                                                                                                                                                                       ccattccggtcaatcgccGCTTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATTCCGTCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGtt                                        <  2:61131/235‑1 (MQ=255)                                                                                                                                                                                                                     gaCCGCCTCTATTCGGGCGAAGAAACTGGCGGATGGTCGGGACGCGCTGGTGATGGTCGTGAAAGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTTCTGGCCACCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACAGGTGag                           <  2:33187/249‑2 (MQ=255)                                                                                                                                                                                                                           ctctATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGCCATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGt                     >  1:35575/1‑250 (MQ=255)                                                                                                                                                                                                                                        aagaaACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGtt       >  1:337035/1‑251 (MQ=255)                                                                                                                                                                                                                                         agaaACTGGCGGAAGGTCTGGTCGCGCTGGTGATGGACGTGAAAGTGGGTATCGGCGCGTTTATGCCGACCGACGAACTTTCTGAAGCCTGTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTTGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCTTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGttt      <  1:353364/251‑1 (MQ=255)                                                                                                                                                                                                                                             aCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGTCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATg  >  1:196819/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        CTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATG  >  NC_000913/4618566‑4619051

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