Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,513,468	T→C	S218P (TCC→CCC)	ydcT →	putative ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,513,468	0	T	C	100.0%	42.8 / NA	13	S218P (TCC→CCC)	ydcT	putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/8);  total (5/8)

AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAC  >  NC_000913/1513337‑1513600                                                                                                                                    |                                                                                                                                     aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  1:581015/139‑1 (MQ=255) aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  1:977363/139‑1 (MQ=255) aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  2:1086615/139‑1 (MQ=255) aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  2:209187/139‑1 (MQ=255)       aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtct                                                                                                                         <  1:55025/139‑1 (MQ=255)                                           ttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                     <  1:699955/139‑1 (MQ=255)                                                             aTCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGt                                                                   >  2:720852/1‑139 (MQ=255)                                                                                      cGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGaaa                                          >  1:1043732/1‑139 (MQ=255)                                                                                            tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgc                                                                                                                              <  1:454061/49‑1 (MQ=255)                                                                                            tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgc                                                                                                                              >  2:454061/1‑49 (MQ=255)                                                                                               ggCGGTTTTCAATAATGGACGCATCGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGc                                 <  1:720852/139‑1 (MQ=255)                                                                                                                gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGc                >  1:325531/1‑139 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAc  >  2:326399/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                     AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAC  >  NC_000913/1513337‑1513600

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