Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,794,549 A→G 24.9% I264I (ATT→ATC plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,794,5490AG24.9% 82.9 / 24.7 72I264I (ATT→ATCplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (12/6);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTT  >  minE/2794482‑2794616
                                                                   |                                                                   
cGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGa                                                                     <  1:2553430/68‑1 (MQ=255)
cGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGa                                                                     <  1:63872/68‑1 (MQ=255)
    aTGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatcat                                                                >  1:2145913/1‑69 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:2470354/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:2425977/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:1138805/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:815590/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:124792/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:2413908/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:2262579/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:1457447/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:3407683/69‑1 (MQ=255)
         cGTTGATGCCCTGGTAAAGTCGGTTCCAGATGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:954093/69‑1 (MQ=255)
           ttGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:1075547/67‑1 (MQ=255)
                gCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAgaga                                                     <  1:2147960/68‑1 (MQ=255)
                gCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAgaga                                                     <  1:3111888/68‑1 (MQ=255)
                gCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAgaga                                                     <  1:41916/68‑1 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                    >  1:1491695/1‑62 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                              >  1:1863487/1‑68 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                              >  1:1313854/1‑68 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                              >  1:3093979/1‑68 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                              >  1:3316383/1‑68 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                              >  1:2216407/1‑68 (MQ=255)
                       tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                              >  1:762632/1‑68 (MQ=255)
                        aaaGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAGTCATCTCGTAAGAGAAATTagcc                                             >  1:396035/1‑68 (MQ=255)
                        aaaGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                             >  1:1283079/1‑68 (MQ=255)
                        aaaGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                             >  1:2136913/1‑68 (MQ=255)
                            tCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                        <  1:1878168/69‑1 (MQ=255)
                            tCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                        <  1:255610/69‑1 (MQ=255)
                            tCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                        <  1:251533/69‑1 (MQ=255)
                              ggTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:2653467/48‑1 (MQ=255)
                                     ggTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGt                                                           <  1:230474/41‑1 (MQ=255)
                                     ggTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCAt                               <  1:3127829/69‑1 (MQ=255)
                                       tGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                             >  1:370626/1‑69 (MQ=255)
                                       tGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                             >  1:1823862/1‑69 (MQ=255)
                                       tGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                             >  1:1544298/1‑69 (MQ=255)
                                                  aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                  <  1:2493952/69‑1 (MQ=255)
                                                  aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                  <  1:49226/69‑1 (MQ=255)
                                                  aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                  >  1:3202963/1‑69 (MQ=255)
                                                  aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                  <  1:1720096/69‑1 (MQ=255)
                                                     cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc               >  1:82901/1‑69 (MQ=255)
                                                     cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc               >  1:3182644/1‑69 (MQ=255)
                                                     cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc               >  1:2514196/1‑69 (MQ=255)
                                                     cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc               >  1:3266505/1‑69 (MQ=255)
                                                     cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc               >  1:1748501/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTTCATCAGCGCAATCGCGTTctgc            >  1:2470639/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCTCGTTctgc            >  1:464302/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:2990191/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:100190/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:2286494/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:438512/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:2257310/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:2166920/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:3386694/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:1968414/1‑69 (MQ=255)
                                                        tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc            >  1:1925577/1‑69 (MQ=255)
                                                        tcagtaagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                              >  1:1294529/1‑51 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc       >  1:1868386/1‑69 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc       >  1:1179197/1‑69 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc       >  1:733123/1‑69 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc       >  1:1867776/1‑69 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc       >  1:3222560/1‑69 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc       >  1:3166649/1‑69 (MQ=255)
                                                             cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTACATCAGTGCAATCGCGTTCTGCTgcgc       >  1:2348052/1‑69 (MQ=255)
                                                                gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt     <  1:1669036/68‑1 (MQ=255)
                                                                gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt     <  1:595794/68‑1 (MQ=255)
                                                                gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt     <  1:165557/68‑1 (MQ=255)
                                                                gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt     <  1:71310/68‑1 (MQ=255)
                                                                gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt     <  1:2503181/68‑1 (MQ=255)
                                                                gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt     <  1:1019032/68‑1 (MQ=255)
                                                                  gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt  <  1:2856228/69‑1 (MQ=255)
                                                                  gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt  <  1:650296/69‑1 (MQ=255)
                                                                  gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt  <  1:266845/69‑1 (MQ=255)
                                                                  gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt  <  1:1057630/69‑1 (MQ=255)
                                                                   |                                                                   
CGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTT  >  minE/2794482‑2794616

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