Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,940,135 G→C 6.6% G258A (GGC→GCC)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,940,1350GC 6.6% 273.2 / 6.6 91G258A (GGC→GCC) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (33/52);  new base C (3/3);  total (36/55)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

CAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATC  >  AP009048/1940066‑1940205
                                                                     |                                                                      
cAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGc                                                                       <  1:3307900/71‑1 (MQ=255)
cAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGc                                                                       <  1:5567976/71‑1 (MQ=255)
cAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGc                                                                       <  1:1600550/71‑1 (MQ=255)
cAATGGATGACATCATCATCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTgg                                                                        <  1:339122/70‑1 (MQ=255)
   tGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGc                                                                       >  1:78518/1‑68 (MQ=255)
   tGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGc                                                                       >  1:3765713/1‑68 (MQ=255)
      aTGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCg                                                                 >  1:2749911/1‑71 (MQ=255)
      aTGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCg                                                                 >  1:1804613/1‑71 (MQ=255)
       tGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCgg                                                                <  1:5690348/71‑1 (MQ=255)
       tGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCgg                                                                <  1:5287549/71‑1 (MQ=255)
       tGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCgg                                                                <  1:177840/71‑1 (MQ=255)
        gACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCgg                                                                <  1:277240/70‑1 (MQ=255)
         aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGa                                                                     <  1:1013543/64‑1 (MQ=255)
         aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGa                                                                     <  1:4595406/64‑1 (MQ=255)
         aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGa                                                                     <  1:5813484/64‑1 (MQ=255)
         aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGaa                                                              <  1:4470842/71‑1 (MQ=255)
         aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGaa                                                              <  1:4036188/71‑1 (MQ=255)
              atcCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGt                                                         >  1:6439065/1‑71 (MQ=255)
               tcCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTc                                                        >  1:5010954/1‑71 (MQ=255)
                  tCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCgg                                                      >  1:5035837/1‑70 (MQ=255)
                     ccTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAtt                                                  <  1:1662652/71‑1 (MQ=255)
                     ccTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAt                                                   <  1:279609/70‑1 (MQ=255)
                      ctctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGaa                                                              >  1:6164115/1‑58 (MQ=255)
                      ctctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGaa                                                              >  1:3258423/1‑58 (MQ=255)
                      ctctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGaa                                                              >  1:1179728/1‑58 (MQ=255)
                      ctctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGaa                                                              >  1:272658/1‑58 (MQ=255)
                       tctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGCCg                                                                      <  1:1248392/49‑1 (MQ=255)
                       tctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGCCg                                                                      <  1:221573/49‑1 (MQ=255)
                       tctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGCCg                                                                      <  1:3354926/49‑1 (MQ=255)
                        ctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAtt                                                  >  1:1702718/1‑68 (MQ=255)
                         tGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGCCGAc                                                                    >  1:3034845/1‑49 (MQ=255)
                         tGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGCCGAc                                                                    >  1:4749333/1‑49 (MQ=255)
                         tGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGCCGAc                                                                    >  1:3441124/1‑49 (MQ=255)
                           aCGTGGTAATGGTTGCAGGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAtt                                                  <  1:5468836/65‑1 (MQ=255)
                           aCGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAg                                               <  1:1286878/68‑1 (MQ=255)
                           aCGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                             <  1:4769566/70‑1 (MQ=255)
                           aCGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                             <  1:3727336/70‑1 (MQ=255)
                           aCGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGa                                              <  1:4522049/69‑1 (MQ=255)
                           aCGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGa                                              <  1:5604967/69‑1 (MQ=255)
                            cGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGaaa                                            >  1:7680062/1‑70 (MQ=255)
                              tGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAg                                           <  1:1210409/69‑1 (MQ=255)
                              tGGTAATGGGTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAg                                           <  1:4273032/69‑1 (MQ=255)
                                gTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTgg                                                          <  1:3197353/52‑1 (MQ=255)
                                gTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCg                                         >  1:561153/1‑69 (MQ=255)
                                gTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCg                                         >  1:6358933/1‑69 (MQ=255)
                                 tAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGa                                                               <  1:5948388/46‑1 (MQ=255)
                                 tAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGa                                                               <  1:4613398/46‑1 (MQ=255)
                                  aaTGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGa                                     >  1:3642252/1‑71 (MQ=255)
                                      gTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGAc                                                                    >  1:2157084/1‑36 (MQ=255)
                                      gTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGAc                                                                    >  1:6241846/1‑36 (MQ=255)
                                      gTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATcc                                  <  1:3458724/70‑1 (MQ=255)
                                      gTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                 <  1:289193/71‑1 (MQ=255)
                                         gCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtc                               <  1:4737986/70‑1 (MQ=255)
                                         gCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtc                               <  1:3258556/70‑1 (MQ=255)
                                          cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAtt                                                  >  1:3944661/1‑50 (MQ=255)
                                          cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcgt                             <  1:2140405/71‑1 (MQ=255)
                                          cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcgt                             <  1:4064002/71‑1 (MQ=255)
                                          cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcg                              <  1:1747972/70‑1 (MQ=255)
                                           aCGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCa                                                <  1:4300644/51‑1 (MQ=255)
                                           aCGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCa                                                <  1:3751339/51‑1 (MQ=255)
                                           aCGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCa                                                <  1:521849/51‑1 (MQ=255)
                                           aCGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCa                                                <  1:5180846/51‑1 (MQ=255)
                                            cGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTgc                           <  1:1317212/71‑1 (MQ=255)
                                               ggCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTgcgc                         >  1:7783103/1‑70 (MQ=255)
                                                gCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGt                       <  1:1355659/71‑1 (MQ=255)
                                                gCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGt                       <  1:2278222/71‑1 (MQ=255)
                                                gCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGt                       <  1:3235683/71‑1 (MQ=255)
                                                gCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGt                       <  1:5161919/71‑1 (MQ=255)
                                                gCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGt                       >  1:1932172/1‑71 (MQ=255)
                                                     cTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcgt                             <  1:4945161/60‑1 (MQ=255)
                                                     cTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCt                  <  1:2740600/71‑1 (MQ=255)
                                                      tCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCt                  >  1:1724618/1‑70 (MQ=255)
                                                       cGGTGTTGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTa                 >  1:5100185/1‑70 (MQ=255)
                                                       cGGTGTGGAAATTGGCGACCCGGACCTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTa                 >  1:701207/1‑70 (MQ=255)
                                                       cGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCt                  >  1:6990490/1‑69 (MQ=255)
                                                       cGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTa                 >  1:2122537/1‑70 (MQ=255)
                                                       cGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCCTCAGCTa                 >  1:1818112/1‑70 (MQ=255)
                                                        ggTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcg                              <  1:6931758/56‑1 (MQ=255)
                                                        ggTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTaa                <  1:1334753/70‑1 (MQ=255)
                                                        ggTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTa                 <  1:6107753/69‑1 (MQ=255)
                                                         gtgtGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTaaa               <  1:2465818/70‑1 (MQ=255)
                                                              gAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAg          <  1:2365136/70‑1 (MQ=255)
                                                               aaaTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCg        >  1:1736479/1‑71 (MQ=255)
                                                                aaTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCg                                         <  1:4429248/37‑1 (MQ=255)
                                                                aaTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCgg       >  1:3489268/1‑71 (MQ=255)
                                                                 aTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGAt                                    >  1:3240681/1‑41 (MQ=255)
                                                                 aTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCgg       >  1:6275824/1‑70 (MQ=255)
                                                                 aTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCgg       >  1:5456618/1‑70 (MQ=255)
                                                                 aTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGt      <  1:3276632/71‑1 (MQ=255)
                                                                  ttGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTa     <  1:5799939/71‑1 (MQ=255)
                                                                  ttGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTa     <  1:5197314/71‑1 (MQ=255)
                                                                  ttGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCCCGTCAg                    >  1:5642888/1‑56 (MQ=255)
                                                                     gCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAAt   <  1:3412031/70‑1 (MQ=255)
                                                                     gCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAAt   <  1:2123184/70‑1 (MQ=255)
                                                                     gCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATc  <  1:6375661/71‑1 (MQ=255)
                                                                     |                                                                      
CAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATC  >  AP009048/1940066‑1940205

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: ‑