Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,371,040 T→C 8.3% intergenic (+36/+345) fadA → / ← fre 3‑ketoacyl‑CoA thiolase/flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,371,0400TC 8.3% 196.2 / 6.9 72intergenic (+36/+345)fadA/fre3‑ketoacyl‑CoA thiolase/flavin reductase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (2/4);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TCAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATT  >  minE/2370971‑2371109
                                                                     |                                                                      
tCAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTc                                                                       <  1:2320286/71‑1 (MQ=255)
tCAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTc                                                                       <  1:135713/71‑1 (MQ=255)
 cAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGt                                                                        >  1:2593628/1‑69 (MQ=255)
 cAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGt                                                                        >  1:2062439/1‑69 (MQ=255)
 cAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTc                                                                       >  1:1554063/1‑70 (MQ=255)
  aGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCt                                                                      <  1:1570156/70‑1 (MQ=255)
  aGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCt                                                                      <  1:2553193/70‑1 (MQ=255)
  aGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTa                                                                     >  1:3869089/1‑71 (MQ=255)
  aGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTa                                                                     <  1:2167443/71‑1 (MQ=255)
  aGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTa                                                                     >  1:1336857/1‑71 (MQ=255)
  aGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTa                                                                     >  1:3171271/1‑71 (MQ=255)
    ggTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTAcc                                                                   <  1:3938615/71‑1 (MQ=255)
     gTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCg                                                                  <  1:166410/71‑1 (MQ=255)
       aTTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGAt                                                                >  1:2384978/1‑71 (MQ=255)
        ttGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGAt                                                                <  1:1291123/70‑1 (MQ=255)
        ttGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATc                                                               >  1:342801/1‑71 (MQ=255)
        ttGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATc                                                               >  1:1706256/1‑71 (MQ=255)
          gCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCg                                                             >  1:1077548/1‑71 (MQ=255)
           cGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTa                                                                     >  1:2245605/1‑62 (MQ=255)
               ggTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAcc                                                       <  1:1825849/71‑1 (MQ=255)
                 tgtTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCaa                                                     <  1:2351709/71‑1 (MQ=255)
                 tgtTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCaa                                                     <  1:2347182/71‑1 (MQ=255)
                   tttGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATg                                                   <  1:1481760/71‑1 (MQ=255)
                     tGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGt                                                  <  1:3644893/70‑1 (MQ=255)
                     tGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGt                                                  <  1:11067/70‑1 (MQ=255)
                      gAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAg                                                >  1:165887/1‑71 (MQ=255)
                      gAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAg                                                >  1:2723594/1‑71 (MQ=255)
                       aGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAgg                                               <  1:2719369/71‑1 (MQ=255)
                       aGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAg                                                <  1:3094135/70‑1 (MQ=255)
                        gCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGc                                              >  1:45117/1‑71 (MQ=255)
                          gggTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCt                                            <  1:353922/71‑1 (MQ=255)
                            gTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGAt                                                                <  1:2970372/50‑1 (MQ=255)
                              ttAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGAt                                         >  1:912343/1‑70 (MQ=255)
                              ttAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGAt                                         >  1:2739348/1‑70 (MQ=255)
                               tAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTa                                                 >  1:1821039/1‑61 (MQ=255)
                                  ccACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGAc                                    <  1:3570427/71‑1 (MQ=255)
                                   cACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGc                                                           <  1:2324633/48‑1 (MQ=255)
                                   cACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGc                                              >  1:2232994/1‑60 (MQ=255)
                                    aCTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAg                                                <  1:2602251/57‑1 (MQ=255)
                                    aCTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGAcgc                                  >  1:3433740/1‑71 (MQ=255)
                                    aCTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGAcg                                   >  1:3468599/1‑70 (MQ=255)
                                       tGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACgcgc                                >  1:1605678/1‑70 (MQ=255)
                                         ccGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGa                                     >  1:186351/1‑63 (MQ=255)
                                         ccGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAg                             <  1:661895/71‑1 (MQ=255)
                                          cGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGc                            <  1:147986/71‑1 (MQ=255)
                                            gATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGt                          >  1:3149180/1‑71 (MQ=255)
                                             aTGCGGCGTGAACGCCTTATCCAGCCTACCAAACCCGCACTGATTGTAGGCATGATAAGACGCGCAAGCGt                          <  1:2634653/71‑1 (MQ=2)
                                              tGCGGCGTGAACGCCTTGTCCGGCCTACAAAACCGAAACGTATGTAGGCCTGATAAGACGCGTCAGCGTc                         >  1:2007651/1‑70 (MQ=9)
                                              tGCGGCGTGAACGCCTTGTCCGGCCTACAAAACCGAAACGTATGTAGGCCTGATAAGACGCGTCAGCGTCg                        <  1:2592090/71‑1 (MQ=9)
                                              tGCGGCGTGAACGCCTTGTCCGGCCTACAAAACCGAAACGTATGTAGGCCTGATAAGACGCGTCAGCGTCg                        <  1:3269882/71‑1 (MQ=9)
                                              tGCGGCGTGAACGCCTTATCCAGCCTACCAAACCCGCACTGATTGTAGGCATGATAAGACGCGCAAGCGTc                         >  1:1553180/1‑71 (MQ=2)
                                                cggcgTGAACGCCTTGTCCGGCCTACAAAACCGAAACGTATGTAGGCCTGATAAGACGCGTCAGCGTCGCa                      <  1:878675/71‑1 (MQ=9)
                                                 ggcgTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGc                       <  1:3178851/69‑1 (MQ=255)
                                                   cgTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                     >  1:2460578/1‑69 (MQ=255)
                                                   cgTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                     >  1:3090557/1‑69 (MQ=255)
                                                   cgTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                     >  1:2869542/1‑69 (MQ=255)
                                                    gTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAg                  >  1:2394320/1‑71 (MQ=255)
                                                     tGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAg                                      <  1:2461296/50‑1 (MQ=255)
                                                     tGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGACAGCGTCGCATCAg                  <  1:576410/70‑1 (MQ=255)
                                                       aaCGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGca               >  1:3861892/1‑71 (MQ=255)
                                                       aaCGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGca               <  1:3238048/71‑1 (MQ=255)
                                                         cGCCTTATCCGGTCTACCGATCCGGCACCAATGTAg                                                >  1:3689881/1‑36 (MQ=255)
                                                         cGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATaa                                       >  1:877349/1‑45 (MQ=255)
                                                         cGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGc                >  1:2480595/1‑68 (MQ=255)
                                                          gCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcatca            <  1:3754189/71‑1 (MQ=255)
                                                           ccTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcatcat           <  1:1729850/71‑1 (MQ=255)
                                                             ttATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATg          >  1:3785551/1‑70 (MQ=255)
                                                              tATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCa        >  1:612392/1‑71 (MQ=255)
                                                              tATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCa        >  1:878056/1‑71 (MQ=255)
                                                              tATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCa        >  1:2640045/1‑71 (MQ=255)
                                                                 ccGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCAcc      >  1:1794831/1‑70 (MQ=255)
                                                                  cGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCaa    <  1:1384032/71‑1 (MQ=255)
                                                                    gTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCg                           <  1:2517061/46‑1 (MQ=255)
                                                                    gTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAAtt  <  1:3476941/71‑1 (MQ=255)
                                                                     tCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                    >  1:432085/1‑52 (MQ=255)
                                                                     |                                                                      
TCAGGGTATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCACCAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATT  >  minE/2370971‑2371109

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: