Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,042 G→A 16.6% D404D (GAC→GAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,0420GA16.6% 308.1 / 44.4 133D404D (GAC→GATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (43/67);  new base A (12/10);  total (55/78)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 135 total aligned reads displayed.

ACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTC  >  Exported/608973‑609108
                                                                     |                                                                  
aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                                                                    <  1:3253912/70‑1 (MQ=255)
 ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                                                                   <  1:1329495/70‑1 (MQ=255)
     cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              >  1:822877/1‑71 (MQ=255)
     cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              <  1:4051844/71‑1 (MQ=255)
      gATCATCGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                                             >  1:3346541/1‑71 (MQ=255)
       accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                                                            <  1:2958484/69‑1 (MQ=255)
        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:947131/1‑64 (MQ=255)
        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:1007848/1‑64 (MQ=255)
        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:4006540/1‑64 (MQ=255)
        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:1557215/1‑64 (MQ=255)
        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                           <  1:2113837/71‑1 (MQ=255)
         catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                          >  1:1482523/1‑71 (MQ=255)
          atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                          <  1:439156/70‑1 (MQ=255)
           ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                        >  1:3290511/3‑71 (MQ=255)
                 ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                   <  1:2953341/70‑1 (MQ=255)
                  tACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                   >  1:683461/1‑69 (MQ=255)
                   aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                   >  1:4535093/1‑68 (MQ=255)
                    cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:2738835/71‑1 (MQ=255)
                     gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                              <  1:4350307/71‑4 (MQ=255)
                     gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  <  1:373776/51‑1 (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:553177/58‑1 (MQ=255)
                        gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                        <  1:564540/58‑1 (MQ=255)
                        gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:2007039/71‑1 (MQ=255)
                         caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:154232/66‑1 (MQ=255)
                         caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:312338/66‑1 (MQ=255)
                          aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                            >  1:1216229/1‑68 (MQ=255)
                          aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                            >  1:3710429/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)
                            gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                                       <  1:4835971/71‑1 (MQ=255)
                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       <  1:419717/71‑1 (MQ=255)
                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       <  1:1461097/71‑1 (MQ=255)
                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       >  1:4337394/1‑71 (MQ=255)
                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       <  1:1435471/71‑1 (MQ=255)
                             caaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATc                                                                  >  1:41501/1‑43 (MQ=12)
                              aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                                     <  1:4618897/71‑2 (MQ=255)
                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:1551045/71‑1 (MQ=255)
                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    >  1:2302400/1‑71 (MQ=255)
                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:5588921/71‑1 (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)
                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                >  1:5055751/1‑40 (MQ=255)
                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                >  1:4247255/1‑40 (MQ=255)
                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                >  1:1150438/1‑40 (MQ=255)
                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                           >  1:2581567/1‑45 (MQ=255)
                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                                 <  1:2751295/71‑1 (MQ=255)
                                  gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                  <  1:940828/70‑1 (MQ=255)
                                      aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                                     <  1:5126632/47‑1 (MQ=255)
                                       gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                                       <  1:2293534/44‑1 (MQ=255)
                                        tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                              >  1:3663707/4‑68 (MQ=255)
                                        tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                              >  1:1970741/4‑68 (MQ=255)
                                        tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                              >  1:2152665/4‑68 (MQ=255)
                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              <  1:4903047/36‑1 (MQ=255)
                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              <  1:3235701/36‑1 (MQ=255)
                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                          <  1:235236/40‑1 (MQ=255)
                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                              <  1:5667610/52‑1 (MQ=255)
                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:3379225/1‑71 (MQ=255)
                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:3179307/1‑71 (MQ=255)
                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:445539/1‑71 (MQ=255)
                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:705335/1‑71 (MQ=255)
                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:1334793/52‑1 (MQ=255)
                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:2456467/52‑1 (MQ=255)
                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:5640415/52‑1 (MQ=255)
                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:3110437/1‑70 (MQ=255)
                                            cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:2280545/51‑1 (MQ=255)
                                             gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                       <  1:845802/70‑1 (MQ=255)
                                             gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                      >  1:5548385/1‑71 (MQ=255)
                                               tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                    >  1:456444/1‑71 (MQ=255)
                                                 ttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                                          <  1:5177378/47‑1 (MQ=255)
                                                  tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                              <  1:2901027/42‑1 (MQ=255)
                                                   gAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                                        >  1:2996215/1‑47 (MQ=255)
                                                        ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt              <  1:1150458/68‑1 (MQ=255)
                                                         cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc           <  1:5337889/70‑1 (MQ=255)
                                                         cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc           <  1:2198931/70‑1 (MQ=255)
                                                           gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        <  1:3077275/71‑1 (MQ=255)
                                                           gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        <  1:5231204/71‑1 (MQ=255)
                                                            aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa       <  1:3909367/71‑1 (MQ=255)
                                                             cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        >  1:1036772/1‑69 (MQ=255)
                                                             cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg        >  1:5175063/1‑69 (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:4737692/71‑1 (MQ=255)
                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc  <  1:2073074/71‑1 (MQ=255)
                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc  <  1:4503201/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)
                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa      >  1:1385118/1‑65 (MQ=255)
                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa       >  1:2873580/1‑64 (MQ=255)
                                                                     |                                                                  
ACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTC  >  Exported/608973‑609108

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: