Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,804 G→A 14.4% F150F (TTC→TTT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,8040GA14.4% 229.5 / 23.6 90F150F (TTC→TTTpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (43/34);  new base A (9/4);  total (52/38)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTC  >  minE/609734‑609869
                                                                      |                                                                 
caTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                                   <  1:1392433/71‑1 (MQ=255)
  tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 <  1:1112119/71‑1 (MQ=255)
  tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 <  1:4442888/71‑1 (MQ=255)
   agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                                 >  1:2467907/4‑70 (MQ=255)
   agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                                 >  1:4223194/4‑70 (MQ=255)
   agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                                 >  1:4424970/4‑70 (MQ=255)
    tGCATCTGGCAGACCGGGCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                               >  1:3822322/1‑71 (MQ=255)
    gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaat                                                                <  1:4228793/68‑2 (MQ=255)
         cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                          >  1:4247352/1‑71 (MQ=255)
          tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                                   >  1:3039087/1‑61 (MQ=255)
          tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                                         >  1:5075076/1‑71 (MQ=255)
          tGGCAGACCGGTCAGAACACCAGATTTACGGAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                                        >  1:2475998/1‑71 (MQ=255)
             cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 >  1:3868989/1‑60 (MQ=255)
             cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 >  1:1376146/1‑60 (MQ=255)
             cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 >  1:4702682/1‑60 (MQ=255)
             cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCGGAGTGTAAACGTCGaa                                                                 >  1:4552518/1‑59 (MQ=255)
              aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                                      >  1:1024750/3‑70 (MQ=255)
                aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                   <  1:183171/71‑1 (MQ=255)
                  cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGcc                                                           <  1:870396/61‑1 (MQ=16)
                  cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatgg                                                 >  1:3710304/1‑70 (MQ=255)
                    gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                                          >  1:4106300/1‑60 (MQ=12)
                    gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                                          >  1:5571530/1‑60 (MQ=12)
                    gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                                          >  1:4581481/1‑60 (MQ=12)
                    gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                                          >  1:4575174/1‑60 (MQ=12)
                    gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                               >  1:1149188/1‑55 (MQ=255)
                    gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAac                                                                >  1:1333273/1‑54 (MQ=255)
                      cAGAACACCAGATTTACGGCAACGCATGATGTCCGGAGTGTAAACGTCg                                                                   <  1:417847/49‑1 (MQ=255)
                       aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                             >  1:4540961/1‑70 (MQ=255)
                        gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTTAACGTCGaa                                                                 >  1:4909086/1‑49 (MQ=255)
                        gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 >  1:2810451/1‑49 (MQ=255)
                        gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 >  1:1354427/1‑49 (MQ=255)
                         aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                               >  1:3661589/1‑50 (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)
                                gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                               >  1:425310/1‑43 (MQ=255)
                                gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                               >  1:3553148/1‑43 (MQ=255)
                                gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                                    >  1:5335648/1‑70 (MQ=255)
                                gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                                   <  1:2987674/71‑1 (MQ=255)
                                 aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGt                                               <  1:1816085/58‑1 (MQ=12)
                                  tttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                  >  1:5605969/1‑38 (MQ=255)
                                  tttAAGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                 >  1:5282951/1‑39 (MQ=255)
                                    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:1039979/1‑50 (MQ=255)
                                             gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                           >  1:708863/1‑50 (MQ=255)
                                             gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                      <  1:1592137/71‑1 (MQ=255)
                                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                         >  1:3635746/1‑51 (MQ=255)
                                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                         >  1:683568/1‑51 (MQ=255)
                                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                         >  1:1739280/1‑51 (MQ=255)
                                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                         >  1:2460056/1‑51 (MQ=255)
                                              cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                         >  1:2670501/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:1694819/39‑1 (MQ=255)
                                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                               <  1:3298447/39‑1 (MQ=255)
                                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                               <  1:166337/39‑1 (MQ=255)
                                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                               <  1:4448383/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:3468678/59‑1 (MQ=255)
                                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat              <  1:1853019/59‑1 (MQ=255)
                                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc  <  1:1604129/71‑1 (MQ=255)
                                                                 aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc  >  1:5094778/1‑71 (MQ=255)
                                                                   gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat              >  1:5603145/1‑57 (MQ=255)
                                                                   gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat              >  1:4820410/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:694077/47‑1 (MQ=255)
                                                                    tCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                       <  1:1833171/47‑1 (MQ=255)
                                                                      |                                                                 
CATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTC  >  minE/609734‑609869

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: