Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,702	T→C	S158S (AGT→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,702	0	T	C	96.0%	54.2 / NA	25	S158S (AGT→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/13);  minor base A (1/0);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

TGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGC  >  NC_000913/4618474‑4618940                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   tgATAAACATTCCACTGGTGGCGTCGGCGATGTGCCGTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATaaa                                                                                                                                                                                                                               <  1:106495/250‑1 (MQ=255)   aTAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGt                                                                                                                                                                                                                            <  1:206594/251‑1 (MQ=255)                                            aTGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCAt                                                                                                                                                                                 <  1:268154/251‑1 (MQ=255)                                                   ggCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGTCCAGCTCACGGGCTcc                                                                                                                                                                                                                                         <  2:192585/190‑1 (MQ=255)                                                          ggTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑cc                                                                                                                                                                   <  2:170274/250‑1 (MQ=255)                                                          ggTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑cc                                                                                                                                                                   <  1:126544/250‑1 (MQ=255)                                                                                         atCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGCATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGCAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGATCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGGGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGACGAAGAAAcacgcgggagg                                                                                                                                   >  1:163296/1‑241 (MQ=255)                                                                                                     ggCCTCGGTCATACTGACGGTACGCTCGACAAACTGGAATCCCTCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCCCCGGCGGATAAACGTTTTTACGCTGACC‑GG‑GATATTAACCCAACCGGGGACCCCCTCCCGCTGATCAACCACC‑CTATTCTGGCGAAGAAACCGGCGGAAGGCCTGGAcgggcc                                                                                                                       >  2:93857/1‑247 (MQ=255)                                                                                                                  cTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGCCGTg                                                                                                          >  1:253679/1‑251 (MQ=255)                                                                                                                             cccGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCACCTCACTGGCCCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTCCCGCACCCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAg                                                                                              <  1:4776/249‑1 (MQ=255)                                                                                                                                     aCTGGAATCCATCCCTGGCATAGGCATTTTCACGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTTTACGC‑GAACCGT‑GAAATTACCGCAACCGGGGACTACATCACGCTGATCA‑CCGCCTCTAATCTGGCGAAGAAACTGGCGGAAGGTCTGGACGAGCTGGGGCTGGACGGG‑GAAGTGGGTAGCGGagcgg                                                                                      >  2:238092/1‑250 (MQ=255)                                                                                                                                          aaTCCATCCCTGGCTTCGACATTTTCCCGGATGACACCCGTTCCCGCGAAATTATTAAAGCCGTCGGCGTGGCGATTATCGGTCAGCCCAGCTCACTGGCTCCGGCGGAAAAACGTTTCTACGC‑GACCCGT‑GATATTCCCGCAACCGTGGACCCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATg                                                                                 <  1:194403/251‑1 (MQ=255)                                                                                                                                                    tGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTCAGC‑GAACCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGGTCA‑CCGCCTCTCTTCTGGCGAAGAAACTGGCGGGAGGTCTGGACGCGCTGGTGATGGACGGG‑AAAGTGGGAAGCGGCGCGTTTATGGCGCCCTACg                                                                       >  1:106300/1‑251 (MQ=255)                                                                                                                                                        ttCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGCTCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTGTATgacggccgccgaacc                                                                   >  2:224464/1‑237 (MQ=255)                                                                                                                                                           gactttttCCCGGATGACAACCTTTCCCGCGAAATTATTAAAGCGGGCGTCGTGGCGATTACCGGTAAGCCCACCTCACTGGCCCCGGCGGATAATCTTTTCTCCGC‑GACCCTT‑GATATTACCGCACCCGTGGACTCCATCCCGCTTATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAActctc                                                                <  1:235125/247‑1 (MQ=255)                                                                                                                                                                          gACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAAAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGAGGGTAGCGGCGCGTTTATGCCGACCTACGAACGCTCTGAAGCCCTTGCCgg                                                 >  1:200242/1‑250 (MQ=255)                                                                                                                                                                                 gTTTCCGCGAACTTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGCCCAGCTAACTGGCTCCGGAGGATAAACGTTTCTACGC‑GACCCGT‑GATATTAACGAAACCGTGGACTCAATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGAG‑AAAGGAGGTAGCGGCGCGTTTAGGCCGACCTACGAACTCTCTAGAGGCCTTGCCGATGCGAtt                                          >  2:41727/1‑251 (MQ=255)                                                                                                                                                                                   ttCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTAGCGTTCAGACCAGCCCACGGGCTCCGGCGGATAACCTTTTCTCCGC‑G‑CCCGTTGATATTACCGCAACCGTGGACCCCACCCCGCTGACCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAttg                                         <  1:133157/250‑1 (MQ=255)                                                                                                                                                                                           aaTTATTAAAGCCGTGGTGGTGGCGATTATCGGTCAGACCAGCTCGCTGGCTCCGGCGTATAAACTTTTCGACGC‑GCCCGGA‑GATATGACCCCCACCGTGGCCTCCATCCCTCTGACCG‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACGCTCTGAAGCCCTTGCCGAAGCGCTTGTTGGCGTgg                                <  2:56203/251‑1 (MQ=255)                                                                                                                                                                                           aaTTATTAAACACGTCGGCGTGGCGATTATCGGTCAGACCAGATCACTGGCTCCGGCGGATAAACGTTTCTAAGC‑GACCCGT‑GATATTCCCGCAACCGTGGACTCCATACCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCGGGACGCGCTGGGGATGGACGTGCAAAG‑GGGAAGCGGCGCGTTTATGCCGACCTACGAACTCTATGAAGCCCTTGCCGCAGCGATGGTTGGCGTgg                                >  2:104550/1‑251 (MQ=255)                                                                                                                                                                                               attaAAGACGGCGGCGTGGCGATTATCGGTCAGACCAGCCCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTaa                            <  1:51356/251‑1 (MQ=255)                                                                                                                                                                                                  aaaGCCGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATCCCCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACg                          <  1:191348/250‑1 (MQ=255)                                                                                                                                                                                                  aaaGACGTCGGCGTGGCGCTTATCGGTCAGACCCGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGGAGCGGCGCGTTTATGCCGACCTTCGAACTCTCTGAAGCCCTTGCCGAAGCGGTTGTTGGCGTGGCTAACgg                         >  1:208980/1‑251 (MQ=255)                                                                                                                                                                                                  aaaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGG‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGGCGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGGCCTACGAACTCTCTGCAGCCCTTGCCGCAGCGCTTGTTGGCGTGGCTAACgg                         >  2:29872/1‑251 (MQ=255)                                                                                                                                                                                                                         cGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAACCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCgcgc  >  1:265862/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   TGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGC  >  NC_000913/4618474‑4618940

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