Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,544,443 G→A 9.7% L427F (CTT→TTT)  yfcY ← beta‑ketoacyl‑CoA thiolase, anaerobic, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,544,4430GA 9.7% 188.5 / 8.3 62L427F (CTT→TTT) yfcYbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit
Reads supporting (aligned to +/- strand):  ref base G (40/16);  new base A (3/3);  total (43/19)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGTCTG  >  minE/1544373‑1544512
                                                                      |                                                                     
cAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAg                                                                       >  1:526835/1‑71 (MQ=255)
cAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAg                                                                       >  1:2879699/1‑71 (MQ=255)
cAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAg                                                                       >  1:1490608/1‑71 (MQ=255)
cAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAg                                                                       >  1:1971234/1‑71 (MQ=255)
   aCGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAAccc                                                                    <  1:406794/71‑1 (MQ=255)
     gAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCAc                                                                  >  1:373839/1‑71 (MQ=255)
     gAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCAc                                                                  >  1:2776007/1‑71 (MQ=255)
        cATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGccc                                                                    <  1:600432/66‑1 (MQ=255)
        cATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGccc                                                                    <  1:2644965/66‑1 (MQ=255)
        cATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGccc                                                                    <  1:1850001/66‑1 (MQ=255)
           tAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGgcag                                                            >  1:933054/1‑71 (MQ=255)
            aaGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGgcag                                                            <  1:716989/70‑1 (MQ=255)
            aaGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGgcag                                                            <  1:2615948/70‑1 (MQ=255)
            aaGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGgcag                                                            <  1:2559764/70‑1 (MQ=255)
            aaGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAACCCACCGgcag                                                            >  1:897164/1‑70 (MQ=255)
                  aaaCGCTGATGTCAATTCCATTAATCCGCCTCCAGAACCATTGCCGCGCCAAAccc                                                                    <  1:3456386/56‑1 (MQ=255)
                   aaCGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCg                                                    >  1:1097775/1‑71 (MQ=255)
                   aaCGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCg                                                    >  1:2556445/1‑71 (MQ=255)
                     cGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGt                                                  <  1:105428/71‑1 (MQ=255)
                     cGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGt                                                  <  1:1077301/71‑1 (MQ=255)
                        tGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGt                                                  <  1:1699040/68‑1 (MQ=255)
                          aTGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACt                                              >  1:3085767/1‑70 (MQ=255)
                          aTGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACt                                              >  1:2883741/1‑70 (MQ=255)
                           tGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTaa                                            <  1:30986/71‑1 (MQ=255)
                            gTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTaaa                                           >  1:884335/1‑71 (MQ=255)
                            gTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTaaa                                           >  1:251226/1‑71 (MQ=255)
                            gTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTaaa                                           >  1:3057932/1‑71 (MQ=255)
                            gTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTaaa                                           >  1:3365250/1‑71 (MQ=255)
                            gTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTaaa                                           >  1:3305340/1‑71 (MQ=255)
                               aTTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAAcc                                         >  1:3374534/1‑70 (MQ=255)
                                 ttCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAAccc                                                                    <  1:2399327/41‑1 (MQ=255)
                                  tCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAAt                                     >  1:3801029/1‑71 (MQ=255)
                                  tCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAAt                                     >  1:691016/1‑71 (MQ=255)
                                  tCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAAt                                     >  1:3305740/1‑71 (MQ=255)
                                    cATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATc                                    >  1:2374563/1‑70 (MQ=255)
                                      ttattCCGCCTCCAGAACCATTGCCGCGCCAAGcc                                                                     >  1:2692643/1‑35 (MQ=255)
                                      ttattCCGCCTCCAGAACCATTGCCGCGCCAAGcc                                                                     >  1:4102959/1‑35 (MQ=255)
                                      ttattCCGCCTCCAGAACCATTGCCGCGCCAAGcc                                                                     >  1:10997/1‑35 (MQ=255)
                                      ttattCCGCCTCCAGAACCATTGCCGCGCCAAGcc                                                                     >  1:3022493/1‑35 (MQ=255)
                                        attCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCAccg                               >  1:3197573/1‑71 (MQ=255)
                                        attCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCAccg                               >  1:2751368/1‑71 (MQ=255)
                                         ttCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGgcagca                                                          >  1:27634/1‑43 (MQ=255)
                                         ttCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGgcagca                                                          >  1:172161/1‑43 (MQ=255)
                                            cGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACcgccgc                           <  1:511555/71‑1 (MQ=255)
                                            cGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACcgccgc                           <  1:3965515/71‑1 (MQ=255)
                                                  cAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGc                     <  1:2826270/71‑1 (MQ=255)
                                                    gAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGa                         >  1:2224652/1‑65 (MQ=255)
                                                    gAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCg                    >  1:1158149/1‑70 (MQ=255)
                                                    gAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCg                    >  1:2334240/1‑70 (MQ=255)
                                                    gAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCg                    >  1:2016463/1‑70 (MQ=255)
                                                    gAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCg                    >  1:2113162/1‑70 (MQ=255)
                                                    gAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCg                    >  1:3124705/1‑70 (MQ=255)
                                                        cATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGtt               <  1:4076791/71‑1 (MQ=255)
                                                        cATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGtt               <  1:4122747/71‑1 (MQ=255)
                                                               gcgcCAAACCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCaa        >  1:3640326/1‑71 (MQ=255)
                                                               gcgcCAAACCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCaa        >  1:439662/1‑71 (MQ=255)
                                                                cgcCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCaa        <  1:2171322/70‑1 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATg      >  1:1958386/1‑70 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGt     >  1:47144/1‑71 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGt     >  1:1752510/1‑71 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGt     >  1:1218858/1‑71 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGt     >  1:2831239/1‑71 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGt     >  1:1719324/1‑71 (MQ=255)
                                                                  ccAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGt     >  1:1894017/1‑71 (MQ=255)
                                                                     aGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGTCTg  >  1:601400/1‑71 (MQ=255)
                                                                     aGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGTCTg  >  1:898419/1‑71 (MQ=255)
                                                                      |                                                                     
CAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCCGCGCCAAGCCCACCGGCAGCACAGGCGGTAACTAAACCAAATCCACCGCCGCGACGGCGAAGTTCATGCAATGTCTG  >  minE/1544373‑1544512

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: