Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,544,643 T→G 33.0% 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,6430TG33.0% 40.3 / 6.6 51Q412P (CAG→CCG) pykAPyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base T (6/26);  new base G (10/7);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 9.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTGGCAGTGTAC  >  ACB122/1544515‑1544774
                                                                                                                                |                                                                                                                                    
tctCACCTTCCTCCACCAGGCCACGTTTGAGCAGTGCGTCGACCGCCGCCTGGCGCACCTTATCCGCCTGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACCTCGAGGCGCGCGCCTgggggg                                                                                                                                <  4:126198/133‑1 (MQ=255)
aggcACCTTCCTCCCCCAGGCCACGTTTGAGCAGTGCGGCGACCGCCGCCTGGCGCACCTTATCGGCCGGCCGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                                <  7:9520/132‑1 (MQ=255)
aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTGCACGCCCCGGAACATCGAGGCGCGCGCCGGGGTgg                                                                                                                                >  2:335467/1‑135 (MQ=255)
aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGCGGCGCGCGCCGgggggg                                                                                                                                >  1:332799/1‑135 (MQ=255)
aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                                <  5:403493/135‑1 (MQ=255)
     cattGCTCCACCAGGCCACGTTTGAGAAGTTCGTAGACCGCCGCCTGGCGCACCTTATCGGCCGGCCGCGCGGCCGGGTCGAAGGCAGACGGGTACACGCGACGGACCGTCGAGGCGCGCGCGTGGGGGGCGCgg                                                                                                                           >  2:296212/3‑135 (MQ=255)
       ttGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTg                                                                                                                         <  3:107464/135‑1 (MQ=255)
           tCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTgcggc                                                                                                                     >  7:268139/1‑135 (MQ=255)
              accaGGCCAGGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTgcggcg                                                                                                                    <  3:136487/133‑1 (MQ=255)
                      aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGTCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                          <  7:59035/135‑1 (MQ=255)
                      aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                          <  2:247023/135‑1 (MQ=255)
                      aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                          <  4:332477/135‑1 (MQ=255)
                      aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                          <  2:419449/135‑1 (MQ=255)
                      aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGCCAGGGc                                                                                                             >  2:47800/1‑132 (MQ=255)
                      aCGTTTGAGCAGTGCGTCGACCGCCGCCTGGCTCAGCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                          <  6:213852/135‑1 (MQ=255)
                          ttGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGGCCCCGGGGGGGGGGGGGCGGCGCCACCGCGGAGATc                                                                                                      >  6:35663/1‑135 (MQ=255)
                                     gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                                >  7:66948/1‑98 (MQ=255)
                                     gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                                <  8:66948/98‑1 (MQ=255)
                                          ccgccgCCTGGCGCACCTTATCGGCCGGCAGCGCGGCCGGGGCGAAGTCTATCGGGTACACGCCACTGAACCTCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGAGCGGGACTTGCGAACGCa                                                                                      <  8:453807/135‑1 (MQ=255)
                                                         cgttATCGGCCGGGAGCGCGGCCGGGTCGATGGCAATCGGGTACACGCCACGGATCCTCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGAgatgat                                                                         <  4:115131/131‑1 (MQ=255)
                                                         ccTTCTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGGGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGAt                                                                         <  7:95127/133‑1 (MQ=255)
                                                         ccTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                       <  4:55511/135‑1 (MQ=255)
                                                         ccTTACCGGCCGGCCGCGCGGCGGGGGGGACGGCAAGGGGGTCCACTCCACGGAACAGCGCGGCGCGCGCCGGGGGGGCGCGGGGCGGCGACAGGGCGAAGAGCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                       <  1:390272/135‑1 (MQ=255)
                                                             aTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGGACATGCGCACGCAGGCCCGCCATGATCAGCgg                                                                   >  7:223120/1‑135 (MQ=255)
                                                                ggccggccGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGgt                                                                <  3:129240/135‑1 (MQ=255)
                                                                       aGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGGGGGGGGGCGCCAGGGCGAAGCTAGGTACATGCGCACGCCGGGGCGCCATGGTCAGCGGGGTGTAGGCg                                                         >  8:85358/1‑135 (MQ=255)
                                                                           cgGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGGGCGGCGCCCGGGGGAAGATCGGGACATGCGAACGCCGGCGCGCCATGATCAGCGGGGTGTTGCCGCtgt                                                     >  4:129240/1‑133 (MQ=255)
                                                                            ggCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTc                                                    <  4:128733/135‑1 (MQ=255)
                                                                              ccGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGCGGAGAGCGGTTCATGCGCACGCAGGCGCGCCATGATCCGCGGGGTGTAGCCGCGTTCgg                                                  >  6:172932/1‑135 (MQ=255)
                                                                                        ggCAATCGGGTACACGCCACGTAACATCGAGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat                                        <  1:300264/135‑1 (MQ=255)
                                                                                           aaTCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATTCCAACCCAGGCGCGACATGATCACCCGGGTGTAGCCGCTTTCGGTCAGGGCgatgat                                     >  8:378359/1‑135 (MQ=255)
                                                                                                ggtccccGCCACGGAACTTCGAGGCGCGCGCCTGGGGTGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑t                               <  6:188301/129‑1 (MQ=255)
                                                                                                          aCGGAACCGCGAGGCGCGCGCCGGGGGGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGGGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCgg                     <  1:192380/135‑1 (MQ=255)
                                                                                                            ggAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGGCAGGGGGAAGATCGGTACATGCGAACGCAGGCGCGCCATGGTGATCGGGGTGTTGCCGCGTTCGGTCAGGGCGGAGATAGCCGTTTCCGccgg                      >  4:179984/1‑131 (MQ=255)
                                                                                                             gAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCCGGCGCGACATGAAAAGCGGGGTGTAGCCGCGTTCGGTCAGGGCGATGATCGCC‑ATTACGCCCgggga                  >  2:360333/1‑133 (MQ=255)
                                                                                                              aaCATCGAGGCGCGCGCCTGGGTGGCGCGGTGGGGCGACAGGGCGAAGATCGGTTCATGCGAACGCAGGCGCGACATGATCCGCGGGGTGTAGACGCTTTCGGTCAGGGCGCTGATCGCC‑GTTACGGCCGGGAAg                 >  3:343612/1‑135 (MQ=255)
                                                                                                              aaCATCGAGGCGCGCGCCAGGGTGGCGCGGTGCGGCGACAGGGCGCAGATCGGTACATGCGAACGCCGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTAGGTCAGGGCGATGATAGCC‑TTTACGCCCGGGAAg                 >  3:230474/1‑135 (MQ=255)
                                                                                                               acacCGCGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGt                <  1:390577/135‑1 (MQ=255)
                                                                                                                  tCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGt             <  4:241035/135‑1 (MQ=255)
                                                                                                                  tCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGt             <  6:86083/135‑1 (MQ=255)
                                                                                                                  tCGAGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGt             <  3:84084/135‑1 (MQ=255)
                                                                                                                    gtggCGCGCGCGGGGGTGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTg           <  6:73550/133‑1 (MQ=255)
                                                                                                                    gAGGCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGTGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTg           <  8:182883/135‑1 (MQ=255)
                                                                                                                     cggCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGGGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTg           <  8:213682/133‑1 (MQ=255)
                                                                                                                     aGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTg           <  4:351070/134‑1 (MQ=255)
                                                                                                                     aGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTg           <  8:268136/134‑1 (MQ=255)
                                                                                                                     aGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTg           <  8:268139/134‑1 (MQ=255)
                                                                                                                      ggCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                       >  8:98789/1‑74 (MQ=255)
                                                                                                                      ggCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                       <  7:98789/74‑1 (MQ=255)
                                                                                                                             gcCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCCGGGCGATTATCGCC‑TTTACGCCCGGGAAGTGGTTGGCAGTGTAc  >  8:252299/1‑135 (MQ=255)
                                                                                                                             gcCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTGGCAGTGTAc  <  7:395554/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                    
AGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCC‑TTTACGCCCGGGAAGTGGTTGGCAGTGTAC  >  ACB122/1544515‑1544774

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