Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,278,574	G→A	G32G (GGG→GGA)	ghxP →	guanine/hypoxanthine transporter GhxP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,278,574	0	G	A	81.8%	19.8 / ‑4.8	11	G32G (GGG→GGA)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (2/7);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTCACTCGAC‑GCCTGGTTTAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGG‑CTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTG  >  NC_000913/4278368‑4278815                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   tGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTCACTCGAC‑GCCTGGTTTAAAATTTCACAACGTTGAAGCACTGTCCGGCAGGAAGTGGTTGCCGGATTAACAACGTTTCTGGCGATgg                                                                                                                                                                                                                              <  1:72265/229‑1 (MQ=255)             aaTTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTCACTCGAC‑GCCTGGTTTAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGATTAACAACGTTTCTGGCTATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTggg                                                                                                                                                                                            <  1:388343/251‑1 (MQ=255)                            ggAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTCACTCGAC‑GCCTGGTTTAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGATTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGG‑CTTcc                                                                                                                                                                             <  1:96501/250‑1 (MQ=255)                                                 ncTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAATTTTCAGGGGAGGTTTTCTT‑TGCTTTA‑GCCCTCAGCCGCGTACCGGAGGTTC‑CTTCACACCCTGGGTTGAAATTTTCCAACCGGGAAGAACTTTCCGCCAGGAAGAGGGTGCCGGGTTAACACCGTTTTTGGGCAGGGGTTCCTCGGGCCTCTTCGTCCCCGGAATGTTGGGTAAA‑CAGGGCTTCCCCCCTGCGGACGTTTTCGTTGc                                                                                                                                                       >  1:85836/2‑251 (MQ=255)                                                                     cACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTAACTCGAA‑GCCTGGTTTAAAATTTTCCAAAGTGGAAGCACTGTCCCGCACGAAGTGGTTGCCGGGATaacaac                                                                                                                                                                                                                                            >  1:390732/1‑147 (MQ=255)                                                                               gACGAAAATTAAACTCTCAGGGGCTGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTCACTCGAC‑GCCTGGTTTAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGGTGCCGGATTAACAACGTTTCTTGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGG‑CTTCCCGCCTGCGGCAGTTTTCGGTGCAACCTGTCTGGTTGCCGGACTCGGTTCTAt                                                                                                                         >  1:5562/1‑251 (MQ=255)                                                                                                                                                                    aaaTTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGATTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGG‑CTTCCCGCCTGCGGCAGTTTTCGttn                                                                                                                                                        <  1:132382/138‑2 (MQ=255)                                                                                                                                                                          cacaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGATTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGG‑CTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTAGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTgggg                                 >  1:98664/1‑251 (MQ=255)                                                                                                                                                                                         cTGTCCGTCAGGAAGTGGTTGCCGTATTAACACCGTTTCTGGCGATGGTCTACTCGGTCGCCGTCGTTCCAGTTATGTTGGGTAAAGCAGG‑CTTCCCGCCTGCGGCAGTTTTCGTTGCAACCGGTCTGGTTGCCTGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTggn                                                                                 <  1:55731/188‑2 (MQ=255)                                                                                                                                                                                                         gggtGCCGGATTAACCACGTTGCGGGCGATGGTCAACTGGGTAATGGGCGTCCCAGGGATGTTGGGTAAAGCAGT‑CTGCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCGGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCTCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTg  <  1:237611/248‑1 (MQ=255)                                                                                                                                                                                                          gTTGCCGGAGTAATAACGTTTCTGCCTATGGTCGACTCGGGCATCTGCGTTCCAGGTATGTCGGGTAAAGCAGG‑CTCGCCGCCTGCGGCAGTTTTCGTTGCAACCTTTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGtn                                                                                                              <  1:39248/142‑2 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   TGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATG‑TCTACGCCATCAGC‑GCGTACCGGCGGTTCACTCGAC‑GCCTGGTTTAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGG‑CTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTG  >  NC_000913/4278368‑4278815

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