Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 466,055 A→T 62.5% M299L (ATG→TTG)  sucA → 2‑oxoglutarate decarboxylase, thiamin‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE466,0550AT62.5% 26.5 / 32.5 40M299L (ATG→TTG) sucA2‑oxoglutarate decarboxylase, thiamin‑requiring
Reads supporting (aligned to +/- strand):  ref base A (9/6);  new base T (9/16);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCC  >  minE/465993‑466119
                                                              |                                                                
tGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCt                                                            >  1:2731412/1‑69 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1040075/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:811839/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:365473/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:340705/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:3163622/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:2678317/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:2607510/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:2025492/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1956525/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1919650/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1823330/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1124570/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1588972/68‑1 (MQ=255)
             tgCCGGTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:1439760/68‑1 (MQ=255)
                  gTAAACATAAAGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:518559/64‑1 (MQ=255)
                           aaGAACACCTCGGCACGGGCGACGTGAAATACCACTTGGGCTTCTCGTCTGACtt                                               <  1:508240/55‑1 (MQ=255)
                                    tCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGc                        <  1:1016058/69‑1 (MQ=255)
                                    tCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGc                        <  1:2305876/69‑1 (MQ=255)
                                    tCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGc                        <  1:2398799/69‑1 (MQ=255)
                                    tCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGc                        <  1:1817537/69‑1 (MQ=255)
                                     cGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGc                        <  1:2850615/68‑1 (MQ=255)
                                     cGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGc                        <  1:1134021/68‑1 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:1592980/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:18899/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:777993/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:725054/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:1678536/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:452163/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:431662/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:1172232/1‑69 (MQ=255)
                                                 cGTGAAATACCACTTGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTTCACctggcgctggc           >  1:3043312/1‑69 (MQ=255)
                                                  gTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGCTGGCGGCCTGGTGCACctggcgctggcg          >  1:2721677/1‑69 (MQ=255)
                                                  gTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACctggcgctggcg          >  1:2199556/1‑69 (MQ=255)
                                                  gTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACctggcgctggcg          >  1:2609439/1‑69 (MQ=255)
                                                  gTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACctggcgctggcg          >  1:1713100/1‑69 (MQ=255)
                                                          ccACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAAccc  >  1:420028/1‑69 (MQ=255)
                                                          ccACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAAccc  >  1:3264609/1‑69 (MQ=255)
                                                          ccACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAAccc  >  1:2507961/1‑69 (MQ=255)
                                                          ccACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAAccc  >  1:1210865/1‑69 (MQ=255)
                                                              |                                                                
TGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCC  >  minE/465993‑466119

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