Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,412 A→G 15.2% G614G (GGT→GGC pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,4120AG15.2% 186.5 / 22.6 79G614G (GGT→GGCpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (7/5);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGT  >  Exported/608342‑608481
                                                                      |                                                                     
gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                                                                       >  1:4342200/1‑71 (MQ=255)
     gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa                                                                  <  1:1332168/71‑1 (MQ=255)
          aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa                                                             >  1:4628740/1‑71 (MQ=255)
            cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                           <  1:3273572/71‑1 (MQ=255)
            cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                           >  1:5220220/1‑71 (MQ=255)
             aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                          <  1:1479695/71‑1 (MQ=255)
              cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                          >  1:2916497/1‑70 (MQ=255)
                cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                       <  1:4386550/71‑1 (MQ=255)
                cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                       <  1:1196942/71‑1 (MQ=255)
                    ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc                                                                     >  1:321984/1‑53 (MQ=255)
                     cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                                   >  1:3633108/1‑70 (MQ=255)
                     cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                                  <  1:899285/71‑1 (MQ=255)
                       aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                       <  1:2862864/64‑1 (MQ=255)
                       aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                       <  1:5010659/64‑1 (MQ=255)
                       aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt                                                >  1:3351370/1‑71 (MQ=255)
                         cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg                                              >  1:3361145/1‑71 (MQ=255)
                           gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                             >  1:1669343/1‑70 (MQ=255)
                           gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                             >  1:2058523/1‑70 (MQ=255)
                           gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                             >  1:2981590/1‑70 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                                <  1:3983295/49‑1 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                                <  1:2340756/49‑1 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                          <  1:210484/55‑1 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                          <  1:4115035/55‑1 (MQ=255)
                              cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                                          >  1:2490241/1‑70 (MQ=255)
                               gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                        >  1:781163/1‑71 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                               <  1:1390336/44‑1 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                                  <  1:392881/57‑1 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                    >  1:1605881/1‑71 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                    >  1:2419006/1‑71 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                    >  1:1599268/1‑71 (MQ=255)
                                    cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                               <  1:3516012/43‑1 (MQ=255)
                                    cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa                                           >  1:2961006/1‑63 (MQ=255)
                                    cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGAc                                   <  1:5056468/71‑1 (MQ=255)
                                      cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                 <  1:985849/71‑1 (MQ=255)
                                      cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                 <  1:1131484/71‑1 (MQ=255)
                                       gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTATaa                                                     >  1:964193/1‑50 (MQ=255)
                                       gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                                  >  1:779657/1‑69 (MQ=255)
                                        acgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                       <  1:5523749/47‑1 (MQ=255)
                                         cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                              >  1:1411661/1‑71 (MQ=255)
                                          gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc                             <  1:2463970/71‑1 (MQ=255)
                                           acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                            <  1:998088/71‑1 (MQ=255)
                                            ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTa                                                        >  1:2208333/1‑42 (MQ=255)
                                            ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                           <  1:1327464/71‑1 (MQ=255)
                                            ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                           >  1:3422674/1‑71 (MQ=255)
                                             cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                          <  1:3416414/71‑1 (MQ=255)
                                             cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                          <  1:3959433/71‑1 (MQ=255)
                                             cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                          <  1:858281/71‑1 (MQ=255)
                                              gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                      >  1:1857874/1‑58 (MQ=255)
                                              gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                      >  1:5010318/1‑58 (MQ=255)
                                                 tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                        <  1:3858501/53‑1 (MQ=255)
                                                 tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                        <  1:3540171/53‑1 (MQ=255)
                                                  ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                       <  1:1368028/37‑1 (MQ=255)
                                                        ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                                   >  1:5545369/1‑35 (MQ=12)
                                                        ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                                 <  1:18036/37‑1 (MQ=12)
                                                        ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta                             >  1:2900815/1‑56 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGt                                                   <  1:3801679/35‑1 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                        <  1:1265460/46‑1 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg               >  1:4274323/1‑71 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt                        <  1:4533783/62‑1 (MQ=255)
                                                         tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                           <  1:4226166/58‑1 (MQ=255)
                                                         tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt                   <  1:3955736/66‑1 (MQ=255)
                                                          aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGAta                                               >  1:4473590/1‑36 (MQ=12)
                                                          aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                                 >  1:1242288/1‑35 (MQ=255)
                                                          aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                                 >  1:3691654/1‑35 (MQ=255)
                                                          aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg              >  1:115970/1‑70 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg              >  1:1062633/1‑66 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg              >  1:3634530/1‑66 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt         >  1:4710746/1‑71 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt         >  1:2879597/1‑71 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt         >  1:2541054/1‑71 (MQ=255)
                                                              gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt         >  1:2102100/1‑71 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                      <  1:38625/41‑1 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                      <  1:2378405/41‑1 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                      <  1:5667625/41‑1 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                <  1:4317281/63‑1 (MQ=255)
                                                                tttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt                                  >  1:4245046/1‑44 (MQ=12)
                                                                 ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt       <  1:4831545/70‑1 (MQ=255)
                                                                 ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg      >  1:1910693/1‑71 (MQ=255)
                                                                 ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                      >  1:3934166/1‑55 (MQ=255)
                                                                  tCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                    <  1:2128709/40‑1 (MQ=255)
                                                                      aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa    <  1:3568046/68‑1 (MQ=255)
                                                                      aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt  >  1:1920407/1‑70 (MQ=255)
                                                                      aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt  <  1:4449503/70‑1 (MQ=255)
                                                                      |                                                                     
GCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGT  >  Exported/608342‑608481

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: