Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	35.7 / NA	11	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/2);  total (9/2)

GATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGC  >  NC_000913/1652140‑1652557                                                                                                                                                                                                |                                                                                                                                                                                                                                   gATGGGAGGACATAGATCTGAAAAATGGAACAATGAGGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCTGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACAGTAACGAGACTGAGCAAAGAGCATATCAATGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAATAAGGCACCTTTGTTGTTCAGCCTGAAgtg                                                                                                                                                                         <  2:174319/251‑1 (MQ=255)            aTAGATCTGACAAATGGAAAAATGATGCTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAAGCAAACGGGTACTTCACCTTATTAAACCAGCAAAAGACGCTCTCCCGAGTCAGAAGTCATTTACGAGAATGAGCCAAAAACATATCATTGATGTTAACCTCCAAGGGAATTGCAGAACAAAAAAACAAAAATGCACCATTGCTTTTCAACCTGAAGAGTCagagagagg                                                                                                                                                              >  2:252287/1‑250 (MQ=255)                                aaTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAAATGACATTAACGAGACTGAGCAAAGAGGATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCAACTTTGTTTTTTAACATGAAGTGTCAGAGAGAGAAAAAAATTATAGTGACCAtt                                                                                                                                          >  1:318327/1‑251 (MQ=255)                                                  tttAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGAAGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAAAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGGGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCaaaa                                                                                                                         >  1:234717/1‑248 (MQ=255)                                                        aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTTACTCAATAAGGc                                                                                                                     >  1:297563/1‑248 (MQ=255)                                                                         tACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGAAGACCTTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGAGTCAGCGAGGGTAAAAAATTATGGTGACCATTTTACCGTTGCCTCAATAAGGCAGCAGTGGGACGCAACGAt                                                                                                  >  1:267521/1‑250 (MQ=255)                                                                          aCCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaa                                                                                                <  1:357229/251‑1 (MQ=255)                                                                                               cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                           >  1:76692/1‑251 (MQ=255)                                                                                                            catCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                               >  1:151415/1‑250 (MQ=255)                                                                                                                          aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACtt                                                >  2:86914/1‑251 (MQ=255)                                                                                                                                                                        aaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc  >  1:33120/1‑251 (MQ=255)                                                                                                                                                                                                |                                                                                                                                                                                                                                   GATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGC  >  NC_000913/1652140‑1652557

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