Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,371,048 T→A 9.1% intergenic (+44/+337) 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,0480TA 9.1% 183.0 / 9.4 66intergenic (+44/+337)fadA/fre3‑ketoacyl‑CoA thiolase/flavin reductase
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base A (2/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

ATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATGC  >  minE/2370978‑2371118
                                                                      |                                                                       
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)
        ggTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAc‑c                                                                <  1:1825849/71‑1 (MQ=255)
          tgtTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑Caa                                                              <  1:2351709/71‑1 (MQ=255)
          tgtTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑Caa                                                              <  1:2347182/71‑1 (MQ=255)
            tttGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATg                                                            <  1:1481760/71‑1 (MQ=255)
              tGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGt                                                           <  1:11067/70‑1 (MQ=255)
              tGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGt                                                           <  1:3644893/70‑1 (MQ=255)
               gAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAg                                                         >  1:2723594/1‑71 (MQ=255)
               gAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAg                                                         >  1:165887/1‑71 (MQ=255)
                aGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAgg                                                        <  1:2719369/71‑1 (MQ=255)
                aGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAg                                                         <  1:3094135/70‑1 (MQ=255)
                 gCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGc                                                       >  1:45117/1‑71 (MQ=255)
                   gggTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCt                                                     <  1:353922/71‑1 (MQ=255)
                     gTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGAt                                                                         <  1:2970372/50‑1 (MQ=255)
                       ttAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGAt                                                  >  1:2739348/1‑70 (MQ=255)
                       ttAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGAt                                                  >  1:912343/1‑70 (MQ=255)
                        tAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTa                                                          >  1:1821039/1‑61 (MQ=255)
                           ccACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGAc                                             <  1:3570427/71‑1 (MQ=255)
                            cACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGc                                                                    <  1:2324633/48‑1 (MQ=255)
                            cACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGc                                                       >  1:2232994/1‑60 (MQ=255)
                             aCTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAg                                                         <  1:2602251/57‑1 (MQ=255)
                             aCTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGAcgc                                           >  1:3433740/1‑71 (MQ=255)
                             aCTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGAcg                                            >  1:3468599/1‑70 (MQ=255)
                                tGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACgcgc                                         >  1:1605678/1‑70 (MQ=255)
                                  ccGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGa                                              >  1:186351/1‑63 (MQ=255)
                                  ccGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAg                                      <  1:661895/71‑1 (MQ=255)
                                   cGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGc                                     <  1:147986/71‑1 (MQ=255)
                                     gATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGt                                   >  1:3149180/1‑71 (MQ=255)
                                      aTGCGGCGTGAACGCCTTATCCAGCCTACCAAACCCGCACTGATTGTAGGCATGATAAGACGCGCAAGCGt                                   <  1:2634653/71‑1 (MQ=2)
                                       tGCGGCGTGAACGCCTTGTCCGGCCTACAAAACCGAAAC‑GTATGTAGGCCTGATAAGACGCGTCAGCGTc                                  >  1:2007651/1‑70 (MQ=9)
                                       tGCGGCGTGAACGCCTTGTCCGGCCTACAAAACCGAAAC‑GTATGTAGGCCTGATAAGACGCGTCAGCGTCg                                 <  1:3269882/71‑1 (MQ=9)
                                       tGCGGCGTGAACGCCTTGTCCGGCCTACAAAACCGAAAC‑GTATGTAGGCCTGATAAGACGCGTCAGCGTCg                                 <  1:2592090/71‑1 (MQ=9)
                                       tGCGGCGTGAACGCCTTATCCAGCCTACCAAACCCGCACTGATTGTAGGCATGATAAGACGCGCAAGCGTc                                  >  1:1553180/1‑71 (MQ=2)
                                         cggcgTGAACGCCTTGTCCGGCCTACAAAACCGAAAC‑GTATGTAGGCCTGATAAGACGCGTCAGCGTCGCa                               <  1:878675/71‑1 (MQ=9)
                                          ggcgTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGc                                <  1:3178851/69‑1 (MQ=255)
                                            cgTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                              >  1:3090557/1‑69 (MQ=255)
                                            cgTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                              >  1:2460578/1‑69 (MQ=255)
                                            cgTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                              >  1:2869542/1‑69 (MQ=255)
                                             gTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAg                           >  1:2394320/1‑71 (MQ=255)
                                              tGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAg                                               <  1:2461296/50‑1 (MQ=255)
                                              tGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGACAGCGTCGCATCAg                           <  1:576410/70‑1 (MQ=255)
                                                aaCGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGca                        >  1:3861892/1‑71 (MQ=255)
                                                aaCGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGca                        <  1:3238048/71‑1 (MQ=255)
                                                  cGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAg                                                         >  1:3689881/1‑36 (MQ=255)
                                                  cGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATaa                                                >  1:877349/1‑45 (MQ=255)
                                                  cGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGc                         >  1:2480595/1‑68 (MQ=255)
                                                   gCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcatca                     <  1:3754189/71‑1 (MQ=255)
                                                    ccTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcatcat                    <  1:1729850/71‑1 (MQ=255)
                                                      ttATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATg                   >  1:3785551/1‑70 (MQ=255)
                                                       tATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCa                 >  1:612392/1‑71 (MQ=255)
                                                       tATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCa                 >  1:2640045/1‑71 (MQ=255)
                                                       tATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCa                 >  1:878056/1‑71 (MQ=255)
                                                          ccGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCAcc               >  1:1794831/1‑70 (MQ=255)
                                                           cGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCaa             <  1:1384032/71‑1 (MQ=255)
                                                             gTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCg                                    <  1:2517061/46‑1 (MQ=255)
                                                             gTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAAtt           <  1:3476941/71‑1 (MQ=255)
                                                              tCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                             >  1:432085/1‑52 (MQ=255)
                                                                tACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATg                   <  1:3685932/60‑1 (MQ=255)
                                                                 aCCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTg          <  1:3366187/68‑1 (MQ=255)
                                                                  ccGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGc         <  1:1577020/68‑1 (MQ=255)
                                                                  ccGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGc         <  1:3083887/68‑1 (MQ=255)
                                                                   cGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGa     >  1:1352357/1‑71 (MQ=255)
                                                                   cGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGa     >  1:3214809/1‑71 (MQ=255)
                                                                     aTCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGt                                   >  1:2423778/1‑39 (MQ=255)
                                                                     aTCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATg   >  1:1832714/1‑71 (MQ=255)
                                                                      tCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATgc  <  1:2800128/71‑1 (MQ=255)
                                                                      |                                                                       
ATTGCGACGGTGTTTGAGCGGGTTTAACCACTTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACCGATCCGGCAC‑CAATGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATGC  >  minE/2370978‑2371118

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: ‑