Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,277,0430TG38.9% 27.3 / 15.7 36G278G (GGT→GGG) hbd3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (6/16);  new base G (14/0);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.32e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGC  >  ACB122/4276910‑4277169
                                                                                                                                     |                                                                                                                              
tGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttc                                                                                                                             <  2:80664/137‑1 (MQ=255)
               cccTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTATTAACCAACTgagga                                                                                                              >  6:221933/1‑137 (MQ=255)
                   gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCGACGGAGGAGCaa                                                                                                          >  6:324305/1‑137 (MQ=255)
                   gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCAAATGAGGAGCaa                                                                                                          >  8:48494/1‑137 (MQ=255)
                   gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGGTATTCCCCTACGTCGGAGCaa                                                                                                          >  4:4078/1‑137 (MQ=255)
                      ggCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAgcg                                                                                                       <  7:256573/137‑1 (MQ=255)
                      ggCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTATTCGACGACTGAGGTGGAAgcc                                                                                                       >  2:309708/1‑136 (MQ=255)
                      agaccgcATGGTGGCCTTCCACGAAGGTTTAAACGACCGCAAATACCGCCCTGCACCACTGCTCAAGGAAAGGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAgcg                                                                                                       <  1:286868/131‑1 (MQ=255)
                             aTGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTTTTGAAATACTGGTGGTCAATCGCCAAGcc                                                                                                >  4:324137/1‑137 (MQ=255)
                             aTGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCTACTGAGGAGCAAGCGCCATGcc                                                                                                >  4:168057/1‑137 (MQ=255)
                             aTGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCGACTGGGGGGCAAGCGcccagcc                                                                                                >  6:187527/1‑132 (MQ=255)
                                                gCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATg                                                                             <  3:151940/137‑1 (MQ=255)
                                                  ttCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGTTCACCTACTGAGGCGCAAGCGCCATGCCCCCGGTCAATGCCGCAATgcg                                                                           >  8:212979/1‑137 (MQ=255)
                                                     aaCGACAGCATATACCGCCCTGCTCAACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtg                                                                        <  4:229854/137‑1 (MQ=255)
                                                             cAAATACCGCCATGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACt                                                                <  3:168057/137‑1 (MQ=255)
                                                              aaaTACCGCCCTGCTCCACTGCTCAAGGAAATGGTCTCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGAACAACTCAGGCGCAACCCCCATGCCCCCGGTAAATGCCTCCAAGCGCTGGGccccagt                                                               >  6:5196/1‑132 (MQ=255)
                                                                        cTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTCGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCTAATGACGAGCAACCGCCAAGCCCCCCGTCAAAGCCGCCATGCGCGGGGCCAACTTTGAGGAGCgc                                                     >  5:269167/1‑136 (MQ=255)
                                                                         tGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGAc                                                    >  6:95743/1‑137 (MQ=255)
                                                                         tGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCCCCGAATGAGGAGCAACCCCCAAGCCCCCGGTGAAAGCCGCCAAGCCCGGTGCCAACGATGACGAGcccc                                                    >  8:289053/1‑134 (MQ=255)
                                                                         tGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCACCTACGGAGGAGCAAGCGCCATGCCACCGGTCAATGCCGCAATGCGCTGTGGCAACTTTTAGGAGCGAc                                                    >  6:356942/1‑137 (MQ=255)
                                                                           cTCCACTGCTCACGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                  <  3:299772/137‑1 (MQ=255)
                                                                           cTCCACTGCTCAAGTAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                  <  7:285048/137‑1 (MQ=255)
                                                                           cTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                  <  8:181212/137‑1 (MQ=255)
                                                                           cTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                  <  8:64299/137‑1 (MQ=255)
                                                                                        ggAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGAACAAAAGAGGAGCAAACCCCCATCCCCCCGTCAAATCCGCCAATCGCGGTGGCCACGTTTAGGAGCGCCCTGACCGGGCCATgg                                     >  6:92113/1‑137 (MQ=255)
                                                                                        ggAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTCACCTACTGAGGAGCAATCGCCCAGCCCCCCGTCAAAGCCGCAATTCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATgg                                     >  4:261468/1‑137 (MQ=255)
                                                                                          aaaTGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTTACCTACGGAGGTGCAAGCGCCCTGCCCCCGGTCACTGCCGCAATGCGCTGTGCCCACTTTGAGGCGCGACCTGCCCGGGCCCTGGcc                                   >  7:146169/1‑136 (MQ=255)
                                                                                           aaTGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCAt                                  >  7:277798/1‑137 (MQ=255)
                                                                                                tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                           >  4:271185/1‑43 (MQ=255)
                                                                                                tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                           <  3:271185/43‑1 (MQ=255)
                                                                                                     gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                        <  5:92113/137‑1 (MQ=255)
                                                                                                     gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                        <  5:221933/137‑1 (MQ=255)
                                                                                                     gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                        <  4:11553/137‑1 (MQ=255)
                                                                                                               tggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCAAAAAGGGCttt              <  7:28486/137‑1 (MQ=255)
                                                                                                                  ggCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGc           <  3:84556/137‑1 (MQ=255)
                                                                                                                           gCGGTCGCGGGTTGTTCACCAACAGAGGAGCAAGCGCCAAGCCCCCGGTCAAAGCCCCCAAGCGCTGTGCCAACTTTGAGGCGCGACGTGGCCGGGCCAAGGCATTTTTTTCCGAAAAGGGCTTTTGCCAGGTAAGc  >  6:181364/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                              
TGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGC  >  ACB122/4276910‑4277169

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