Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,048 G→A 16.1% N402N (AAC→AAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,0480GA16.1% 292.5 / 40.0 124N402N (AAC→AATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (38/66);  new base A (10/10);  total (48/76)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 131 total aligned reads displayed.

CGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTT  >  minE/608978‑609118
                                                                      |                                                                      
cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                        >  1:822877/1‑71 (MQ=255)
 gATCATCGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                                                       >  1:3346541/1‑71 (MQ=255)
  accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                                                                      <  1:2958484/69‑1 (MQ=255)
    catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                    >  1:5542797/1‑71 (MQ=255)
    catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                    >  1:1482523/1‑71 (MQ=255)
    catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                    >  1:4750025/1‑71 (MQ=255)
     atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                    <  1:439156/70‑1 (MQ=255)
      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                  >  1:5684195/3‑71 (MQ=255)
      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                  >  1:3290511/3‑71 (MQ=255)
      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                  >  1:4556021/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)
                gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                                        <  1:3895082/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:1710426/58‑1 (MQ=255)
                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                  <  1:564540/58‑1 (MQ=255)
                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                  <  1:553177/58‑1 (MQ=255)
                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                  <  1:3215043/58‑1 (MQ=255)
                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                     <  1:2007039/71‑1 (MQ=255)
                    caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                         <  1:154232/66‑1 (MQ=255)
                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                      >  1:3710429/1‑68 (MQ=255)
                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                      >  1:515667/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)
                       gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                                                 <  1:4835971/71‑1 (MQ=255)
                       gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                 <  1:1435471/71‑1 (MQ=255)
                       gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                 <  1:419717/71‑1 (MQ=255)
                         aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                                               <  1:4618897/71‑2 (MQ=255)
                          aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                              <  1:1082422/71‑1 (MQ=255)
                           tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                         <  1:1803017/59‑1 (MQ=255)
                           tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                         <  1:1820379/59‑1 (MQ=255)
                           tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                             <  1:5383566/71‑1 (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:3663707/4‑68 (MQ=255)
                                   tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                        <  1:5290414/36‑1 (MQ=255)
                                   tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                        <  1:3235701/36‑1 (MQ=255)
                                   tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                                        <  1:4903047/36‑1 (MQ=255)
                                   tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                    <  1:235236/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:2901027/42‑1 (MQ=255)
                                             tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                        <  1:4411910/42‑1 (MQ=255)
                                              gAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                                                  >  1:2996215/1‑47 (MQ=255)
                                                   ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt                        <  1:1150458/68‑1 (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)
                                                        cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                     >  1:1559307/1‑50 (MQ=12)
                                                        cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg                  >  1:5175063/1‑69 (MQ=255)
                                                          cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                              >  1:1500114/1‑55 (MQ=255)
                                                          cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg              <  1:1224327/71‑2 (MQ=255)
                                                            tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc            <  1:3847743/71‑1 (MQ=255)
                                                            tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                              <  1:16782/37‑1 (MQ=255)
                                                            tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                              <  1:611434/37‑1 (MQ=255)
                                                            tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                  <  1:2152486/65‑1 (MQ=255)
                                                            tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc            <  1:4503201/71‑1 (MQ=255)
                                                            tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc            <  1:4737692/71‑1 (MQ=255)
                                                            tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc            <  1:2073074/71‑1 (MQ=255)
                                                              aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                              <  1:3710638/35‑1 (MQ=255)
                                                              aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                            >  1:5544640/1‑37 (MQ=255)
                                                              aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                >  1:1385118/1‑65 (MQ=255)
                                                              aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                >  1:2111126/1‑65 (MQ=255)
                                                              aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                 >  1:2873580/1‑64 (MQ=255)
                                                                gTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                 <  1:5636241/46‑1 (MQ=255)
                                                                 tCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg       <  1:4179864/71‑1 (MQ=255)
                                                                  cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                   >  1:3452629/1‑42 (MQ=255)
                                                                  cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg       >  1:5694392/1‑70 (MQ=255)
                                                                    tCGTTCTCATTCTGCAGAGAAGAGGTGTCGATGGACACTTTAGGGGCGAATttgttg                  <  1:3695301/57‑1 (MQ=255)
                                                                    tCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAg    <  1:1416829/71‑1 (MQ=255)
                                                                    tCGTTCTCATACTGCAGAGAAGAGCTGTCGATGGa                                        <  1:4657500/35‑1 (MQ=255)
                                                                      gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGGGGCAGtt  >  1:145938/1‑71 (MQ=255)
                                                                      |                                                                      
CGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTT  >  minE/608978‑609118

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: