Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,544,643 T→G 36.2% Q412P (CAG→CCG)  pykA ← Pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,544,6430TG36.2% 26.3 / 3.1 41Q412P (CAG→CCG) pykAPyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base T (3/22);  new base G (7/8);  total (11/30)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCCAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTGGCAGT  >  ACB122/1544509‑1544770
                                                                                                                                       |                                                                                                                               
tCACCCAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGGCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑t                                                                                                                                 <  6:192742/135‑1 (MQ=255)
      aggccccTTGCGCCCCCCGGCCACGTTTGAGCAGGGCGTCGCCCGCCGCCGGGCGCACCTTCTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTgg                                                                                                                           <  5:121238/130‑1 (MQ=255)
      aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGGGCC‑Cgggggg                                                                                                                           >  1:129717/1‑135 (MQ=255)
           ccTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACCCGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGt                                                                                                                     >  7:98877/1‑136 (MQ=255)
           ccTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCGTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGACGCGCGCCCC‑GGGGTGGCGCgg                                                                                                                      >  2:16719/1‑135 (MQ=255)
                 tCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCCGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCCCGGAACATCGAGGCGCGCGCC‑GGGGTGGCGGGGTggggc                                                                                                                >  6:159144/1‑135 (MQ=255)
                 tCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGCACATCGACGCGCGGGCCGGGGGTGGCggggggggg                                                                                                                 >  5:250375/1‑130 (MQ=255)
                      caGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAg                                                                                                           <  2:6091/135‑1 (MQ=255)
                         gCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCg                                                                                                       <  8:244253/136‑1 (MQ=255)
                           cACGTTTGAGCAGGTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                     <  8:116681/136‑1 (MQ=255)
                            acggtggAGCAGTGCGTCGACCGCCGCCTGGCGCACCTTCTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                     <  3:144889/129‑1 (MQ=255)
                            aCGTTTGAGCCGTGCGGCGCCCGCCGCCTGGCTCACGTTAGCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                     <  1:305270/135‑1 (MQ=255)
                            aCGTTTGAGCAGTTCGGCGCCCGCCGCCTGGCTCACCTTAGCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                     <  3:298890/135‑1 (MQ=255)
                                  gAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCg                                                                                                <  3:328901/134‑1 (MQ=255)
                                  gAGCAGTGCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑GGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCg                                                                                                <  5:344859/134‑1 (MQ=255)
                                           gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGCAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGGC‑GGGGGGGCGCGGGGGGGGGCCCGGGCGCAGCTCGGGACATGCgc                                                                                      >  1:99799/1‑134 (MQ=255)
                                           gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑GGGGTGGCGCGGTGCGGCGCCCGGGCGCAGCTCCGGACAAGCGa                                                                                      >  3:199260/1‑135 (MQ=255)
                                                ccgccgCCTGGCTCACCTTCTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑GGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCa                                                                                 <  4:67605/135‑1 (MQ=255)
                                                     gccTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGGCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGgcgc                                                                            <  3:84068/135‑1 (MQ=255)
                                                        gggCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGAc                                                                         <  5:159144/134‑1 (MQ=255)
                                                              aCCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                  <  8:311896/136‑1 (MQ=255)
                                                               ccTTCTCGGCGGGCAGCGCGGGCGGGTCGAAGGCAATCGGGTACACGCCACGGAACCTCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                  <  6:302031/135‑1 (MQ=255)
                                                               ccTTATGGGCCGGCAGCGCGGCCGGGTCGAAGGAAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGAGCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                  <  4:294648/135‑1 (MQ=255)
                                                               ccTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAAGCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                  <  3:47963/135‑1 (MQ=255)
                                                               ccTTAGCGGCCGGCCGCGCGGCCGGGGCGAAGGCAACCGGGTCCACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAg                                                                 <  8:349398/136‑1 (MQ=255)
                                                                       gccggcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGgtgt                                                         <  7:182184/136‑1 (MQ=255)
                                                                           gcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCCCC‑GGGGTGGCGCGGTGCGGCGAGACGGCGAAGAAACGTACAAGCCAACCCAGGCGCGACCAGAATATCGGGGTGTatc                                                      >  6:121635/1‑133 (MQ=255)
                                                                            cAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGcc                                                     <  3:56378/135‑1 (MQ=255)
                                                                            cAGCGCGGCCGGGTCGAAGGCAATCGGGGACACGCCACGGAACATCGCGGCGCGCGCC‑GGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGcc                                                     <  4:256011/135‑1 (MQ=255)
                                                                             aGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGCCAGGGCGAAGCTCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCg                                                    >  3:90266/1‑135 (MQ=255)
                                                                             aGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑GGGGGGGCGCGGTGGGGCGCCCGGGGGAAGCTGGGTACCTGCGCACGCCGGCGCGCCATGATCAGCGGGGGGTAGCCg                                                    >  1:23496/1‑135 (MQ=255)
                                                                                  ggCCGGGTGGACGGCAGTCGGGTACACGCCACGGTACATCGCGGGGCGCGCC‑GGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTc                                               <  5:23417/135‑1 (MQ=255)
                                                                                              ggcacccGGGTCCACGCCCCGGCACATCGAGGCGCGCGCC‑GGGGTGGCGCGGGGCGGCGACAGGGCGTAGATCGGGACATCCGCACGCAGGCGCGACATGATCTGCGGGCTGTAGCCGCTTTCGGTCAGGGCgat                                   <  3:318679/129‑1 (MQ=255)
                                                                                              ggCAATCGGTTACACGCCACGGAACATCGATGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat                                   <  2:152696/135‑1 (MQ=255)
                                                                                                              ccACGGAACATCGAGGCGCGCGCC‑GGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGccc                   <  3:130203/135‑1 (MQ=255)
                                                                                                                     aCATCGCGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt            <  1:16719/135‑1 (MQ=255)
                                                                                                                        tCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt         <  2:235534/135‑1 (MQ=255)
                                                                                                                          gAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg       <  8:43498/135‑1 (MQ=255)
                                                                                                                          gAGGCGCGCGCC‑GGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg       <  8:49963/135‑1 (MQ=255)
                                                                                                                          gAGGCGCGCGCC‑GGGGTGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg       <  1:380879/135‑1 (MQ=255)
                                                                                                                            ggCGCGCGCC‑TGGGGGGCGCGGTGCGGCGACAGGGCGCAGATCGGGTCCTGCGCAACCAGGCGCGACCTGATCAGCGGGGTGTAGCCGCGTTGGGTCATGGCGCTGATCGGCTTTACGCCCGGGAAGGGGTTg       >  4:171227/1‑133 (MQ=255)
                                                                                                                               gcgggcC‑TGGGGGGCGCGGTGCGTCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTGGCAgt  <  1:360380/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                               
TCACCCAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCC‑TGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTGGCAGT  >  ACB122/1544509‑1544770

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑