Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,822 A→G 7.6% T144T (ACT→ACC pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,8220AG 7.6% 224.3 / 6.7 79T144T (ACT→ACCpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (27/46);  new base G (3/3);  total (30/49)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAAC  >  Exported/609757‑609892
                                                                 |                                                                      
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:718323/53‑1 (MQ=255)
             tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                                        <  1:5670167/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:3635746/1‑51 (MQ=255)
                       cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                                >  1:2670501/1‑51 (MQ=255)
                       cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                                >  1:1739280/1‑51 (MQ=255)
                       cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                                                >  1:683568/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:4448383/39‑1 (MQ=255)
                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                                      <  1:1694819/39‑1 (MQ=255)
                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                                      <  1:166337/39‑1 (MQ=255)
                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                                      <  1:3298447/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)
                                                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:1416889/37‑1 (MQ=255)
                                                            ttGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                                         <  1:4693357/37‑1 (MQ=255)
                                                            ttGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                                         <  1:1054771/37‑1 (MQ=255)
                                                               tgAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAgg     <  1:1130080/70‑1 (MQ=255)
                                                               tgAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAgg     <  1:3785648/70‑1 (MQ=255)
                                                                 aGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAAc  >  1:4633532/1‑71 (MQ=255)
                                                                 aGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAAc  >  1:3035167/1‑71 (MQ=255)
                                                                 |                                                                      
AGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAAC  >  Exported/609757‑609892

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: