Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,037 A→G 12.2% R739R (CGT→CGC pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,0370AG12.2% 207.0 / 17.0 82R739R (CGT→CGCpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (43/29);  new base G (4/6);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

CTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCA  >  Exported/607967‑608107
                                                                       |                                                                      
cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                                        >  1:982748/1‑71 (MQ=255)
cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                                        >  1:2170650/1‑71 (MQ=255)
cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                                        >  1:4477486/1‑71 (MQ=255)
   aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                                        <  1:4419698/68‑1 (MQ=255)
   aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                                      <  1:1962975/70‑1 (MQ=255)
   aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                                      >  1:2298041/1‑70 (MQ=255)
   aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                                     <  1:3026259/71‑1 (MQ=255)
    ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                                     >  1:4396549/1‑70 (MQ=255)
    ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:2300988/1‑71 (MQ=255)
    ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                                     >  1:3494523/1‑70 (MQ=255)
       cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACt                                                                  <  1:4815956/70‑1 (MQ=255)
         tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGGCAGCGAGTTGAAACGTACTGc                                                                <  1:3749529/70‑1 (MQ=255)
          aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                                               >  1:3536577/1‑70 (MQ=255)
             ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                           <  1:2617571/71‑1 (MQ=255)
              tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                                      >  1:1801633/1‑59 (MQ=255)
              tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                                      >  1:1953076/1‑59 (MQ=255)
               gAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                                                          <  1:83408/70‑1 (MQ=255)
                 gtaaaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa                                                       <  1:350587/68‑1 (MQ=255)
                  tGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                                     <  1:1294776/56‑1 (MQ=255)
                    aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                                     <  1:2969174/70‑3 (MQ=255)
                    aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                                     <  1:404675/70‑3 (MQ=255)
                    aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:227011/1‑70 (MQ=255)
                    aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:5315911/1‑70 (MQ=255)
                     aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                                      >  1:805418/1‑52 (MQ=255)
                     aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                                      >  1:5010636/1‑52 (MQ=255)
                        tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                                                  >  1:3279435/1‑69 (MQ=255)
                        tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                 >  1:398507/1‑70 (MQ=255)
                            aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg                                            <  1:87811/71‑1 (MQ=255)
                            aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                                            <  1:3814349/71‑1 (MQ=255)
                              taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                <  1:5158368/65‑1 (MQ=255)
                              taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                <  1:3326983/65‑1 (MQ=255)
                              taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                                 <  1:366473/64‑1 (MQ=255)
                               aataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa                                         <  1:569736/71‑1 (MQ=255)
                                ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:531529/1‑65 (MQ=255)
                                ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:2715948/1‑65 (MQ=255)
                                ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:5538430/1‑65 (MQ=255)
                                ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:3593189/1‑65 (MQ=255)
                                ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                        <  1:472346/71‑1 (MQ=255)
                                ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                        >  1:1064665/1‑71 (MQ=255)
                                ataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg                                        >  1:4040379/1‑71 (MQ=255)
                                 taaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                         <  1:3450180/69‑1 (MQ=255)
                                  aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:3235411/1‑41 (MQ=255)
                                  aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:2849456/1‑41 (MQ=255)
                                  aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:739180/1‑41 (MQ=255)
                                  aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                                            <  1:1961296/49‑1 (MQ=255)
                                   aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                                              >  1:1115134/1‑46 (MQ=255)
                                   aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                >  1:2975706/1‑60 (MQ=255)
                                   aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                         >  1:2340301/1‑67 (MQ=255)
                                     gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                   >  1:3563115/1‑71 (MQ=255)
                                     gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                   >  1:5527442/1‑71 (MQ=255)
                                      tCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                           <  1:5225139/62‑1 (MQ=255)
                                       ccTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTAc                                                                   >  1:2986568/1‑37 (MQ=255)
                                        ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                         >  1:753132/1‑62 (MQ=255)
                                        ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:2966344/1‑70 (MQ=255)
                                        ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:3375138/1‑70 (MQ=255)
                                        ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:1178530/1‑70 (MQ=255)
                                        ctgctgcagTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:562899/1‑70 (MQ=255)
                                         tgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                      <  1:1232043/64‑1 (MQ=255)
                                           ctgctgGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                                >  1:4867945/1‑37 (MQ=38)
                                              ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                          <  1:1052531/39‑1 (MQ=255)
                                              ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                          <  1:1119205/71‑1 (MQ=255)
                                                  ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                        <  1:1488244/70‑1 (MQ=255)
                                                   tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                      >  1:2144533/1‑71 (MQ=255)
                                                   tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                      >  1:5385737/1‑71 (MQ=255)
                                                   tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                      >  1:3291109/1‑71 (MQ=255)
                                                    cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                                       >  1:5045808/1‑37 (MQ=255)
                                                    cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                                       >  1:4501310/1‑37 (MQ=255)
                                                    cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:2000940/1‑39 (MQ=255)
                                                    cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:2851405/1‑39 (MQ=255)
                                                    cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                        >  1:4205704/1‑52 (MQ=255)
                                                        aGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                          <  1:5332446/62‑1 (MQ=255)
                                                         gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                      >  1:3904211/1‑65 (MQ=9)
                                                         gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                >  1:4713787/1‑71 (MQ=255)
                                                         gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                <  1:2381112/71‑1 (MQ=255)
                                                          tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc               <  1:1067665/71‑1 (MQ=255)
                                                          tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc               <  1:4341197/71‑1 (MQ=255)
                                                            aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                   <  1:5478999/65‑1 (MQ=255)
                                                            aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                   <  1:862411/65‑1 (MQ=255)
                                                            aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt             >  1:4465790/1‑71 (MQ=255)
                                                             gCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                           <  1:2426693/40‑1 (MQ=255)
                                                              cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcg           >  1:2599153/1‑71 (MQ=255)
                                                                cgTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                        <  1:5500462/55‑1 (MQ=12)
                                                                  ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                         <  1:1169838/53‑1 (MQ=255)
                                                                  ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                         <  1:1855915/53‑1 (MQ=255)
                                                                    gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCa              >  1:3090896/1‑62 (MQ=255)
                                                                    gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg      <  1:1841906/70‑1 (MQ=255)
                                                                      aaCGTACTGCGTAGCCAGATACAGGGATGGTCAGCTGCGGATATTTTTcc                        <  1:4784597/50‑1 (MQ=255)
                                                                       aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca  <  1:2574704/71‑1 (MQ=255)
                                                                       aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag    >  1:4723853/1‑69 (MQ=255)
                                                                       |                                                                      
CTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCA  >  Exported/607967‑608107

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: