Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,544,6430TG55.2% 11.3 / 26.1 31Q412P (CAG→CCG) pykAPyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base T (2/11);  new base G (9/7);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 5.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-03
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTG  >  KT2440/1544515‑1544765
                                                                                                                                |                                                                                                                          
aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCTGCCGGCAGCGAGGCCGGGTCGAAGGCAATAGGGTACACGCCCCGGCACATCGAGGCGCGCGCCGGGGTgg                                                                                                                      >  6:253631/1‑135 (MQ=255)
aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGCAGGCAATCGGGTACCCGCCACGGAACCTCGAGGCGCCCGCCCgggggg                                                                                                                      >  7:108083/1‑135 (MQ=255)
     ccTTGCTCCACCAGGCCAAGTTTGAGCAGTTCGTCGACCCCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAACGAAATCCGGTACACGCCCCGGAACATCCAGGCGCGCGCCAGGGGGgggccg                                                                                                                 >  8:234208/1‑133 (MQ=255)
           tCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTgcgg                                                                                                            <  2:10203/134‑1 (MQ=255)
                 aGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAgg                                                                                                     <  3:273634/135‑1 (MQ=255)
                   gCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGc                                                                                                   >  7:228786/1‑135 (MQ=255)
                                aGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGCTCGGTACa                                                                                      >  7:250072/1‑135 (MQ=255)
                                     gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGGGGCGCCAGGGCGAAGCTCGGTACAAGCgc                                                                                 >  8:20776/1‑134 (MQ=255)
                                     gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGCGGTGCGGCGGCAGGGCGCAGATCGGTTCATGCGa                                                                                 >  3:103151/1‑135 (MQ=255)
                                     gtcgACCGCCGCCTGGCTCAACTTATCGGCCGGCCGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCCCCGAACATCGAGGCGCGCGCCGTGGTTGCGCGGGGGCGCGCCACGGCGAACATCGGGGCCTGCGa                                                                                 >  8:177619/1‑135 (MQ=255)
                                               gccTGGCTCACGTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGAACATGCGAACGCAGgcgc                                                                       <  6:71294/135‑1 (MQ=255)
                                                         ccTTTTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATGGGCTACACGCCACGGAACAGCGAGGCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCGGGCGCGACATGATCa                                                             <  6:12079/135‑1 (MQ=255)
                                                                 gccggcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCCAAGAACCGTTCATGGGAACGCAGGCGCGAAATTATACACGGGgt                                                      >  8:229642/1‑134 (MQ=255)
                                                                 gccggcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGCGGTGCGGCGCCACGGGGAAGAACCGTTCATGCGAACGCAGGGGCCGCAAGATAACCGGGgtg                                                     >  8:175922/1‑135 (MQ=255)
                                                                       aGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCCCGGAACATCGAGGCGCGCGCCGGGGGGGGGCGGGGCGGCGCCCGGGCGCAGATCGGTAACTGCGCCCGCAGGCGCGACAAGATCAGCGGGGTGTAGCCg                                               >  2:101974/1‑135 (MQ=255)
                                                                            ggCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGCGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTc                                          <  2:273429/135‑1 (MQ=255)
                                                                             gCCGGGTCGTAGTCAATGGGGTAGAGGCGACGGAAGATCTAGGCCCGCGCGTGTGTGGCGCGGTGCGGGGACTGGGCGAAGATCGGTTCATGCGAACGCAGGCTCGACATGATCATCGGGGTGTAGCCGCTTTCg                                         <  2:276904/135‑1 (MQ=255)
                                                                                        ggcatgcGGGAACACGCCACGCACCAGCGCGGCGCGCGCCGGGGGGGGGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat                              <  1:253020/129‑1 (MQ=255)
                                                                                        ggCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat                              <  7:26580/135‑1 (MQ=255)
                                                                                        ggCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGAGCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat                              <  8:122802/135‑1 (MQ=255)
                                                                                        ggCAATCGGGAACACGCCAGGGAACTTCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGATGTTCGGTACATGCGAACGCAGGCGCGTCATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGTCgat                              <  7:177619/135‑1 (MQ=255)
                                                                                         gCAAGCGGGAACCCGCCACGGAACATCGCGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgatg                             <  5:243052/135‑1 (MQ=255)
                                                                                               gggTACACGCCACGGAACATCGAGGCGCGCGCCGGGGGGGCGCGGGGGGGCGCCAGGGCGAAGATCGGTTCATGCGAACCCAGGCGCGACATGATAAGCGGGGTGTAGCCGCGTTCGGGGACGGCGATGATCGcc                       >  2:47801/1‑135 (MQ=255)
                                                                                                ggTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCt                      <  8:228786/135‑1 (MQ=255)
                                                                                                       gCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGcc               >  7:194638/1‑135 (MQ=255)
                                                                                                               ccagcGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt       <  1:116470/131‑1 (MQ=255)
                                                                                                               acctCGCGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt       <  5:186683/132‑1 (MQ=255)
                                                                                                               aCATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt       <  7:20776/135‑1 (MQ=255)
                                                                                                                  tCGAGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt    <  2:285880/135‑1 (MQ=255)
                                                                                                                    gcggggcgcgcGGGGGGGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt    <  2:297150/128‑1 (MQ=255)
                                                                                                                     aGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCTAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg  <  8:44788/134‑1 (MQ=255)
                                                                                                                                |                                                                                                                          
AGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTG  >  KT2440/1544515‑1544765

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: