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

TGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATTTA  >  NC_000913/3852769‑3853024                                                                                                                           |                                                                                                                                       tgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgac                                                                                                                       <  2:826758/139‑1 (MQ=255) tgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgac                                                                                                                       >  2:757544/1‑139 (MQ=255)    ccGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCgacgac                                                                                                                    <  1:78403/139‑1 (MQ=255)    ccGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCgacgac                                                                                                                    <  1:771347/139‑1 (MQ=255)    ccGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCgacgac                                                                                                                    <  2:232608/139‑1 (MQ=255)       aaCTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                                 <  1:281561/139‑1 (MQ=255)          ttGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccaccac                                                                                                              <  1:14233/139‑1 (MQ=255)                             ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCAc                                                                                            >  2:1171729/1‑139 (MQ=255)                             ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACt                                                                                           >  2:26386/1‑139 (MQ=255)                                  ctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcga                                                                                                                        >  1:1199337/1‑105 (MQ=255)                                  ctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcga                                                                                                                        <  2:1199337/105‑1 (MQ=255)                                     ggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAA‑CGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGTGGCGGAt                                                                                  >  2:711184/1‑139 (MQ=255)                                       ccGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATg                                                                                 <  1:26386/139‑1 (MQ=255)                                          aCTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                              >  1:12527/1‑139 (MQ=255)                                          aCTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                              >  2:706482/1‑139 (MQ=255)                                               ccGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGtt                                                                         <  2:130336/139‑1 (MQ=255)                                                cGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGtt                                                                         <  1:1171729/139‑1 (MQ=255)                                                 gggAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTgg                                                                       <  2:518410/139‑1 (MQ=255)                                                      gTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGgtagta                                                                  >  2:23684/1‑139 (MQ=255)                                                      gTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGgtagta                                                                  >  2:48255/1‑139 (MQ=255)                                                        gCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGTGCAGTTGGTagtagt                                                                >  2:1234857/1‑139 (MQ=255)                                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                               >  2:86504/1‑139 (MQ=255)                                                             aaCGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGc                                                           >  2:1014743/1‑139 (MQ=255)                                                              aCGAAAAACCCCGCCTGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                          <  1:1053892/139‑1 (MQ=255)                                                              aCGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                          <  2:1182376/139‑1 (MQ=255)                                                              aCGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                          <  1:23684/139‑1 (MQ=255)                                                              aCGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                                          <  2:325605/139‑1 (MQ=255)                                                                gAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGTGCAGTTGGTAGTAGTTTTGCGtt                                                        <  1:1234857/139‑1 (MQ=255)                                                                   aaaCCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAg                                                     <  1:732112/139‑1 (MQ=255)                                                                     aCCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTTCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                                   <  2:866171/139‑1 (MQ=255)                                                                     aCCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                                   <  1:757544/139‑1 (MQ=255)                                                                      ccccGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGc                                                                                                                           >  2:1136469/1‑66 (MQ=255)                                                                      ccccGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGc                                                                                                                           <  1:1136469/66‑1 (MQ=255)                                                                            ccgGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagt                                            <  1:706482/139‑1 (MQ=255)                                                                                tttGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTc                                        >  2:938878/1‑139 (MQ=255)                                                                                  tGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCAt                                      <  1:938878/139‑1 (MQ=255)                                                                                     gcCGGC‑GGGGTTTTGGAA‑CGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGTGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTggg                                  <  1:711184/139‑1 (MQ=255)                                                                                               gTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCttgt                          >  2:1091250/1‑139 (MQ=255)                                                                                                      aaTCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGtttt                   <  1:1091250/139‑1 (MQ=255)                                                                                                          gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                              >  2:606190/1‑76 (MQ=255)                                                                                                          gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                              <  1:606190/76‑1 (MQ=255)                                                                                                          gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTgtgt               >  2:9785/1‑139 (MQ=255)                                                                                                            gtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGttt             <  2:12527/139‑1 (MQ=255)                                                                                                               ttgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAAc          <  1:86504/139‑1 (MQ=255)                                                                                                                tgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACa         <  1:9785/139‑1 (MQ=255)                                                                                                                  ttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                                   >  2:17392/1‑95 (MQ=255)                                                                                                                  ttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                                   <  1:17392/95‑1 (MQ=255)                                                                                                                    ccAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAatat     >  2:385068/1‑139 (MQ=255)                                                                                                                       gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATTta  >  2:47301/1‑139 (MQ=255)                                                                                                                           |                                                                                                                                       TGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGC‑GGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATATTTA  >  NC_000913/3852769‑3853024

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