Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,399,875	A→G	E208E (GAA→GAG)	miaA →	tRNA dimethylallyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,399,875	0	A	G	100.0%	159.1 / NA	45	E208E (GAA→GAG)	miaA	tRNA dimethylallyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/21);  total (24/21)

ACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCTTTT  >  NC_000913/4399798‑4399952                                                                              |                                                                              aCGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCt                                                                             <  2:838257/80‑1 (MQ=255)  cGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGc                                                                           <  2:838258/81‑1 (MQ=255)   ggAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAgctgct                                                                          >  1:838306/1‑81 (MQ=255)     aaCTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTTGCCATCGCCCCGGCGAGCCGTGAGCTGCTc                                                                         >  2:838304/1‑80 (MQ=255)       cTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCAt                                                                      <  1:838260/81‑1 (MQ=255)         gACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATc                                                                     >  1:838312/1‑80 (MQ=255)             cAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACg                                                                 >  1:838293/1‑80 (MQ=255)              aaaCGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCa                                                               >  1:838323/1‑81 (MQ=255)              aaaCGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCa                                                               >  2:838315/1‑81 (MQ=255)                  gTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCAtt                                                             >  2:838328/1‑79 (MQ=255)                    cAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGag                                                          >  2:838294/1‑80 (MQ=255)                        agaCGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCg                                                     >  1:838302/1‑81 (MQ=255)                             ctctACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTtc                                                <  1:838261/81‑1 (MQ=255)                               ctACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTtcat                                              >  2:838321/1‑81 (MQ=255)                                  ccGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTtcatc                                             <  1:838262/79‑1 (MQ=255)                                   cGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTtcat                                              >  2:838317/1‑77 (MQ=255)                                   cGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGa                                          <  2:838262/81‑1 (MQ=255)                                    gTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGAt                                         >  2:838306/1‑81 (MQ=255)                                    gTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGAt                                         >  1:838331/1‑81 (MQ=255)                                    gTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGAt                                         <  2:838263/81‑1 (MQ=255)                                      aTCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGt                                       >  2:838308/1‑81 (MQ=255)                                       tCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGtt                                      <  2:838264/81‑1 (MQ=255)                                         aGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTgg                                    >  2:838312/1‑81 (MQ=255)                                              cATCAGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTc                                <  1:838266/80‑1 (MQ=255)                                                  aGTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGtt                           >  2:838318/1‑81 (MQ=255)                                                   gTTCGCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGttt                          <  2:838266/81‑1 (MQ=255)                                                       gCCATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGa                       <  1:838268/80‑1 (MQ=255)                                                         cATCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGc                    <  1:838269/81‑1 (MQ=255)                                                           tCGCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAg                  <  2:838269/81‑1 (MQ=255)                                                             gCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAg                  >  2:838336/1‑79 (MQ=255)                                                             gCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAg                  >  1:838317/1‑79 (MQ=255)                                                             gCCCCGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGaa                <  2:838270/81‑1 (MQ=255)                                                              ccccGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAg               <  2:838271/81‑1 (MQ=255)                                                               cccGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGt              >  2:838319/1‑81 (MQ=255)                                                                ccGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTc             >  2:838334/1‑81 (MQ=255)                                                                 cGGCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTcc            <  1:838273/81‑1 (MQ=255)                                                                  ggCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCg           <  2:838273/81‑1 (MQ=255)                                                                   gCGAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCgg          <  1:838275/81‑1 (MQ=255)                                                                     gAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGgc        <  2:838275/81‑1 (MQ=255)                                                                     gAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGgc        >  1:838333/1‑81 (MQ=255)                                                                     gAGCCGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGgc        >  1:838342/1‑81 (MQ=255)                                                                        ccGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCt     <  1:838277/81‑1 (MQ=255)                                                                        ccGTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCt     >  2:838327/1‑81 (MQ=255)                                                                          gTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCttt   >  2:838320/1‑81 (MQ=255)                                                                          gTGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCttt   >  2:838325/1‑81 (MQ=255)                                                                           tGAGCTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCtttt  <  1:838278/81‑1 (MQ=255)                                                                              |                                                                              ACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTGGCTTCAGGTTTTGAAGCAGAAGTCCGGGCGCTTTT  >  NC_000913/4399798‑4399952

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