Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,064 G→T 17.8% T730T (ACC→ACA)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,0640GT17.8% 153.2 / 14.9 56T730T (ACC→ACA) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (28/18);  new base T (4/6);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTA  >  minE/607996‑608130
                                                                    |                                                                  
gaaataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGt                                                                  >  1:1494929/3‑71 (MQ=255)
gTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                                  >  1:964003/1‑71 (MQ=255)
gTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                                  >  1:580131/1‑71 (MQ=255)
 taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                                  >  1:2108558/1‑70 (MQ=255)
   ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                               >  1:1002818/1‑71 (MQ=255)
   ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                               >  1:1282906/1‑71 (MQ=255)
      aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                             <  1:2694977/70‑1 (MQ=255)
      aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                             <  1:264757/70‑1 (MQ=255)
      aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                             <  1:2502981/70‑1 (MQ=255)
          ccTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                        >  1:145123/1‑71 (MQ=255)
           ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                        <  1:910165/70‑1 (MQ=255)
               tgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                                   >  1:3792193/1‑71 (MQ=255)
               tgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                                   >  1:1300684/1‑71 (MQ=255)
                 ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATg                                                                    >  1:2330520/1‑52 (MQ=255)
                    ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                               <  1:2198588/70‑1 (MQ=255)
                    ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                               <  1:1404730/70‑1 (MQ=255)
                    ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                               <  1:2531698/70‑1 (MQ=255)
                         tAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                                  >  1:3532970/1‑46 (MQ=255)
                         tAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt                                         >  1:941510/1‑71 (MQ=255)
                            tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                               <  1:2401492/62‑1 (MQ=255)
                              aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                    >  1:1828542/1‑71 (MQ=255)
                                cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTg                                                            <  1:2570363/45‑1 (MQ=255)
                                   gatGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTACCGGGtttt                                        >  1:1958951/3‑62 (MQ=255)
                                    ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                          >  1:2265872/1‑43 (MQ=255)
                                      gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                                 >  1:1997915/1‑50 (MQ=255)
                                      gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTc                              >  1:2861031/1‑69 (MQ=255)
                                      gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                             <  1:2574956/70‑1 (MQ=255)
                                       aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                             <  1:248729/37‑1 (MQ=255)
                                       aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCCTCGag                           <  1:589931/71‑1 (MQ=255)
                                        aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGat                                                       <  1:1339895/42‑1 (MQ=255)
                                        aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                             >  1:1936877/1‑68 (MQ=255)
                                         aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                           >  1:566372/1‑69 (MQ=255)
                                         aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                           >  1:529629/1‑69 (MQ=255)
                                           gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                         >  1:3640/1‑69 (MQ=255)
                                             aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAt                     >  1:3378394/1‑71 (MQ=255)
                                              cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                    <  1:1265563/69‑1 (MQ=12)
                                              cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                    <  1:3645971/69‑1 (MQ=12)
                                                gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                   >  1:2362056/1‑70 (MQ=255)
                                                 cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                 >  1:2963865/1‑71 (MQ=255)
                                                     gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt             <  1:3744660/71‑1 (MQ=255)
                                                     gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                             >  1:75985/1‑39 (MQ=9)
                                                     gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                             >  1:931519/1‑39 (MQ=9)
                                                     gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGtgc                                         >  1:1021577/1‑41 (MQ=9)
                                                       cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga              <  1:2729242/68‑2 (MQ=9)
                                                        agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg               <  1:263141/66‑1 (MQ=9)
                                                        aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                          <  1:871208/39‑1 (MQ=255)
                                                        aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca          <  1:385533/71‑1 (MQ=255)
                                                        aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca          >  1:396873/1‑71 (MQ=255)
                                                         taTACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt                                         >  1:454934/2‑39 (MQ=255)
                                                         gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat         <  1:803100/68‑1 (MQ=9)
                                                         gATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                   <  1:303918/61‑1 (MQ=255)
                                                          aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat         >  1:3380607/1‑70 (MQ=255)
                                                            acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc      >  1:1149639/1‑71 (MQ=255)
                                                             cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc      <  1:251882/70‑1 (MQ=255)
                                                              acGGATGGTCAGCTGCGGAAATTTTTCCGGGTTTTCCATCGCGTCGagcagc                       >  1:25087/1‑52 (MQ=255)
                                                                ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg  <  1:1812077/71‑2 (MQ=255)
                                                                ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt   >  1:2612834/1‑70 (MQ=255)
                                                                    |                                                                  
GTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTA  >  minE/607996‑608130

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