Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,160 T→C 22.3% A698A (GCA→GCG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,1600TC22.3% 244.8 / 62.3 121A698A (GCA→GCGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (60/34);  new base C (16/11);  total (76/45)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 136 total aligned reads displayed.

CCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  Exported/608093‑608229
                                                                   |                                                                     
ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:1309634/1‑71 (MQ=255)
ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:2694238/1‑71 (MQ=255)
     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                                                               >  1:5125089/1‑71 (MQ=255)
     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgt                                                                <  1:4659684/70‑1 (MQ=255)
       gTCGAGAAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                                             <  1:2142944/71‑1 (MQ=255)
         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                                                            >  1:1158104/1‑70 (MQ=255)
         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                                                            >  1:229330/1‑70 (MQ=255)
          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:1019380/70‑1 (MQ=255)
          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2140849/70‑1 (MQ=255)
          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2369404/70‑1 (MQ=255)
          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2385207/70‑1 (MQ=255)
          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:3282078/70‑1 (MQ=255)
          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:4259562/70‑1 (MQ=255)
          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                                                           >  1:5280240/2‑70 (MQ=255)
          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                                                           <  1:4463507/69‑1 (MQ=255)
          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                                                          <  1:2301901/70‑1 (MQ=12)
            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2948280/1‑70 (MQ=12)
            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2978923/1‑70 (MQ=12)
             cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                                                                 <  1:3371798/61‑1 (MQ=255)
              agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATa                                                      >  1:4622590/1‑71 (MQ=12)
               gcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATac                                                     >  1:5666752/1‑70 (MQ=12)
                caTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                                                    <  1:3871174/71‑3 (MQ=9)
                caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                    >  1:1334170/1‑71 (MQ=255)
                 aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                                                    >  1:1072629/1‑68 (MQ=9)
                  tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                  >  1:4915494/1‑71 (MQ=255)
                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                   >  1:4274676/1‑70 (MQ=255)
                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                   >  1:4281342/1‑70 (MQ=255)
                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAaca                                                  >  1:3366388/1‑71 (MQ=255)
                   ttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2831214/1‑63 (MQ=255)
                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                                                >  1:701406/1‑71 (MQ=255)
                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat                                                 <  1:700995/70‑1 (MQ=255)
                      aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc                                               <  1:2767935/70‑1 (MQ=255)
                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:3110125/71‑1 (MQ=255)
                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:3082681/71‑1 (MQ=255)
                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:4829242/71‑1 (MQ=255)
                        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:3961175/1‑70 (MQ=255)
                        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:5401046/1‑70 (MQ=255)
                        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:5186034/1‑70 (MQ=255)
                            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                        >  1:1950331/1‑69 (MQ=255)
                            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                        >  1:5406626/1‑69 (MQ=255)
                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                        >  1:5157747/1‑71 (MQ=255)
                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                        >  1:4958019/1‑71 (MQ=255)
                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:5129726/64‑1 (MQ=9)
                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:2189330/64‑1 (MQ=9)
                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                      <  1:5048742/69‑2 (MQ=255)
                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                      <  1:3477946/69‑2 (MQ=255)
                              tcaCGTTAACGTTCAGGTCCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAg                                          <  1:244548/67‑1 (MQ=255)
                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                                         >  1:3474597/1‑50 (MQ=255)
                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:3737946/1‑61 (MQ=255)
                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:1697490/1‑61 (MQ=255)
                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:1947077/1‑61 (MQ=255)
                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                     <  1:1231279/70‑1 (MQ=255)
                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:4196643/47‑1 (MQ=255)
                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                                    >  1:1497237/1‑70 (MQ=255)
                                     gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:847201/2‑59 (MQ=9)
                                     gaCGTTGAGGTGTTGACAGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:2651620/2‑59 (MQ=255)
                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcc                                       >  1:1661064/1‑63 (MQ=255)
                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg                                <  1:2450950/70‑1 (MQ=255)
                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:3673025/1‑71 (MQ=255)
                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:3487062/1‑71 (MQ=255)
                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:2440528/1‑71 (MQ=255)
                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:4374389/1‑71 (MQ=255)
                                       cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGATCAGCCAGGTTg                             >  1:1063503/1‑71 (MQ=255)
                                        gTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTgg                            >  1:1924003/1‑71 (MQ=255)
                                          tgaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtt                          >  1:3347066/3‑70 (MQ=255)
                                           cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:4177110/1‑71 (MQ=255)
                                              gtgtTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                                             <  1:5696337/48‑2 (MQ=255)
                                              gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      <  1:1212158/71‑1 (MQ=255)
                                              gTGATGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                                       >  1:4922195/1‑38 (MQ=255)
                                               tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                  >  1:4328242/1‑42 (MQ=255)
                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      <  1:4439662/69‑1 (MQ=255)
                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:2530585/1‑70 (MQ=255)
                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:5183963/1‑70 (MQ=255)
                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:1493020/1‑70 (MQ=255)
                                                  tGACCGCCTTCGACATACGCTTCGTGGTGGAAATACACATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                  >  1:1881386/1‑71 (MQ=255)
                                                   gACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     <  1:4104449/67‑1 (MQ=255)
                                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:935576/1‑66 (MQ=255)
                                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:3098408/1‑66 (MQ=255)
                                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                 <  1:1572301/70‑1 (MQ=255)
                                                    accaccTACGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt                >  1:4563970/1‑71 (MQ=255)
                                                        ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa                    >  1:1552807/1‑63 (MQ=255)
                                                        ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCAAGGTTGGTCTTACGa                    >  1:3175673/1‑63 (MQ=255)
                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:5352547/1‑50 (MQ=255)
                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:2331445/1‑50 (MQ=255)
                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:5594752/1‑50 (MQ=255)
                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:4045519/1‑50 (MQ=255)
                                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt          <  1:4803447/69‑1 (MQ=255)
                                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt        >  1:4551955/1‑71 (MQ=255)
                                                              atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt      <  1:1041323/71‑1 (MQ=255)
                                                              atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt       >  1:1830649/1‑70 (MQ=255)
                                                               tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt      <  1:5088166/70‑1 (MQ=255)
                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                           <  1:571495/47‑1 (MQ=255)
                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   <  1:3715514/71‑1 (MQ=255)
                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   >  1:4544237/1‑71 (MQ=255)
                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   >  1:756104/1‑71 (MQ=255)
                                                                  atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:3626318/71‑1 (MQ=255)
                                                                  atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:260989/71‑1 (MQ=255)
                                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:3246150/1‑47 (MQ=255)
                                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:3721636/1‑47 (MQ=255)
                                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:4405462/1‑47 (MQ=255)
                                                                   |                                                                     
CCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  Exported/608093‑608229

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: