Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,306,736	T→G	S325A (TCC→GCC)	oppF →	murein tripeptide ABC transporter/oligopeptide ABC transporter ATP binding subunit OppF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,306,736	0	T	G	90.0%	24.6 / ‑5.5	10	S325A (TCC→GCC)	oppF	murein tripeptide ABC transporter/oligopeptide ABC transporter ATP binding subunit OppF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (3/6);  minor base C (1/0);  total (4/6)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGAACTGGGGACCTATGATGAGGTCTACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCA  >  NC_000913/1306492‑1306926                                                                                                                                                                                                                                                     |                                                                                                                                                                                               ggTAGAACTGGGGACCTATGATGAGGTCTCCCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGttt                                                                                                                                                                                          <  1:46607/251‑1 (MQ=255)                                cacaATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCAGACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTGGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACa                                                                                                                                                           <  1:13551/251‑1 (MQ=255)                                                                          gTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCg                                                                                                                 >  1:1194/1‑251 (MQ=255)                                                                                       aTCCGGATCTGGAGAAGAACAAAACCATCCAGTTACGGGAAGGGTAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTcc                                                                                                     <  1:181913/250‑1 (MQ=255)                                                                                                                       ttACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                                                                    >  2:38644/1‑251 (MQ=255)                                                                                                                              aaGGGGAGTTACCGTCGCCGATCAACCCGCCTTACGGTTGTGTTTTCCGTACCCGTTGCCCGATTGACGGTCAAGAGTGCGCCACAACACGTCTTATTCTGGAGGGGACTTTCAGACACCCCGTTTATTGCCTGACACTCGATCTGCTTTAAAAACAATATCGTCTGACAGTTGTAAGCCCTTATTCACGCGAAATGCGATTATTTATTCTCGCCTGATGCGAAGGCAAAGTTTGGTATCTTTTTcgtgcc                                                             >  2:151859/1‑247 (MQ=255)                                                                                                                                  ggAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATa                                                          >  2:152901/1‑250 (MQ=255)                                                                                                                                   gAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAAc                                                        <  1:38644/251‑1 (MQ=255)                                                                                                                                                       cccGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATc                                    <  1:152901/251‑1 (MQ=255)                                                                                                                                                                                         ccgtttgccgttcagaagTGCCCCAAAACACGTCCGGTTCTGGAGTGGATTTTAAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCTGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCa  <  2:108359/235‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                               GGTAGAACTGGGGACCTATGATGAGGTCTACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCA  >  NC_000913/1306492‑1306926

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