Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,632,970	A→G	K421R (AAA→AGA)	livM →	leucine/isoleucine/valine transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,632,970	0	A	G	100.0%	157.0 / NA	44	K421R (AAA→AGA)	livM	leucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/22);  total (22/22)

GATCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCCTGC  >  minE/2632902‑2633037                                                                     |                                                                    gatCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:918530/1‑71 (MQ=255)  atCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:2592053/1‑70 (MQ=255)  atCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:135491/1‑70 (MQ=255)  atCTGGCGGCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   <  1:3051026/70‑1 (MQ=255)    cTGGCGTCCGCAGGGCTTGCTGCCCATGCCGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:291827/1‑68 (MQ=255)    cTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:1095325/1‑68 (MQ=255)    cTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:2197963/1‑68 (MQ=255)    cTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                   >  1:2127507/1‑68 (MQ=255)      ggCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCgagaggagag                                                              >  1:2986231/1‑71 (MQ=255)       gCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCgagaggagag                                                              <  1:2559182/70‑1 (MQ=255)       gCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCgagaggagag                                                              >  1:237023/1‑70 (MQ=255)       gCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCgagaggagag                                                              <  1:1065737/70‑1 (MQ=255)       gCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGgaga                                                               <  1:2720562/69‑1 (MQ=255)            cGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGc                                                         <  1:2817877/70‑1 (MQ=255)             gCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCAt                                                       <  1:739498/71‑1 (MQ=255)              cAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCAt                                                       >  1:775934/1‑70 (MQ=255)              cAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCAt                                                       <  1:1593365/70‑1 (MQ=255)                 tgCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCATCGagag                                                                   <  1:2218216/54‑1 (MQ=255)                        tgcCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCat                                            >  1:273245/1‑71 (MQ=255)                        tgcCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCat                                            >  1:2371822/1‑71 (MQ=255)                           ccATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCattat                                         >  1:2724080/1‑71 (MQ=255)                            cATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGc                                                             <  1:552034/50‑1 (MQ=255)                            cATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCattatt                                        <  1:1612714/71‑1 (MQ=255)                            cATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCattat                                         >  1:2868565/1‑70 (MQ=255)                             aTGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCa                                                            <  1:804877/50‑1 (MQ=255)                                aCGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCt                                    >  1:1346475/1‑71 (MQ=255)                                 cgcgCCCGGAACTGTAGCTGAAAAAGGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTg                                   <  1:2762355/71‑1 (MQ=255)                                   cgcCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGtt                                 <  1:907758/71‑1 (MQ=255)                                    gcCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCa                                             <  1:327419/58‑1 (MQ=255)                                     cccGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTaa                               <  1:1045612/71‑1 (MQ=255)                                      ccGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTaa                               >  1:30038/1‑70 (MQ=255)                                       cGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCatt                                           <  1:2753752/57‑1 (MQ=255)                                         cAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGc                           <  1:242833/71‑1 (MQ=255)                                         cAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGc                           >  1:2084243/1‑71 (MQ=255)                                         cAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGc                           <  1:2727950/71‑1 (MQ=255)                                          aaCTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGcc                          >  1:2098411/1‑71 (MQ=255)                                            cTGAAGCTGAAAAACGGCGCAGCGAGAGAAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCtg                        <  1:1207200/71‑1 (MQ=255)                                               aaGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCtgat                      >  1:915026/1‑70 (MQ=255)                                                  cTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGatgatg                  <  1:1268966/71‑1 (MQ=255)                                                  cTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGatgatg                  <  1:71159/71‑1 (MQ=255)                                                   tGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGtt                                 <  1:2415506/55‑1 (MQ=255)                                                    gAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATgcg                >  1:2807673/1‑71 (MQ=255)                                                     aaaaaCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCt                         <  1:1891979/61‑1 (MQ=255)                                                      aaaaCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATtattat                                      >  1:254684/1‑47 (MQ=255)                                                      aaaaCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCtgatga                    >  1:660913/1‑65 (MQ=255)                                                       aaaCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCtt             >  1:1544114/1‑71 (MQ=255)                                                       aaaCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCtt             <  1:1131175/71‑1 (MQ=255)                                                         aCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTaa                               >  1:1651321/1‑51 (MQ=255)                                                         aCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTcg           <  1:1182019/71‑1 (MQ=255)                                                              gcAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCggcggc      >  1:2442445/1‑71 (MQ=255)                                                               cAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGcc     >  1:733075/1‑71 (MQ=255)                                                               cAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGcc     >  1:2073276/1‑71 (MQ=255)                                                                 gCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACgg                            <  1:658587/46‑1 (MQ=255)                                                                 gCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCctg   >  1:2342183/1‑71 (MQ=255)                                                                  cGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGatgatg                  >  1:685482/1‑55 (MQ=255)                                                                  cGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGatgatg                  >  1:2844866/1‑55 (MQ=255)                                                                  cGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCctgc  >  1:2147283/1‑71 (MQ=255)                                                                     |                                                                    GATCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCCTGC  >  minE/2632902‑2633037

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