Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,626,1120TG52.6% 11.4 / 19.4 21M176R (ATG→AGG) hcf136_3Ycf48‑like protein
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base G (0/10);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGC  >  KT2440/3625994‑3626236
                                                                                                                        |                                                                                                                            
gAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAACCCGTTGTTGGCCCTGAAGTTCACCGATGCACGCAACGGCCTGGTCGTAGGTGCCGGCCGGGTG‑G‑CCAAGCCCCCCGACGAccgc                                                                                                              >  8:315045/1‑135 (MQ=255)
     gcGCCGGGGGGGGGCCAAACGGCGGCGGGCCGACGGCGCGGACAAGCCGTTGTGGGCCCTTAATTACCCCGAGCCCCGCACCGGCCTGGGCGGCGGGGCCTCCGGGGGG‑G‑CCAGGCCCACCGACGACggcgggcg                                                                                                         <  2:302291/134‑1 (MQ=255)
                                ggCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTG‑G‑CCATGCGCACCGCCGCCGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACt                                                                              >  1:6075/1‑135 (MQ=255)
                                ggCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACt                                                                              >  6:150937/1‑135 (MQ=255)
                                   ggcgggCGCGTACACGCCGTTGTTGCCCCTTATTTCCACCGAGGCCCGCACCGGCCTGGGCGTGGGGGCCTAGGGGGGG‑G‑CCTGGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                           <  8:29086/131‑1 (MQ=255)
                                   cgccGGCGCGCACACGCGGTTGTGGCCCCTCAATTCCCCCGCGCCCCCCACCGGCGGGGGGGGGGGT‑CCCACGGGGGG‑GCCCAGGCCCCCCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                           <  8:28493/132‑1 (MQ=255)
                                   cGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGAGGCCCGCAACGGCCTGGTCGGCGGTGCCTACGGGGTG‑G‑CAATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                           <  8:283908/135‑1 (MQ=255)
                                         ggcgcACACGCGGTTGTGGGCCCTCAATTCCCCCTCGCCCCCCACCGGCGGGGGCGGCGGGGCCTCCGGGGGG‑C‑CCAGGCCCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACgg                                                                     <  6:93336/134‑1 (MQ=255)
                                         cgcgGACAAGCCGTTGTTGGCCATGAATTTCCCCGATCCCCGCACGGGCGTGGGCGTGGGTTCATACGGGGTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACgg                                                                     <  8:101409/135‑1 (MQ=255)
                                         cgcgCACACGGCGTGGTTGCCCCTGAATTCCCCCGAGCCCCCCACCGCCGGGGGCGGCGGGGCCTCGGGGGGG‑C‑CCAGGCGCCCCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACgg                                                                     <  6:98484/135‑1 (MQ=255)
                                             gACAAGCGGTTGGGGGCCCTGAATGTCCCCCAGGCCGGCACGGGCGGGGCGGGCGGGGCCTCCGGGGGG‑C‑CCAGGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGt                                                                 <  7:128706/135‑1 (MQ=255)
                                                acgCCGTTGGTGGCCATCAATTTCCCGGTGGCCTGCAAGGGCGTGTTCGTCGGTCCCTCCGGGGGG‑G‑CCTTGGGCACGGACGACGGCGGGCGCAGCTGTCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGa                                                              <  8:351680/133‑1 (MQ=255)
                                                acgCCGTGTTTGCCCCTGAATTGCCCCCCGCCCCCCACCGGCGTGGGGGTGGGTCCCGCGGGGGGG‑G‑CCAGGCCCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGa                                                              <  3:84033/133‑1 (MQ=255)
                                                                 gAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACCGGGTG‑G‑CCAAGCGCCCCCACGACCGCGGGCGCCGCAGGCAGTCGGAGATCGGCCCACTGCCGAACCGGCCTGGCCAGCACCAGTATGccca                                             >  2:332968/1‑134 (MQ=255)
                                                                              ccccGCACCGCCGGGGGCGTCGGGGCCTA‑GGGGTGCG‑CCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                                <  8:200490/134‑1 (MQ=255)
                                                                               cccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCa                               >  7:277747/1‑135 (MQ=255)
                                                                                                    gTGCCTCCGGGGGG‑G‑CCATGCCCACCGACGACGGCGGGCGCAGCTGGCAGGCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCg          <  3:223545/135‑1 (MQ=255)
                                                                                                      gCCTACGGGGTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGCCTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCCAGGCCACCCGGTGTACATCGCCGGCgcg        >  1:164727/1‑133 (MQ=255)
                                                                                                       ccTCCGGGGGG‑G‑CATGGCCCACCGCCCACGGCGTGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGc       <  7:202409/135‑1 (MQ=255)
                                                                                                            ggggTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc  >  1:130247/1‑135 (MQ=255)
                                                                                                            ggggTG‑C‑CCTGGCGCACCGAGGACGGCGGGCGCAGCTGGCAGTCGGAGAGCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc  <  6:127698/135‑1 (MQ=255)
                                                                                                                        |                                                                                                                            
GAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTG‑G‑CCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGC  >  KT2440/3625994‑3626236

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