Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,807 C→T 17.5% V149V (GTG→GTA pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,8070CT17.5% 174.6 / 23.5 74V149V (GTG→GTApflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base T (9/4);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTT  >  Exported/609738‑609875
                                                                     |                                                                    
tGCATCTGGCAGACCGGGCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                     >  1:3822322/1‑71 (MQ=255)
     cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                >  1:4247352/1‑71 (MQ=255)
      tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                                               >  1:5075076/1‑71 (MQ=255)
      tGGCAGACCGGTCAGAACACCAGATTTACGGAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                                              >  1:2475998/1‑71 (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:4575174/1‑60 (MQ=12)
                gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                                                >  1:4581481/1‑60 (MQ=12)
                gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                     >  1:1149188/1‑55 (MQ=255)
                gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAac                                                                      >  1:1333273/1‑54 (MQ=255)
                   aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                   >  1:4540961/1‑70 (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)
                                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:708863/1‑50 (MQ=255)
                                         gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                 >  1:1039979/1‑50 (MQ=255)
                                         gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                            <  1:1592137/71‑1 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                               >  1:683568/1‑51 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                               >  1:3635746/1‑51 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                               >  1:2670501/1‑51 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                               >  1:2460056/1‑51 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                               >  1:1739280/1‑51 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                           >  1:1047864/1‑71 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                           >  1:2801026/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:3298447/39‑1 (MQ=255)
                                                gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                     <  1:166337/39‑1 (MQ=255)
                                                gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGt                                                     <  1:1694819/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:5094778/1‑71 (MQ=255)
                                                             aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc        <  1:1604129/71‑1 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                    >  1:5603145/1‑57 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                    >  1:4820410/1‑57 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                  >  1:2703148/1‑59 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                  >  1:3815200/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)
                                                                     |                                                                    
TGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTT  >  Exported/609738‑609875

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: