Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,391	+C	coding (474/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,391	1	.	C	100.0%	20.1 / NA	10	K158N (AAG→AAC)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (2/8);  total (2/8)

GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTTGTG  >  NC_000913/4008317‑4008432                                                                              |                                                 gTCGACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAg                                                   <  2:5578868/76‑1 (MQ=255) gTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTTA‑AAAAGCCATGCGCTGGGAg                                                   >  1:4762621/1‑76 (MQ=255) gTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAg                                                   <  1:3531462/76‑1 (MQ=255)  tCATCGTCAGTTGCTCTACCTTACCCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgc                                                  <  1:3964125/76‑2 (MQ=255)  tCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgc                                                  <  1:4875979/76‑2 (MQ=255)   cAACGTCAGTTTCTCTACCTTAGCCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcg                                                 <  2:3148274/76‑3 (MQ=255)   cAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcg                                                 <  1:2302911/76‑3 (MQ=255)   cAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcg                                                 <  1:4163196/76‑3 (MQ=255)   cAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAgcg                                                 <  2:2371275/76‑3 (MQ=255)    aaCGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAAGCCATCCCCCGGAAg‑ggg                                               >  2:586117/1‑74 (MQ=255)            tttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                    <  2:4415745/76‑2 (MQ=255)            tttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                    >  2:2872272/1‑75 (MQ=255)             ttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                   <  1:5258628/76‑3 (MQ=255)             ttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                   <  1:1105603/76‑3 (MQ=255)              tctcTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                  <  2:2142419/76‑4 (MQ=255)              tctcTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                  <  2:4149562/76‑4 (MQ=255)                                                cggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt    <  1:3155074/76‑2 (MQ=255)                                                 ggcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg   >  2:2714591/1‑74 (MQ=255)                                                  gcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  <  1:4395452/76‑4 (MQ=255)                                                  gcg‑TTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  <  1:1469210/76‑4 (MQ=255)                                                                              |                                                 GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCG‑TTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTTGTG  >  NC_000913/4008317‑4008432

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