Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,393	C→G	A159G (GCT→GGT)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,393	0	C	G	100.0%	23.4 / NA	11	A159G (GCT→GGT)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (7/4);  total (7/4)

GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTTGTG  >  NC_000913/4008327‑4008432                                                                     |                                               gTTTCTCTGCCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGCTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                     <  2:1973290/76‑1 (MQ=255) gTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                     <  2:4634551/76‑1 (MQ=255)  tttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGCTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                    >  2:4349596/1‑75 (MQ=255)   ttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                   <  1:4796033/76‑3 (MQ=255)   ttCTCTACCGTAACCGGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                   <  2:4594680/76‑3 (MQ=255)                                       ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg   >  1:2069/1‑74 (MQ=255)                                       ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg   >  1:538722/1‑74 (MQ=255)                                        gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  1:3856558/1‑73 (MQ=255)                                        gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  2:2019842/1‑73 (MQ=255)                                        gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  2:3326156/1‑73 (MQ=255)                                        gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  2:3764795/1‑73 (MQ=255)                                                                     |                                               GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTTGTG  >  NC_000913/4008327‑4008432

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