Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,026,074 A→G 18.2% R56R (CGT→CGC plsC ← 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,026,0740AG18.2% 131.6 / 18.8 87R56R (CGT→CGCplsC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/27);  new base G (14/2);  total (58/29)
Fisher's exact test for biased strand distribution p-value = 7.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TCATAGTTGTTCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGT  >  minE/2026009‑2026133
                                                                 |                                                           
tCATAGTTGTTCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACgg                                                          >  1:1831078/1‑69 (MQ=255)
tCATAGTTGTTCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACgg                                                          >  1:1784877/1‑69 (MQ=255)
      ttgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTc                                                     >  1:1765111/1‑68 (MQ=255)
      ttgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:922403/1‑69 (MQ=255)
      ttgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2492909/1‑69 (MQ=255)
      ttgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2473377/1‑69 (MQ=255)
      ttgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2290745/1‑69 (MQ=255)
      ttgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:3171847/1‑69 (MQ=255)
       tgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    <  1:697058/68‑1 (MQ=255)
       tgttCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    <  1:2863510/68‑1 (MQ=255)
           cTGGTGGTTAGCAATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGCAGGTTTGCGGCACTCCACtt                                                <  1:2849073/68‑1 (MQ=255)
            tggtggTTAGCAATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGCAGGTTTGCGGCACTCCACtt                                                <  1:1086279/67‑1 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCa                                             >  1:2435182/1‑68 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAg                                            >  1:258952/1‑69 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAg                                            >  1:89299/1‑69 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAg                                            >  1:2263169/1‑69 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAg                                            >  1:3216169/1‑69 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAg                                            >  1:2891435/1‑69 (MQ=255)
              gtggtTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAg                                            >  1:3139856/1‑69 (MQ=255)
                             gatagCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAgcggc                               >  1:70398/1‑67 (MQ=255)
                             gatagCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAgcggc                               >  1:2184120/1‑67 (MQ=255)
                             gatagCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAgcggc                               >  1:207070/1‑67 (MQ=255)
                             gatagCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAgcggc                               >  1:1434321/1‑67 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTc                                                     >  1:2923946/1‑42 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2033353/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2600205/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:71915/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2403317/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:2384089/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:3264506/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:1400889/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:219496/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:1724531/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:420788/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:650791/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCa                                                    >  1:1458153/1‑43 (MQ=255)
                                agCATTGCCGTAGCTTTCAGCGTCAGTAGGTTTACGGCACTc                                                     >  1:2927818/1‑42 (MQ=255)
                                   aTTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACg                       <  1:3248247/69‑1 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:589093/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:3328055/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:3347250/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:2524675/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:2696406/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:34105/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:1597981/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:1616941/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:1131464/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:1329084/1‑68 (MQ=255)
                                      gCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGc                     >  1:2109869/1‑68 (MQ=255)
                                         gTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAATCAGCGGCGCAAGACGGCCaaa                 <  1:1610365/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGGCAAACa               <  1:1361266/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:211478/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:3352703/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:1336899/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:1725450/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:3056169/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:992567/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:948585/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:250744/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:2631474/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:2705186/69‑1 (MQ=255)
                                           aGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:1188259/69‑1 (MQ=255)
                                            gCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAAct               <  1:9166/68‑2 (MQ=255)
                                            gCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:2168344/68‑1 (MQ=255)
                                            gCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:927853/68‑1 (MQ=255)
                                            gCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:1807576/68‑1 (MQ=255)
                                            gCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:2857404/68‑1 (MQ=255)
                                            gCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACa               <  1:2232495/68‑1 (MQ=255)
                                              tttCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACata             <  1:713901/68‑1 (MQ=255)
                                              tttCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACata             <  1:1475626/68‑1 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCaaa     >  1:2529612/4‑66 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAgg   >  1:881776/4‑68 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAgg   >  1:1634744/4‑68 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAgg   >  1:2651517/4‑68 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:34158/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:1627967/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:1675382/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:1676702/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:853215/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:2507776/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:2675137/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:2514644/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTCAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:1073609/4‑69 (MQ=255)
                                                        tgcaGGTTTGCGGCACTCCACTTTAAGGCCGAAAAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  >  1:587672/4‑69 (MQ=255)
                                                        tGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  <  1:2644831/69‑1 (MQ=255)
                                                        tGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  <  1:2588707/69‑1 (MQ=255)
                                                          tAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGt  <  1:162485/67‑1 (MQ=255)
                                                                 |                                                           
TCATAGTTGTTCTGGTGGTTAGCGATATAGATAGCATTGCCGTAGCTTTCAGCGTCTGTAGGTTTACGGCACTCAACTTTCAGGCCAAACAGCGGCGCAAGACGGCCAAACATATGCCCAAAGGT  >  minE/2026009‑2026133

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: