Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,045 A→G 16.3% D403D (GAT→GAC pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,0450AG16.3% 294.0 / 43.7 129D403D (GAT→GACpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (41/67);  new base G (11/10);  total (52/77)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 131 total aligned reads displayed.

AACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGC  >  minE/608976‑609114
                                                                     |                                                                     
aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                      <  1:2096741/71‑1 (MQ=255)
  cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                    >  1:822877/1‑71 (MQ=255)
  cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                    <  1:4051844/71‑1 (MQ=255)
    accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                                                                  <  1:2958484/69‑1 (MQ=255)
      catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                >  1:4750025/1‑71 (MQ=255)
      catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                >  1:5542797/1‑71 (MQ=255)
       atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                <  1:439156/70‑1 (MQ=255)
        ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                              >  1:3290511/3‑71 (MQ=255)
        ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                              >  1:4556021/3‑71 (MQ=255)
        ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                              >  1:5684195/3‑71 (MQ=255)
          gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                                             <  1:4952891/70‑1 (MQ=255)
          gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                                            >  1:230515/1‑71 (MQ=255)
              ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                         <  1:2953341/70‑1 (MQ=255)
                aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                         >  1:4535093/1‑68 (MQ=255)
                 cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                     <  1:2738835/71‑1 (MQ=255)
                  gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                                    <  1:4350307/71‑4 (MQ=255)
                   cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                                                   <  1:2118396/71‑1 (MQ=255)
                    agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGATCTCATACTGCAGagaa                                                  <  1:1570117/71‑1 (MQ=255)
                     gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag                                                 <  1:2366258/71‑1 (MQ=255)
                     gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                              <  1:564540/58‑1 (MQ=255)
                     gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                              <  1:553177/58‑1 (MQ=255)
                     gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                              <  1:1710426/58‑1 (MQ=255)
                     gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                              <  1:3215043/58‑1 (MQ=255)
                      caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                     <  1:154232/66‑1 (MQ=255)
                       aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                  >  1:515667/1‑68 (MQ=255)
                       aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                  >  1:3710429/1‑68 (MQ=255)
                       aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                  >  1:1216229/1‑68 (MQ=255)
                       aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                               <  1:316184/71‑1 (MQ=255)
                       aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                               >  1:3078264/1‑71 (MQ=255)
                        ggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                              >  1:905889/2‑55 (MQ=255)
                        agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                               <  1:3330683/70‑1 (MQ=255)
                         gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                             <  1:419717/71‑1 (MQ=255)
                         gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                             <  1:1461097/71‑1 (MQ=255)
                         gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                             <  1:1435471/71‑1 (MQ=255)
                           aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                                           <  1:4618897/71‑2 (MQ=255)
                            aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCa                                                                       <  1:4323708/42‑1 (MQ=255)
                            aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                          <  1:1082422/71‑1 (MQ=255)
                            aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                          <  1:1551045/71‑1 (MQ=255)
                             tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                              <  1:1385924/50‑1 (MQ=255)
                             tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                     <  1:1820379/59‑1 (MQ=255)
                             tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                         <  1:5383566/71‑1 (MQ=255)
                               gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                      >  1:5055751/1‑40 (MQ=255)
                               gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                      >  1:1150438/1‑40 (MQ=255)
                               gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                      >  1:4247255/1‑40 (MQ=255)
                               gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGtt                                                                  >  1:938943/1‑44 (MQ=255)
                               gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                                 >  1:2581567/1‑45 (MQ=255)
                               gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                        <  1:940828/70‑1 (MQ=255)
                               gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt                                       >  1:2137280/1‑71 (MQ=255)
                                  tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                                                   <  1:1850512/40‑1 (MQ=255)
                                  tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa                                    <  1:3301044/71‑1 (MQ=255)
                                   aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                                           <  1:2223316/47‑1 (MQ=255)
                                    gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                                             <  1:2293534/44‑1 (MQ=255)
                                     tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                 <  1:1115463/68‑1 (MQ=255)
                                     tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                    >  1:1970741/4‑68 (MQ=255)
                                     tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                    <  1:5290414/36‑1 (MQ=255)
                                     tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                    <  1:4903047/36‑1 (MQ=255)
                                     tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                <  1:235236/40‑1 (MQ=255)
                                     tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                <  1:3550239/40‑1 (MQ=255)
                                       aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                                            >  1:770698/1‑42 (MQ=255)
                                       aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa                                                       <  1:4137198/47‑1 (MQ=255)
                                       aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:3379225/1‑71 (MQ=255)
                                       aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:3179307/1‑71 (MQ=255)
                                       aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:705335/1‑71 (MQ=255)
                                        tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                 <  1:5640415/52‑1 (MQ=255)
                                        tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                 <  1:1334793/52‑1 (MQ=255)
                                        tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:3110437/1‑70 (MQ=255)
                                        tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:3116851/1‑70 (MQ=255)
                                         cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                 <  1:2280545/51‑1 (MQ=255)
                                         cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:2781174/1‑69 (MQ=255)
                                          gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                            >  1:5548385/1‑71 (MQ=255)
                                              ttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                                                <  1:5177378/47‑1 (MQ=255)
                                               tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                    <  1:4411910/42‑1 (MQ=255)
                                               tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                    <  1:2901027/42‑1 (MQ=255)
                                                gAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                                              >  1:2996215/1‑47 (MQ=255)
                                                      cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                 <  1:2198931/70‑1 (MQ=255)
                                                      cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                 <  1:5337889/70‑1 (MQ=255)
                                                        gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg              <  1:3077275/71‑1 (MQ=255)
                                                        gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGAAACTTTAGCGGCGAATTTCtt               >  1:249047/1‑70 (MQ=255)
                                                         aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa             <  1:3909367/71‑1 (MQ=255)
                                                          cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg              >  1:1036772/1‑69 (MQ=255)
                                                          cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg              >  1:5175063/1‑69 (MQ=255)
                                                          cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGa             <  1:1731496/70‑1 (MQ=255)
                                                            cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                          >  1:1500114/1‑55 (MQ=255)
                                                            cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg          <  1:1224327/71‑2 (MQ=255)
                                                             aTCAGATCGTCGTTTTCATACTGCAAAGAGGATGTGACGATAGACACCTGCGCGGCATATTTTTTGaa            >  1:3717520/1‑68 (MQ=255)
                                                              tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc        <  1:3847743/71‑1 (MQ=255)
                                                              tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                          <  1:611434/37‑1 (MQ=255)
                                                              tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                          <  1:16782/37‑1 (MQ=255)
                                                              tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg              <  1:2152486/65‑1 (MQ=255)
                                                              tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc        <  1:2073074/71‑1 (MQ=255)
                                                              tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc        <  1:4503201/71‑1 (MQ=255)
                                                              tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc        <  1:4737692/71‑1 (MQ=255)
                                                               cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt         >  1:4818782/1‑69 (MQ=255)
                                                                aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                          <  1:3710638/35‑1 (MQ=255)
                                                                aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                        >  1:5544640/1‑37 (MQ=255)
                                                                aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa            >  1:2111126/1‑65 (MQ=255)
                                                                aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa             >  1:2873580/1‑64 (MQ=255)
                                                                  gTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                             <  1:5636241/46‑1 (MQ=255)
                                                                    cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                               >  1:3452629/1‑42 (MQ=255)
                                                                    cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg   >  1:5694392/1‑70 (MQ=255)
                                                                     |                                                                     
AACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGC  >  minE/608976‑609114

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