Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,819 G→A 7.5% H145H (CAC→CAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,8190GA 7.5% 237.0 / 6.2 80H145H (CAC→CATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (27/47);  new base A (3/3);  total (30/50)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTG  >  minE/609750‑609885
                                                                     |                                                                  
aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                                   <  1:183171/71‑1 (MQ=255)
  cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatgg                                                                 >  1:3710304/1‑70 (MQ=255)
       aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                             >  1:4540961/1‑70 (MQ=255)
         aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                           >  1:1366394/1‑70 (MQ=255)
             ccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGtt                                                              <  1:1475341/63‑1 (MQ=16)
             ccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACgcag                                                      <  1:5526446/71‑4 (MQ=255)
              cAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                                         >  1:4448415/1‑67 (MQ=255)
               aGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                                         >  1:5020873/1‑66 (MQ=255)
               aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                                                    >  1:4220959/1‑71 (MQ=255)
                gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                                                    >  1:5335648/1‑70 (MQ=255)
                gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                                                   <  1:2987674/71‑1 (MQ=255)
                 aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGt                                                               <  1:1816085/58‑1 (MQ=12)
                    tACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                                         >  1:624521/1‑61 (MQ=16)
                    tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                                 <  1:5670167/53‑1 (MQ=255)
                    tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                                 <  1:718323/53‑1 (MQ=255)
                     aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                               >  1:4481659/1‑70 (MQ=255)
                        gCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                                 <  1:257130/49‑1 (MQ=255)
                         cAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                                            >  1:1238475/1‑53 (MQ=255)
                          aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt                                         >  1:3581036/1‑71 (MQ=255)
                           aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                               >  1:2324190/1‑64 (MQ=255)
                           aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                               <  1:2607159/64‑1 (MQ=255)
                            cGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                                       <  1:1350574/71‑1 (MQ=255)
                             gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                           >  1:708863/1‑50 (MQ=255)
                             gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                           >  1:1039979/1‑50 (MQ=255)
                             gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                                      <  1:1592137/71‑1 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                         >  1:1739280/1‑51 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                         >  1:683568/1‑51 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                         >  1:3635746/1‑51 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                         >  1:2670501/1‑51 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                         >  1:2460056/1‑51 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                                     >  1:2801026/1‑71 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                                     >  1:1047864/1‑71 (MQ=255)
                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                                      <  1:5433830/70‑1 (MQ=255)
                                ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAg                                                                <  1:3102581/42‑1 (MQ=255)
                                ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCGGTGAAGATTTTTTTg                                   <  1:1513943/71‑1 (MQ=255)
                                ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                                   <  1:3189165/71‑1 (MQ=255)
                                 gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                           <  1:3490918/46‑1 (MQ=255)
                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                               <  1:166337/39‑1 (MQ=255)
                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                               <  1:4448383/39‑1 (MQ=255)
                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                               <  1:3298447/39‑1 (MQ=255)
                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                               <  1:1694819/39‑1 (MQ=255)
                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                               <  1:2950033/39‑1 (MQ=255)
                                      ccGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                                      <  1:4760888/62‑1 (MQ=255)
                                       cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                             <  1:1222441/38‑1 (MQ=255)
                                       cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                                                   <  1:4623315/48‑1 (MQ=255)
                                       cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                                      <  1:1991010/61‑1 (MQ=255)
                                         gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt                                    <  1:2310444/61‑1 (MQ=255)
                                          aGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt                                    <  1:210451/60‑1 (MQ=255)
                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGa                                              >  1:4853562/1‑46 (MQ=255)
                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatc                                >  1:4175961/1‑60 (MQ=255)
                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                                  >  1:1071532/1‑58 (MQ=255)
                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCa                      <  1:4763167/70‑1 (MQ=255)
                                               aaaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg                           <  1:64306/64‑1 (MQ=255)
                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                              <  1:1853019/59‑1 (MQ=255)
                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                              <  1:3468678/59‑1 (MQ=255)
                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc                  >  1:5094778/1‑71 (MQ=255)
                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc                  <  1:1604129/71‑1 (MQ=255)
                                                   gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                              >  1:4820410/1‑57 (MQ=255)
                                                   gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                              >  1:5603145/1‑57 (MQ=255)
                                                   gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                            >  1:3815200/1‑59 (MQ=255)
                                                   gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                            >  1:2703148/1‑59 (MQ=255)
                                                    tCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                                                   <  1:2157885/35‑1 (MQ=255)
                                                    tCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                                       <  1:1833171/47‑1 (MQ=255)
                                                    tCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                                       <  1:694077/47‑1 (MQ=255)
                                                       aaCACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGtt            <  1:2846060/71‑1 (MQ=255)
                                                        acacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCa                      >  1:1178933/1‑60 (MQ=255)
                                                         cacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                                   <  1:1811448/46‑1 (MQ=255)
                                                          acGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa                          <  1:1515411/54‑1 (MQ=255)
                                                           cGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGt             <  1:5306097/66‑1 (MQ=255)
                                                              ccTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCt     <  1:2460983/71‑1 (MQ=255)
                                                               cTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                                      <  1:3202469/37‑1 (MQ=255)
                                                               cTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa                          <  1:577646/49‑1 (MQ=255)
                                                               cTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGtt            <  1:4684239/63‑1 (MQ=255)
                                                               cTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGtt            <  1:3093267/63‑1 (MQ=255)
                                                               cTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTAc        <  1:5336736/67‑1 (MQ=255)
                                                               cTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTAc        <  1:2848391/67‑1 (MQ=255)
                                                                 ggTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatc                                <  1:1107521/41‑1 (MQ=255)
                                                                 ggTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTg  >  1:2483856/1‑71 (MQ=255)
                                                                  gTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                                   <  1:3959453/37‑1 (MQ=255)
                                                                   ttGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                                  <  1:4693357/37‑1 (MQ=255)
                                                                   ttGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                                  <  1:1416889/37‑1 (MQ=255)
                                                                   ttGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                                  <  1:1054771/37‑1 (MQ=255)
                                                                     |                                                                  
ACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTG  >  minE/609750‑609885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: