Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,521,815	A→G	V69A (GTT→GCT)	yjhF ←	putative transporter YjhF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,521,815	0	A	G	90.0%	21.2 / NA	10	V69A (GTT→GCT)	yjhF	putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (5/4);  minor base T (0/1);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGAC  >  NC_000913/4521575‑4522047                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         aCCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAATTTTATTTTTCCCGAAATAATTCATCATGGTCAGCGCTATTCGCTGCGAACCGCCCGCATCGGCCAGCAAGCGACCGAGCATAGCACCGAAGcc                                                                                                                                                                                                                                <  1:286978/251‑1 (MQ=255)     gCAAACCAGAAAGCATCGGCACACACATAAACAGCGGTGATATTTTTGCTTCACGGGCAATAGCTAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCTCATCGGCCAGCAAGCGACCGAGCATATCACCGAAGCCAAat                                                                                                                                                                                                                            <  1:291251/251‑1 (MQ=255)                                                  ttGCTTCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGAc                                                                                                                                                                                <  1:250865/250‑1 (MQ=255)                                                                                      gaagaaTTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCGAGTCAAGTTTATTTTTCCAGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCAACCGCATCGGCCAGCAAGAGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACtggcg                                                                                                                                           >  1:183406/1‑251 (MQ=255)                                                                                                   gACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATACTTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGAGACAGCa                                                                                                                               >  1:281435/1‑250 (MQ=255)                                                                                                         gAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCTGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACtggcgatggcgan                                                                                                                                   <  1:82808/240‑2 (MQ=255)                                                                                                                                                       cAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCTGCGGCCAGCAAGCGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATTGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGAGACGACTTTACTCAAATCCATACCACTGGCGATGGCGAc                                                                                                                                   <  1:322276/194‑1 (MQ=255)                                                                                                                                                                                                                               aGCAAGCGACCGAGCATAGAACCGAAGCCAAATATCAAACCAATATGGCCAAGCCTACCGCCAAGTCAGGATTCGACTGAGGTGACGACTTTACTCAACTCCATACCACTGGAGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAACGAAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGAc  >  1:225376/1‑251 (MQ=255)                                                                                                                                                                                                                                gCAAGAGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACaata                                                                                                                   >  1:267383/1‑137 (MQ=255)                                                                                                                                                                                                                                   aGCGACCGAGCATAGCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATNATTNNNNnaann                                                                                                       >  2:198506/1‑144 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         ACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGAC  >  NC_000913/4521575‑4522047

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