Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,852,890	C→T	*33* (TAG→TAA)	ivbL ←	ilvB operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,852,890	0	C	T	100.0%	46.7 / NA	16	*33* (TAG→TAA)	ivbL	ilvB operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/7);  total (9/7)

GCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATT  >  NC_000913/3852785‑3853022                                                                                                          |                                                                                                                                     gCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                                     >  2:588868/1‑139 (MQ=255)                                      gTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                                     >  1:108336/1‑102 (MQ=255)                                      gTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                                     <  2:108336/102‑1 (MQ=255)                                        gCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTagtagt                                                              >  2:802566/1‑139 (MQ=255)                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                             <  1:21229/139‑1 (MQ=255)                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                             >  2:296741/1‑139 (MQ=255)                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                             >  2:41485/1‑139 (MQ=255)                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                             >  2:95660/1‑139 (MQ=255)                                                    aaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        <  1:487676/133‑1 (MQ=255)                                                    aaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        >  2:487676/1‑133 (MQ=255)                                                    aaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGc                                                  >  2:313925/1‑139 (MQ=255)                                                             cgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagta                                         <  2:602536/139‑1 (MQ=255)                                                              ggTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTTGTTTTGCGTTGAGCATGGAagtag                                        <  1:296741/139‑1 (MQ=255)                                                                                              ttgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAAc        <  1:588868/139‑1 (MQ=255)                                                                                               tgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACa       <  2:90926/139‑1 (MQ=255)                                                                                                    cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATAtt  >  2:66734/1‑139 (MQ=255)                                                                                                          |                                                                                                                                     GCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATT  >  NC_000913/3852785‑3853022

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