Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,104,294 A→T 10.2% G339G (GGA→GGT ydiF → fused predicted acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit and beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,104,2940AT10.2% 130.5 / 7.7 49G339G (GGA→GGTydiFfused predicted acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit and beta subunit
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base T (2/3);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACG  >  minE/1104232‑1104357
                                                              |                                                               
gCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTc                                                         <  1:3803820/71‑1 (MQ=255)
gCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTc                                                         <  1:3234401/71‑1 (MQ=255)
      gaagaaGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGTATTACTTca                                                        <  1:274938/66‑1 (MQ=255)
      gaagaaGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGTATTACTTca                                                        <  1:1577789/66‑1 (MQ=255)
      gaagaaGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGTATTACTTca                                                        <  1:3675666/66‑1 (MQ=255)
        agaagGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGTATTACTTcaca                                                      >  1:3293747/1‑66 (MQ=255)
        agaagGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGTATTACTTcaca                                                      >  1:1600240/1‑66 (MQ=255)
          aagGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATc                                               <  1:3440463/71‑1 (MQ=255)
            ggTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGc                                             <  1:98194/71‑1 (MQ=255)
                gtgCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGaa                                                                >  1:3685791/1‑48 (MQ=255)
                   cTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGgc                                      >  1:608984/1‑71 (MQ=255)
                     gatgaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCttt                                         <  1:4049897/66‑1 (MQ=255)
                       tgaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGcgcg                                   >  1:2134421/1‑70 (MQ=255)
                       tgaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGcgcg                                   >  1:1388435/1‑70 (MQ=255)
                        gaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTcaca                                                      >  1:2114004/1‑50 (MQ=255)
                        gaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAg                                                     >  1:2292029/1‑51 (MQ=255)
                         aCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAggg                                                   >  1:2674863/1‑52 (MQ=255)
                           tttATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGt                              >  1:1749346/1‑71 (MQ=255)
                           tttATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGt                              >  1:1243626/1‑71 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGTATCGCCTTTGGCGCGAACGTg                             <  1:1137005/70‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACTTg                             <  1:602060/70‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             <  1:4057134/70‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             <  1:229905/70‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             <  1:3209840/70‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             >  1:296514/1‑70 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             <  1:2890179/70‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             >  1:1893764/1‑70 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTg                             >  1:98627/1‑70 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGAGCGCCTTTGGCGCGAACGt                              <  1:2077400/69‑1 (MQ=255)
                             tATTCTGACGGTAGAAACAGGTCCGATTGGGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGt                              <  1:1670817/68‑1 (MQ=255)
                              aTTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGc                                             <  1:931635/53‑1 (MQ=255)
                               ttCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGaa                           <  1:4094874/70‑1 (MQ=255)
                               ttCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGaa                           <  1:1581941/70‑1 (MQ=255)
                                   gACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCtt                                          >  1:1681277/1‑51 (MQ=255)
                                   gACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAccc                      >  1:3214547/1‑71 (MQ=255)
                                   gACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAccc                      >  1:1233271/1‑71 (MQ=255)
                                   gACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAcc                       >  1:3455194/1‑70 (MQ=255)
                                       gTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGaa                                 <  1:4255199/56‑1 (MQ=255)
                                         aGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCa                >  1:453735/1‑71 (MQ=255)
                                         aGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCa                >  1:2410247/1‑71 (MQ=255)
                                                ggTCCGATTGGCGGAATTACTTCACACGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGa         <  1:3666696/71‑1 (MQ=255)
                                                 gTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTg           >  1:3874960/1‑68 (MQ=255)
                                                 gTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTg           >  1:3295561/1‑68 (MQ=255)
                                                    cGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGatat      <  1:4183701/70‑1 (MQ=255)
                                                    cGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGatat      <  1:3072258/70‑1 (MQ=255)
                                                    cGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGatat      <  1:1707606/70‑1 (MQ=255)
                                                      aTTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGa    >  1:406936/1‑70 (MQ=255)
                                                       ttGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACg  >  1:4023923/1‑71 (MQ=255)
                                                       ttGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACg  <  1:1489000/71‑1 (MQ=255)
                                                         ggCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGa         >  1:3886266/1‑62 (MQ=255)
                                                              |                                                               
GCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACG  >  minE/1104232‑1104357

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: