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

TGTGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGT  >  NC_000913/3852767‑3853011                                                                                                                            |                                                                                                                          tgtgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            <  2:494261/139‑1 (MQ=255) tgtgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            <  2:487244/139‑1 (MQ=255) tgtgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            <  1:1310749/139‑1 (MQ=255)    gCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgacg                                                                                                         <  2:1744271/139‑1 (MQ=255)     cccGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCATTCTTTACGGCGCATTGCcgccga                                                                                                        >  1:1612515/1‑139 (MQ=255)                             ttttCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            <  2:275283/111‑1 (MQ=255)                             ttttCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            >  1:275283/1‑111 (MQ=255)                             ttttCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGGTTGCGCCGGCGGGGGTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCa                                                                                <  2:1612515/139‑1 (MQ=255)                              tttCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCAc                                                                               >  1:1648730/1‑139 (MQ=255)                                    ctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTACCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcggc                                                                         >  2:1816508/1‑139 (MQ=255)                                    ctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTACCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcggc                                                                         <  1:1816508/139‑1 (MQ=255)                                      gggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCacac                                                                                            <  1:964914/118‑1 (MQ=255)                                      gggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCacac                                                                                            >  2:964914/1‑118 (MQ=255)                                      gggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCa                                                                                <  2:70928/130‑1 (MQ=255)                                      gggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCa                                                                                >  1:70928/1‑130 (MQ=255)                                      gggCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGcggcgg                                                                       >  1:1277737/1‑139 (MQ=255)                                               ttCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACAGGCACGACGACCACTGCGGCGGATGGCGCAg                                                              <  2:1648730/139‑1 (MQ=255)                                               ttCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAg                                                              <  1:1169135/139‑1 (MQ=255)                                                tCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGt                                                             <  1:1580513/139‑1 (MQ=255)                                                 ccGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcg                                                                           <  1:169275/124‑1 (MQ=255)                                                 ccGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcg                                                                           >  2:169275/1‑124 (MQ=255)                                                        gTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGgtagta                                                     >  2:1636749/1‑139 (MQ=255)                                                                     aaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                 >  1:551667/1‑114 (MQ=255)                                                                     aaaCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                 <  2:551667/114‑1 (MQ=255)                                                                       aCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                      <  1:470869/139‑1 (MQ=255)                                                                             cccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                          >  1:53961/1‑81 (MQ=255)                                                                             cccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                          <  2:53961/81‑1 (MQ=255)                                                                             cccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCa                                                               <  1:14533/108‑1 (MQ=255)                                                                             cccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCa                                                               >  2:14533/1‑108 (MQ=255)                                                                                 gTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                             >  2:1779860/1‑122 (MQ=255)                                                                                 gTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                             <  1:1779860/122‑1 (MQ=255)                                                                                      cgcCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            <  2:188526/54‑1 (MQ=255)                                                                                      cgcCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                            >  1:188526/1‑54 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                    >  1:1264590/1‑56 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                    <  2:1264590/56‑1 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                      <  1:1598719/118‑1 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                      >  2:1598719/1‑118 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTggg                     <  1:1798274/135‑1 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTggg                     >  2:1798274/1‑135 (MQ=255)                                                                                            cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGAcc                 >  1:702777/1‑139 (MQ=255)                                                                                                     tggaacGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGtt        <  2:639320/134‑1 (MQ=255)                                                                                                          cGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                    <  2:1016153/42‑1 (MQ=255)                                                                                                          cGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccac                                                                                                    >  1:1016153/1‑42 (MQ=255)                                                                                                          cGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGgtag                                                       >  1:1816913/1‑87 (MQ=255)                                                                                                          cGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGgtag                                                       <  2:1816913/87‑1 (MQ=255)                                                                                                          cGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTtgtg   >  1:984691/1‑139 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCa                                                                                >  1:1703737/1‑61 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCa                                                                                <  2:1703737/61‑1 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                 >  1:1186801/1‑76 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                 >  1:658213/1‑76 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                 <  2:658213/76‑1 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                 <  2:1186801/76‑1 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATgg                                   >  1:1359466/1‑106 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATgg                                   <  2:1359466/106‑1 (MQ=255)                                                                                                           gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTCCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTgtgt  >  1:743058/1‑139 (MQ=255)                                                                                                                   ttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCttgttg           <  1:1615906/122‑1 (MQ=255)                                                                                                                   ttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCttgttg           >  2:1615906/1‑122 (MQ=255)                                                                                                                      cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                             >  1:1484334/1‑85 (MQ=255)                                                                                                                      cAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                             <  2:1484334/85‑1 (MQ=255)                                                                                                                        gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCa                                                               >  1:1456927/1‑65 (MQ=255)                                                                                                                        gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCa                                                               <  2:1456927/65‑1 (MQ=255)                                                                                                                        gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                             >  1:949724/1‑83 (MQ=255)                                                                                                                        gTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCg                                             <  2:949724/83‑1 (MQ=255)                                                                                                                            |                                                                                                                          TGTGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGT  >  NC_000913/3852767‑3853011

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