Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 67,673 G→A 39.2% I199I (ATC→ATT leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE67,6730GA39.2% 40.1 / 51.2 79I199I (ATC→ATTleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (10/38);  new base A (18/13);  total (28/51)
Fisher's exact test for biased strand distribution p-value = 1.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGG  >  minE/67608‑67741
                                                                 |                                                                    
ccGGCATGTACCTCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGAt                                                                    <  1:2529707/68‑1 (MQ=255)
ccGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGAt                                                                    <  1:2669528/68‑1 (MQ=255)
ccGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGAt                                                                    <  1:1833472/68‑1 (MQ=255)
ccGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGAt                                                                    <  1:253081/68‑1 (MQ=255)
         aCCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCttt                                                           >  1:1718168/1‑68 (MQ=255)
         aCCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTg                                                          >  1:1381416/1‑69 (MQ=255)
            gccgccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCg                                                       <  1:1699105/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:1015221/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:3142185/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:137285/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:303418/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:2732571/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:1885546/69‑1 (MQ=255)
                ccAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:2360187/69‑1 (MQ=255)
                  aGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:508867/67‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:2564201/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:1845691/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:2058131/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:2440299/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:510729/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:326565/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:2652634/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:997249/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:3098197/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:3261161/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:1816628/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:1445171/50‑1 (MQ=255)
                  aGTGAGTTGCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAAAt                                                                    <  1:15429/50‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGAt                                                                    <  1:1876641/48‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAgg                                               <  1:1001280/69‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAgg                                               <  1:1153859/69‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAgg                                               <  1:2951574/69‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAgg                                               <  1:2700013/69‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAgg                                               <  1:1810484/69‑1 (MQ=255)
                    tGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAgg                                               <  1:1756095/69‑1 (MQ=255)
                      aGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGt                                                   <  1:2867524/63‑1 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATTGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:1290412/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCAt                                    >  1:150569/1‑60 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:1048628/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:1080007/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:591725/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:503052/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:1855736/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:305895/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:751922/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:2038096/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:82989/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:2005381/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:162305/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:2320084/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:89100/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:245435/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACGCGTTCATACAGGCCg                            >  1:2223609/1‑68 (MQ=255)
                                        ccAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTGGGTACACGTTCATACAGGCCg                            >  1:2757782/1‑68 (MQ=255)
                                             tcgtcgTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCtgatg                      >  1:939388/1‑69 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCAcgc                                          >  1:215848/1‑37 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:1981967/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:1920646/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:2556623/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:1663497/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:2288463/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:92317/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGtt                                       <  1:1206193/40‑1 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAg                                >  1:3407171/1‑47 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAg                                >  1:2059776/1‑47 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAg                                >  1:649631/1‑47 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAg                                >  1:2359463/1‑47 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAg                                >  1:1532439/1‑47 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAg                                >  1:2486837/1‑47 (MQ=255)
                                                         tGTACGGAGATAATGGCTTTGTCGATGTTAGGCACACGtt                                       <  1:1251236/40‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGTCGAACTCa       <  1:986783/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCaa      <  1:2062526/69‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:541021/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:508036/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:3338685/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:3038067/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:1631975/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:2237079/68‑1 (MQ=255)
                                                             cGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCa       <  1:954653/68‑1 (MQ=255)
                                                               gagaTAATGGCTTTGTCGATGTTAGGCACTCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAAc    >  1:22351/1‑69 (MQ=255)
                                                               gagaTAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAAc    >  1:510134/1‑69 (MQ=255)
                                                               gagaTAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAAc    >  1:742100/1‑69 (MQ=255)
                                                                 gaTAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACgg  <  1:1389715/69‑1 (MQ=255)
                                                                 gaTAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACgg  <  1:1567606/69‑1 (MQ=255)
                                                                 gaTAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACgg  <  1:1751032/69‑1 (MQ=255)
                                                                 |                                                                    
CCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGG  >  minE/67608‑67741

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: