Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,862,357	G→A	P39L (CCG→CTG)	glvB ←	putative PTS enzyme II component GlvB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,357	0	G	A	90.0%	20.3 / ‑3.7	12	P39L (CCG→CTG)	glvB	putative PTS enzyme II component GlvB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (2/7);  total (3/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

AAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACCCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC  >  NC_000913/3862137‑3862584                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    aaaGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACCCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCGTTGTATTCGGCTTTTGAGTAGAGTTTCACTGCCGCATCTGCACGTCCCAGTGTGCACATGTTGAACTGCAGAATCAgtg                                                                                                                                                                                                       <  1:114463/251‑1 (MQ=255)                                                                                                                 ccGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTCTTAGGCTTTTGAGTAGAGTTTCACTTCAGCATCTTCACGTCCCAGCGTGCCCATGTTGAACTGAAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTCAAGCACAAc                                                                                                                                                                 >  1:224374/1‑177 (MQ=255)                                                                                                                 ccGGCAGCCTGATCCAGCTCTTTTTGTGGCTCGGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCCNCTCNGCATNTTNACGTCCCAGCGTGCANATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGNNGNNNANAGCACAAc                                                                                                                                                                 <  2:224374/177‑1 (MQ=255)                                                                                                                       nccTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCCGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAAc                                                                                                                                                                                             >  1:67031/2‑143 (MQ=255)                                                                                                                             gtgtgctcttccgttcgCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCAGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATAGCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGAcc                                                                           <  1:228479/235‑1 (MQ=255)                                                                                                                              ccAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTGCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGGAAAGCACAACCCAATTGCGATTTGGGTCAGCATCATATCCGCTTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGAcc                                                                           <  1:165975/250‑1 (MQ=255)                                                                                                                                                     gcgGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTGTCACTTCCGCATCTTCACGTCCCAGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGTTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCcaccacacn                                                   <  1:30363/251‑2 (MQ=255)                                                                                                                                                               ggCCTCGCGAGGCTTTGTATTCGGCGTTTGAGTAGAGTTTCACTGCCGCATCTTCACTTCCCTGCTTGCTCATGTTGAACTGCAGAATCATTGTGCGGAAAACCGCTAAGTACAGCACGGTAAAGCACAACCCAATGGCGATTTGGGTCAGcatca                                                                                                                                        <  1:286506/156‑1 (MQ=255)                                                                                                                                                                gCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCACAGCGTGGACATGTTGAACTACAGCATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACCACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACAAAAGAGATAc                                         >  1:228994/1‑250 (MQ=255)                                                                                                                                                                             ttGTATTCGGCTTTTGAGTAGAGTTTCGCTTCCGCATCTTCACGTCCCAGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGan                             <  1:261537/249‑2 (MQ=255)                                                                                                                                                                               gTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCAGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTTAACATCGGGATCCATTTTTTCGATAAAACCTAGGCAATCAGACCTCCGCCCATTTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTg                         <  1:48705/251‑1 (MQ=255)                                                                                                                                                                                                       cACTTCCGCATCTTCACGTCCCAGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGAAAATCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCGTAACATCGGTATCCAGTTTTTCGGTAAACTCTAGTCAATCATACCTCCGCCCGGGTTGCCCTCCACACCTAAGAGATACATTACAGTCGACATTGAGGCCGCCAGCACGGCGTGTAcc  <  1:307719/250‑1 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    AAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACCCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACC  >  NC_000913/3862137‑3862584

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