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%	24.1 / NA	11	K158N (AAG→AAC)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (7/4);  total (7/4)

GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑A‑AAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTTGTG  >  NC_000913/4008317‑4008432                                                                              |                                                 gTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAg                                                   <  1:3090775/76‑1 (MQ=255)  tCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTC‑ACAAGGCCATGCGCCGGGAg‑g                                                 >  2:4793629/1‑76 (MQ=255)  tCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAgc                                                  <  1:2762800/76‑2 (MQ=255)  tCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAgc                                                  <  2:1270364/76‑2 (MQ=255)  tCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAgc                                                  <  1:3251439/76‑2 (MQ=255)   cAACGTCAGTTTCTCTACCTTCACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAA‑AAGGCCATGCGCTGGAAgcg                                                 >  1:584667/1‑74 (MQ=255)   cAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAgcg                                                 >  1:582257/1‑74 (MQ=255)           gTTTCTCTGCCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGCTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                     <  2:1973290/76‑1 (MQ=255)           gTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                     <  2:4634551/76‑1 (MQ=255)            tttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAgcgg                                                <  1:3729064/68‑4 (MQ=255)            tttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGCTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                    >  2:4349596/1‑75 (MQ=255)             ttCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                   <  1:4796033/76‑3 (MQ=255)             ttCTCTACCGTAACCGGTCTGGAAGTGATGGCGGGCGGCGTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                   <  2:4594680/76‑3 (MQ=255)                                                 ggcgTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg   >  1:538722/1‑74 (MQ=255)                                                 ggcgTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg   >  1:2069/1‑74 (MQ=255)                                                  gcgTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  1:3856558/1‑73 (MQ=255)                                                  gcgTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  2:2019842/1‑73 (MQ=255)                                                  gcgTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  2:3326156/1‑73 (MQ=255)                                                  gcgTTTCAA‑AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg  >  2:3764795/1‑73 (MQ=255)                                                                              |                                                 GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑A‑AAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTTGTG  >  NC_000913/4008317‑4008432

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