Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,124 C→T 7.0% V710V (GTG→GTA pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,1240CT 7.0% 233.7 / 4.8 71V710V (GTG→GTApflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base T (3/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAG  >  minE/608055‑608194
                                                                     |                                                                      
tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatca                                                                       >  1:2134546/1‑71 (MQ=255)
  cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                      <  1:954018/70‑1 (MQ=255)
     ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTa                                                                  <  1:2867888/71‑1 (MQ=255)
     ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTa                                                                  <  1:2862875/71‑1 (MQ=255)
      gATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttga                                                                 <  1:2025127/71‑3 (MQ=255)
      gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                                 <  1:3732417/71‑1 (MQ=255)
          gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                                        >  1:2968005/1‑60 (MQ=255)
          gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                             <  1:1895377/71‑1 (MQ=255)
           tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                            <  1:2209672/71‑1 (MQ=255)
            cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                             >  1:341219/1‑69 (MQ=255)
             aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                           <  1:428259/70‑1 (MQ=255)
                 taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                                      >  1:2042885/3‑70 (MQ=255)
                  cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACg                                                               <  1:1592085/61‑1 (MQ=255)
                  cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                     <  1:2993827/71‑1 (MQ=255)
                  aggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgtt                                                     >  1:1240154/2‑69 (MQ=255)
                   ggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTg                                                    <  1:808871/71‑1 (MQ=12)
                    gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                   >  1:87146/1‑71 (MQ=255)
                     atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                             <  1:2751272/60‑1 (MQ=255)
                     atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                                  >  1:2458871/1‑71 (MQ=255)
                          tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                           >  1:1845417/1‑57 (MQ=255)
                          tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                     >  1:3168975/1‑63 (MQ=255)
                          tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                              >  1:747842/1‑70 (MQ=255)
                             ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                           <  1:2698739/70‑1 (MQ=255)
                             ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                           <  1:1257222/70‑1 (MQ=255)
                              cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                       <  1:2694289/57‑1 (MQ=255)
                              cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacc                                                 <  1:1602284/63‑1 (MQ=255)
                                ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                       <  1:2205812/71‑1 (MQ=255)
                                ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                       <  1:662534/71‑1 (MQ=255)
                                ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                       <  1:757415/71‑1 (MQ=255)
                                 gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                             <  1:178005/48‑1 (MQ=255)
                                 gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                             <  1:3973057/48‑1 (MQ=255)
                                 gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                       >  1:1728911/1‑70 (MQ=255)
                                  ttttCCATCGCGTCGAGCATCATTTCACGGTTCATCACGTTAACGtt                                                             <  1:3039242/47‑1 (MQ=255)
                                  ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                            <  1:2731694/48‑1 (MQ=255)
                                   tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                      <  1:1306678/37‑1 (MQ=255)
                                   tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                      <  1:2931425/37‑1 (MQ=255)
                                   tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                      <  1:2926582/37‑1 (MQ=255)
                                    ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                   >  1:319676/1‑71 (MQ=255)
                                    ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                   >  1:1585647/1‑71 (MQ=255)
                                    ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                    >  1:3293199/1‑70 (MQ=255)
                                    ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                       <  1:1791464/67‑1 (MQ=255)
                                     tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                  >  1:994847/1‑71 (MQ=255)
                                      ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                 >  1:1545759/1‑71 (MQ=255)
                                         tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                       <  1:637657/46‑1 (MQ=255)
                                           gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                         >  1:1281253/1‑58 (MQ=255)
                                           gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                            >  1:789870/1‑71 (MQ=255)
                                             gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                          >  1:3071801/1‑71 (MQ=255)
                                              tCGATCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                           <  1:1402645/37‑1 (MQ=255)
                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                     >  1:2311773/1‑57 (MQ=255)
                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                            >  1:1126238/1‑66 (MQ=255)
                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                        >  1:3085286/1‑70 (MQ=255)
                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                        >  1:1333066/1‑70 (MQ=255)
                                                    agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                     >  1:591115/1‑68 (MQ=12)
                                                     gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                     >  1:856523/1‑52 (MQ=255)
                                                     gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                     >  1:65986/1‑52 (MQ=255)
                                                     gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                     >  1:58018/1‑52 (MQ=255)
                                                     gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                  <  1:91753/71‑1 (MQ=255)
                                                        tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                           <  1:2349467/43‑1 (MQ=255)
                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                                >  1:1106761/1‑36 (MQ=255)
                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                                >  1:511709/1‑36 (MQ=255)
                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                   >  1:3355235/1‑49 (MQ=255)
                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                   >  1:3546612/1‑49 (MQ=255)
                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                          >  1:1086844/1‑58 (MQ=255)
                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc             >  1:1453375/1‑71 (MQ=255)
                                                           cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc             >  1:924633/1‑70 (MQ=255)
                                                             cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat          <  1:2671877/71‑1 (MQ=255)
                                                              ggctcatcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat          <  1:3003433/67‑1 (MQ=255)
                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                     >  1:581999/1‑43 (MQ=255)
                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat          >  1:2027532/1‑70 (MQ=255)
                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat          >  1:1273188/1‑70 (MQ=255)
                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat          >  1:1601783/1‑70 (MQ=255)
                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat          >  1:354624/1‑70 (MQ=255)
                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa        <  1:1962647/71‑1 (MQ=255)
                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa        <  1:585938/71‑1 (MQ=255)
                                                                 tcatcaCGTTAACGTTCAGGTGCTGACCACCTTCgatg                                       >  1:2193590/1‑38 (MQ=255)
                                                                 tcatcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc             >  1:1791889/1‑64 (MQ=255)
                                                                 tcatcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCCTCGTGGTGGAAGTAAc                 >  1:3146625/1‑60 (MQ=255)
                                                                  catcaCGTTAAAGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGa      <  1:2986580/70‑1 (MQ=255)
                                                                     caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                      >  1:1169785/1‑51 (MQ=255)
                                                                     caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                >  1:400568/1‑57 (MQ=255)
                                                                     caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag  >  1:4059787/1‑71 (MQ=255)
                                                                     caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag  >  1:2253541/1‑71 (MQ=255)
                                                                     |                                                                      
TACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAG  >  minE/608055‑608194

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: