Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,242,914	C→T	Q201Q (CAG→CAA)	malG ←	maltose ABC transporter membrane subunit MalG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,242,914	0	C	T	92.9%	32.1 / NA	14	Q201Q (CAG→CAA)	malG	maltose ABC transporter membrane subunit MalG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/3);  minor base G (0/1);  total (10/4)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CGGGTTGAGGTATTGCTGCATCCCCACGGCCAGGGTGTAGCTGTTTACGTCACGCAGTAACAGCGACGCGACCGGAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCGCCGTGAGTATTCAGGCCAATGAATGGAATGTACTCACCCAGACGATCAAACAACGCATACAACGCGAC  >  NC_000913/4242752‑4243105                                                                                                                                                                   |                                                                                                                                                                                                ngggTTGAGGTATTGCTGCATCCCCACGGCCAGGGTGTAGCTGTTTACGTCACGCAGTAACAGCGACGCGACCGGAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCTTGCCACg                                                                                                                                                                                          >  1:117817/2‑170 (MQ=255)                  gcATCCCCACGGCCAGGGTGTAGCTGTTTACGTCACGCAGTAACAGCGACGCGACCGGAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACCGCGGCAACAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAgc                                                                                                                                                                         >  1:254468/1‑170 (MQ=255)                              ccAGGGTGTAGCTGTTTACGTCACGCAGTAACAGCGACGCGACCGGAACTTCAGTAAGGGCAGCGATAAACGACAGGATGAATACCTCCGCCAGAATCGGGACTGACAGCGGCAACAGGACAAGGCGGAAGGCGTGCCACGGTGTCGCACCATCCAGCGCAGCAGcttcn                                                                                                                                                             <  1:40573/170‑2 (MQ=255)                                  ggTGTAGCTGTTTACGCCACGCAGTAACAGCGACGCGACCGGAACTTCAGTAATGGCAGCGATAAACGACAGGATGACTACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCTTGCCACGGTGTCg                                                                                                                                                                                    >  1:365240/1‑143 (MQ=255)                                                                      nacCGGAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAAAAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCACACCCAGATACGCGaaa                                                                               >  1:91749/2‑208 (MQ=255)                                                                          ggAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGTTTGCCACTGTGTCGCACCAGCCAGCGCAGCAGCTTCTTCCAGCGAACGGTCGATGGTTTCGAAATAGCCTGTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATAcgcg                                                                                  <  1:254886/201‑1 (MQ=255)                                                                           gAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAAt                                                                             >  1:34922/1‑205 (MQ=255)                                                                                          cAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTCCTGACAGCGGCAACAGGACACGGCGGAAGGCTTGCCCCGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTAc                                                                          >  1:206472/1‑193 (MQ=255)                                                                                                         ggATGAATACCACCGCCAGAATCGGTACTGACAGCGGCCACAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAgcgc                                                                                                     >  1:396877/1‑151 (MQ=255)                                                                                                                       gCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCACCGTGAGTATTCAGGCCAATGAATGGAATGTACTCACCCAGACGATCAAACAACGCATACAACGCGAc  >  1:100304/1‑236 (MQ=255)                                                                                                                                        cAGCGGCAACAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCGTCCAGCGAACTGTCGATGGTTTCTAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTAcgc                                                                        <  1:252285/149‑1 (MQ=255)                                                                                                                                             gCAACAGGACAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGATCCAGACATGCAGCGCAATCCCACCCAGATACGCGCAAATTACGCCGCCGTGAGTATTCa                                                         >  1:322712/1‑159 (MQ=255)                                                                                                                                                     aCAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCGCCGTGAGTATTCAGGCCAATgaatggaatg                                        >  1:139165/1‑168 (MQ=255)                                                                                                                                                      cAAGGCGGAAGGCTTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATcc                                                                                                <  1:174579/111‑1 (MQ=255)                                                                                                                                                                   |                                                                                                                                                                                                CGGGTTGAGGTATTGCTGCATCCCCACGGCCAGGGTGTAGCTGTTTACGTCACGCAGTAACAGCGACGCGACCGGAACTTCAGTAATGGCAGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCGCCGTGAGTATTCAGGCCAATGAATGGAATGTACTCACCCAGACGATCAAACAACGCATACAACGCGAC  >  NC_000913/4242752‑4243105

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