Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,279,0740TG48.1% 24.4 / 10.6 27G278G (GGT→GGG) hbd3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/14);  new base G (13/0);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.99e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCA  >  KT2440/4278942‑4279193
                                                                                                                                    |                                                                                                                       
gagcGGCCTGGACACCCTGCTGGCAATCATGGCGGCCTCCCACGACGGCTTAAACGACAGAAAATACCCCCCTGCACCACTGCTCACGCAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGttt                                                                                                                       <  7:19097/132‑1 (MQ=255)
   cGGCCTGGACACCCTGCTGGAATTCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGTAAATGGCCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                    <  7:245924/135‑1 (MQ=255)
    ggCCTGGACAGCATGCTGGCATTCTTGCAGGCCTGCGACGAAGGCTTAAACGACACCTATTACCGGCCTGGCCCAGTGCTGAAGTTAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                    <  8:285684/134‑1 (MQ=255)
    ggCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGgttgttg                                                                                                                   >  3:285651/1‑134 (MQ=255)
    ggCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCGGGGGCGCAAGAGCGGTCGCGGgttgttg                                                                                                                   >  7:128005/1‑134 (MQ=255)
         ggACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTccccg                                                                                                               >  3:48656/1‑133 (MQ=255)
         ggACACACTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGgttgttgagc                                                                                                                >  5:319470/1‑129 (MQ=255)
         gcacacCCGGATGGAAAACATCGCGCCCGCCCACGAAGGCTTAAACCACACCAAAAACCGCCCTGCACCACTGCTCACGTAATTGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCTGTTTCTTCACCTa                                                                                                              <  8:323699/133‑1 (MQ=255)
                  gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGGTGACCTAACGAGGGGc                                                                                                     >  3:313169/1‑135 (MQ=255)
                  gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCGTCACCGACAgaggcg                                                                                                      >  2:203272/1‑132 (MQ=255)
                  gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTATTCACCGACTgaggag                                                                                                      >  6:140936/1‑134 (MQ=255)
                     ggCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCTACAGAGGAGCAAg                                                                                                  >  4:136402/1‑135 (MQ=255)
                       cAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGGTGGTGACCAACTgaggagacaccc                                                                                                >  8:364702/1‑129 (MQ=255)
                              ggAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTATTGACCACATGTGGGGCAAGCCCCaagac                                                                                         >  7:346334/1‑131 (MQ=255)
                                        cACGAAGGCTCCAACGACAGCAAATACCCCCCTGCTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCa                                                                                 <  7:118860/133‑1 (MQ=255)
                                         aCGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCTAAGGAGGAGCAAACGCCCATCCCCCCGTCAATg                                                                              >  8:42904/1‑135 (MQ=255)
                                                 ttCAACGACACCAAATACCGCCCGCGTCCACTCCTCACGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATg                                                                      <  4:140821/135‑1 (MQ=255)
                                                 tgcAACGACAGCAAACACCGCCCGGCTCCACTGCTCACGTAATTGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATg                                                                      <  1:38853/133‑1 (MQ=255)
                                                      cGACACCAAATCCCGCCCGCCTCCACTGCCCACGGAAGGGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtg                                                                 <  1:113935/135‑1 (MQ=255)
                                                            cAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCGACTGGGGGCGCAGCGCCCTGCCCCCCGGCAATGGCGCCAAGCGCGGGGccc                                                            >  7:225981/1‑133 (MQ=255)
                                                                   cGCCCTGCGCCACTGCTCAAGGAAAGGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTgag                                                    <  8:284460/135‑1 (MQ=255)
                                                                    gCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGCCCGACGGAGGCGCCAGCGCCCTGCCCCCCGTCCATGCCGCCAAGCGCGGTGCCAACGTTgtgg                                                   >  7:224775/1‑135 (MQ=255)
                                                                    gCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCACCGACTGAGGAGCAAGCGCCAAGCCCCCCGTCCATGCCGCCATGCGCTGTGCCAACGTTgagg                                                   >  2:213185/1‑135 (MQ=255)
                                                                       cTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCTACGGTGGGGCAAACGCCCAGCCCCCGGGGAATGCCGCCATGCGCGGGGCCCACTTTGAGgcgc                                                >  3:132881/1‑135 (MQ=255)
                                                                            gaccggCTCAAGCATAGGGACGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                           <  7:45329/130‑1 (MQ=255)
                                                                            ccAGTGCTCAAGGAATGGGCCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                           <  1:326221/135‑1 (MQ=255)
                                                                                               tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCGACTGAGGGGCAAGCGCCCAGCCCCCGGGGAATTCCGCCAAGCGCTGTGGCCACTTTGAGGGGCGCCGGGGCCGGGGCCTgggaattt                        >  6:143948/1‑129 (MQ=255)
                                                                                                      cggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                 <  1:105490/135‑1 (MQ=255)
                                                                                                             ctggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAgg            <  7:364702/133‑1 (MQ=255)
                                                                                                             ctggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAgg            <  1:159371/133‑1 (MQ=255)
                                                                                                                     cAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCa  <  4:285651/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                       
GATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCA  >  KT2440/4278942‑4279193

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