Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,118 G→A 8.2% N712N (AAC→AAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,1180GA 8.2% 234.1 / 7.5 73N712N (AAC→AATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (3/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCA  >  minE/608049‑608188
                                                                     |                                                                      
gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                        <  1:4059636/70‑1 (MQ=255)
gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                        <  1:3759950/70‑1 (MQ=255)
gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                        <  1:2782407/70‑1 (MQ=255)
gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                        <  1:1554197/70‑1 (MQ=255)
gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                                       >  1:2263107/1‑71 (MQ=255)
 ccAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                                                      >  1:678018/1‑71 (MQ=255)
      tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatca                                                                 >  1:2134546/1‑71 (MQ=255)
        cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                <  1:954018/70‑1 (MQ=255)
          cGGATGGTCAGCT‑‑GGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                                       >  1:3782096/1‑59 (MQ=38)
           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)
                   aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                        >  1:873134/1‑51 (MQ=255)
                   aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                     <  1:428259/70‑1 (MQ=255)
                       taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                                >  1:2042885/3‑70 (MQ=255)
                        cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                                   <  1:3068170/51‑1 (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)
                          gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                                                      <  1:3237633/46‑1 (MQ=255)
                          gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                             >  1:87146/1‑71 (MQ=255)
                          gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGAtt                                                                      <  1:2344523/46‑1 (MQ=12)
                           atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                       <  1:2751272/60‑1 (MQ=255)
                           atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                            >  1:2458871/1‑71 (MQ=255)
                              tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                                   <  1:2823193/45‑1 (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:1257222/70‑1 (MQ=255)
                                   ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                     <  1:2698739/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:757415/71‑1 (MQ=255)
                                      ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                 <  1:662534/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:2926582/37‑1 (MQ=255)
                                         tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                <  1:2931425/37‑1 (MQ=255)
                                         tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                <  1:1306678/37‑1 (MQ=255)
                                          ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                             >  1:1585647/1‑71 (MQ=255)
                                          ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                             >  1:319676/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:1333066/1‑70 (MQ=255)
                                                      gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                  >  1:3085286/1‑70 (MQ=255)
                                                          agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa               >  1:591115/1‑68 (MQ=12)
                                                           gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                               >  1:856523/1‑52 (MQ=255)
                                                           gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                               >  1:58018/1‑52 (MQ=255)
                                                           gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                               >  1:65986/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)
                                                                    ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                               >  1:581999/1‑43 (MQ=255)
                                                                    ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    >  1:1601783/1‑70 (MQ=255)
                                                                    ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    >  1:2027532/1‑70 (MQ=255)
                                                                    ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    >  1:1273188/1‑70 (MQ=255)
                                                                    ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat    >  1:354624/1‑70 (MQ=255)
                                                                     gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa  <  1:585938/71‑1 (MQ=255)
                                                                     gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa  <  1:1962647/71‑1 (MQ=255)
                                                                     |                                                                      
GCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCA  >  minE/608049‑608188

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: