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%	92.6 / NA	30	*33* (TAG→TAA)	ivbL	ilvB operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/15);  total (15/15)

CCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATTTA  >  NC_000913/3852771‑3853024                                                                                                                        |                                                                                                                                        cccGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgacga                                                                                                                      >  2:1566149/1‑139 (MQ=255) cccGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgacga                                                                                                                      >  2:396090/1‑139 (MQ=255)   cGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGacc                                                                                                                    >  2:546968/1‑139 (MQ=255)     aaCTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                                  <  1:546968/139‑1 (MQ=255)         tGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccacca‑cc                                                                                                             <  1:1688864/139‑1 (MQ=255)         tGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccacca‑cc                                                                                                             <  1:85450/139‑1 (MQ=255)            cATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCTCCA‑Ccaca                                                                                                          <  2:1575739/139‑1 (MQ=255)                           ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgac                                                                                                                        <  1:2335/111‑1 (MQ=255)                           ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgac                                                                                                                        >  2:2335/1‑111 (MQ=255)                             cttctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTgc                                                                                         >  2:1079173/1‑139 (MQ=255)                                   ggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCACATTGCCGACGACCACCA‑CCAcacgc                                                                                                       >  2:403668/1‑119 (MQ=255)                                   ggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCACATTGCCGACGACCACCA‑CCAcacgc                                                                                                       <  1:403668/119‑1 (MQ=255)                                                 gAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGgta                                                                     <  1:396090/139‑1 (MQ=255)                                                       cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACgacgac                                                                                               >  2:280617/1‑107 (MQ=255)                                                       cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACgacgac                                                                                               <  1:280617/107‑1 (MQ=255)                                                       cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                               >  2:171913/1‑139 (MQ=255)                                                            aCGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                          <  1:1079173/139‑1 (MQ=255)                                                              gAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGt                                                         >  2:1060356/1‑139 (MQ=255)                                                               aaaaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTgcggc                                                                                      >  2:1195975/1‑108 (MQ=255)                                                               aaaaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTgcggc                                                                                      <  1:1195975/108‑1 (MQ=255)                                                               aaaaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTg                                                       >  2:35231/1‑139 (MQ=255)                                                               aaaaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTg                                                       >  2:840578/1‑139 (MQ=255)                                                                    ccccGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCAt                                                  <  1:1566149/139‑1 (MQ=255)                                                                         cccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAag                                             <  1:840578/139‑1 (MQ=255)                                                                              tttGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagtagt                                         <  1:1060356/139‑1 (MQ=255)                                                                                  cgcCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagtag                                          >  2:1253588/1‑133 (MQ=255)                                                                                  cgcCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagtag                                          <  1:1253588/133‑1 (MQ=255)                                                                                   gcagtcggGGTCTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTgg                                   >  2:981989/6‑139 (MQ=255)                                                                                                            ttgttCCAGTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAAc          <  1:873042/139‑1 (MQ=255)                                                                                                                    gTCCTTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATTta  >  2:739893/1‑139 (MQ=255)                                                                                                                        |                                                                                                                                        CCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCA‑CCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATTTA  >  NC_000913/3852771‑3853024

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