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

TGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATA  >  NC_000913/3852769‑3853020                                                                                                                          |                                                                                                                                   tgCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcgac                                                                                                                   <  1:1144882/139‑1 (MQ=255)            gCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGAccaccacca                                                                                                        <  2:1082451/139‑1 (MQ=255)                       gTCCTTTTTTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCacgacg                                                                                             >  1:766463/1‑139 (MQ=255)                       gTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCacgacg                                                                                             >  2:1257439/1‑139 (MQ=255)                             ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACg                                                                                                                                <  2:687518/98‑1 (MQ=255)                             ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACg                                                                                                                                >  1:687518/1‑98 (MQ=255)                             ttcttcTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACt                                                                                       >  1:653171/1‑139 (MQ=255)                                 tctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                                     <  2:859237/105‑1 (MQ=255)                                 tctGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                                     >  1:859237/1‑105 (MQ=255)                                               ccGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                                     >  1:993843/1‑91 (MQ=255)                                               ccGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCcg                                                                                                                     <  2:993843/91‑1 (MQ=255)                                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccaccac                                                                                                          <  2:81956/92‑1 (MQ=255)                                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGaccaccac                                                                                                          >  1:81956/1‑92 (MQ=255)                                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACcaca                                                                                                      <  2:931970/96‑1 (MQ=255)                                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACcaca                                                                                                      >  1:931970/1‑96 (MQ=255)                                                         cTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGtt                                                           >  1:787010/1‑139 (MQ=255)                                                                     aCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                               <  2:787010/139‑1 (MQ=255)                                                                     aCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                               <  1:1256052/139‑1 (MQ=255)                                                                     aCCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCa                                               <  1:1257439/139‑1 (MQ=255)                                                                        ccgcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcg                                                                                    <  1:740522/99‑1 (MQ=255)                                                                        ccgcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcg                                                                                    >  2:740522/1‑99 (MQ=255)                                                                         cgcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcggcg                                                                                 <  1:360510/101‑1 (MQ=255)                                                                         cgcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTgcggcg                                                                                 >  2:360510/1‑101 (MQ=255)                                                                          gcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACgacgac                                                                                           >  1:401888/1‑90 (MQ=255)                                                                          gcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACgacgac                                                                                           <  2:401888/90‑1 (MQ=255)                                                                          gcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                          >  2:161302/1‑107 (MQ=255)                                                                          gcccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgcg                                                                          <  1:161302/107‑1 (MQ=255)                                                                           cccgGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAag                                         <  1:1281175/139‑1 (MQ=255)                                                                                tttGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagtag                                      >  2:1087734/1‑137 (MQ=255)                                                                                tttGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAagtag                                      <  1:1087734/137‑1 (MQ=255)                                                                                         gcggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGAc                           >  2:336088/1‑139 (MQ=255)                                                                                          cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                                   >  1:844896/1‑66 (MQ=255)                                                                                          cggGGTTTTGGAATCGTGTGTTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCAcacgc                                                                                                   <  2:844896/66‑1 (MQ=255)                                                                                                         gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgc                                                                           >  1:1057159/1‑75 (MQ=255)                                                                                                         gtgtgtTGTTCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGgc                                                                           <  2:1057159/75‑1 (MQ=255)                                                                                                              ttgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAAc      <  2:653171/139‑1 (MQ=255)                                                                                                               tgttCCAGTCCTTACGGCGCATTGCCGACGACCACCGCCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACa     <  2:625963/139‑1 (MQ=255)                                                                                                               tgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACa     <  1:336088/139‑1 (MQ=255)                                                                                                               tgttCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACa     <  1:1223765/139‑1 (MQ=255)                                                                                                                  tCCAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAata  <  2:1139201/139‑1 (MQ=255)                                                                                                                   ccAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCt                         >  2:359861/1‑115 (MQ=255)                                                                                                                   ccAGTCCTTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCt                         <  1:359861/115‑1 (MQ=255)                                                                                                                          |                                                                                                                                   TGCCCGAACTTGCCATGCTCCAGTCCTTTTCTTCTGGGCCGACTTTCCGGGAGGTGCTTAAACGAAAAACCCCGCCCGGTTTGCGCCGGCGGGGTTTTGGAATCGTGTGTTGTTCCAGTCCCTACGGCGCATTGCCGACGACCACCACCACACGCACGACGACCACTGCGGCGGATGGCGCAGTTGGTAGTAGTTTTGCGTTGAGCATGGAAGTAGTCATTGGGGACCTTGTTGGTTTTGTGTTTAACAATA  >  NC_000913/3852769‑3853020

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: ‑