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

TGTGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATT  >  NC_000913/3852767‑3853022                                                                                                                            |                                                                                                                                     tgtgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                                       <  2:520376/139‑1 (MQ=255)   tgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgac                                                                                                                     >  2:770069/1‑139 (MQ=255)            ttGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccaccac                                                                                                            <  2:143818/139‑1 (MQ=255)                    cTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccacc                                                                                                              <  2:454589/129‑1 (MQ=255)                    cTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccacc                                                                                                              >  1:454589/1‑129 (MQ=255)                    cTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgca                                                                                                    >  2:695353/1‑139 (MQ=255)                        aGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCacgac                                                                                                <  2:459876/139‑1 (MQ=255)                               ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACt                                                                                         >  1:544728/1‑139 (MQ=255)                                       ggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGa                                                                                 <  2:544728/139‑1 (MQ=255)                                            aCTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                            >  2:579367/1‑139 (MQ=255)                                                 ccGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGtt                                                                       <  1:492683/139‑1 (MQ=255)                                                                aCGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        <  1:770069/139‑1 (MQ=255)                                                                   aaaaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTg                                                     >  2:812537/1‑139 (MQ=255)                                                                        ccccGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCAt                                                <  1:812537/139‑1 (MQ=255)                                                                           cgcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGa                                             <  1:579367/139‑1 (MQ=255)                                                                                 gTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagtagt                                       <  1:1060244/139‑1 (MQ=255)                                                                                  tttGCGCCGGCGTGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGaa                                            <  1:884065/133‑1 (MQ=255)                                                                                  tttGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGaa                                            >  2:884065/1‑133 (MQ=255)                                                                                           gcggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGAc                             <  1:695353/139‑1 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                               >  1:121395/1‑56 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                               <  2:121395/56‑1 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        >  2:34675/1‑111 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        <  1:34675/111‑1 (MQ=255)                                                                                                 ttttGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATg                                               >  2:309785/1‑115 (MQ=255)                                                                                                 ttttGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATg                                               <  1:309785/115‑1 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        >  2:195601/1‑96 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        >  2:905317/1‑96 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        <  1:905317/96‑1 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                        <  1:195601/96‑1 (MQ=255)                                                                                                             gtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGttt           >  2:1212625/1‑139 (MQ=255)                                                                                                              tgttgtTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTa          >  2:157647/1‑139 (MQ=255)                                                                                                                ttgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAAc        <  1:157647/139‑1 (MQ=255)                                                                                                                 tgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACa       <  1:144719/139‑1 (MQ=255)                                                                                                                   ttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTggg                                >  2:992715/1‑112 (MQ=255)                                                                                                                   ttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTggg                                <  1:992715/112‑1 (MQ=255)                                                                                                                      cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                            <  1:511221/65‑1 (MQ=255)                                                                                                                      cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                            >  2:511221/1‑65 (MQ=255)                                                                                                                      cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTTGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATAtt  <  2:473720/139‑1 (MQ=255)                                                                                                                      cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATAtt  >  2:349096/1‑139 (MQ=255)                                                                                                                        gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTgggg                               >  2:653329/1‑108 (MQ=255)                                                                                                                        gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTgggg                               <  1:653329/108‑1 (MQ=255)                                                                                                                            |                                                                                                                                     TGTGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATT  >  NC_000913/3852767‑3853022

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