Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,626,0790AC52.2% 9.8 / 22.1 23N165T (AAC→ACC) hcf136_3Ycf48‑like protein
Reads supporting (aligned to +/- strand):  ref base A (8/3);  new base C (2/10);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAAGCGCGCCGCCGCCCTGGAGGTGCAGGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGC  >  ACB122/3625945‑3626206
                                                                                                                                      |                                                                                                                               
cAAGCGCGCCGCCGCCCTGGAGGTGCAGGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACTGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCaa                                                                                                                                 >  6:72548/1‑135 (MQ=255)
   gcgcgcCGCCGCCCTGGAGGTGCAGGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGCCCAGCCGTTGTTGGCCCGGAATTTCACCGATGCCCGCCACgg                                                                                                                              >  5:359854/1‑135 (MQ=255)
    cgcgcCGCGGCCGTGTTGGGGCAGGCCGCGGTGGCGGCGTCAGACGCGCGGCGCGGGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGc                                                                                                                             <  3:202668/135‑1 (MQ=255)
     gcgcCGCCGCCCTGGAGGTGCAGGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCCCCGGcc                                                                                                                            >  5:122850/1‑135 (MQ=255)
         cgccgccCTGGAGGGGCTGGCAGCCGAGGCGGCGGAAGACGAGCGTCGCCGGGGGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGgt                                                                                                                        <  4:409907/135‑1 (MQ=255)
                                   aGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCCTgc                                                                                              >  4:385889/1‑135 (MQ=255)
                                                                    cgcAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCCCCGATGCCCGCACCGGCCTGGTCGTCGGTGCCTCCGGGGTGGCCATGCGCACCGCCGACGGCGGGCGCCGCTGGCAGTc                                                               >  2:333410/1‑133 (MQ=255)
                                                                     caaccggcggCTGGCCGCCGGCGCGCACACGCCGTTGTTGGCCCTTAATTTCCCCGAGGCCCGCACCGGCCGGGGGGGCGGGGCCTCCGGGGGGCCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGga                                                            <  4:304666/131‑1 (MQ=255)
                                                                                 ggCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACt                                                >  6:211855/1‑135 (MQ=255)
                                                                                 ggCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCGACGGGGTGGCCATGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGCCt                                                >  3:430278/1‑135 (MQ=255)
                                                                                   ccgcgggCGCGGACAAGCCGTGGTGGGCCCTCAATTCCCCCCCGCCCCGCACCGGCCGGGGCGGCGGGGCCTCGGGGGGGGCCAGGCCCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                             <  7:332804/131‑1 (MQ=255)
                                                                                   ccGACGGCGCGCACACGCCGTTGTGGCCCCTCAATTCCCCCGAGGCCCGCACCGGCCGGGGCGGCGGGGCCTCGGGGGGGGCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                             <  8:269886/136‑1 (MQ=255)
                                                                                    cgccggcggGGCCAAGCGGTTGTTGGCCCTCATTTACCCCGAGCCCCCCACCGGCGGGGCCGGGGGGGCCTCGGGGGTGGCCAGCCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                             <  4:265869/132‑1 (MQ=255)
                                                                                    cGACGGCGCGGACACGCCGTTGTTGGCCCTGAATTCCCCCGAGGCCCGCACCGGCGGGGCCGGCGGGGCCTCCGGGGTGCCCAGCCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                             <  3:328290/135‑1 (MQ=255)
                                                                                        ggCGCGCACACGCCGTGGTGGCCCCTCAATTTCCCCGAGCCCCGCACCGGCCGGGGCGGCGGTGCCTCGGGGGGGGCCATGCCCCCCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACg                                        <  8:24131/136‑1 (MQ=255)
                                                                                                 aaGCGGTGGTTGGCCCTTATTTGCCCCGGTCCCCGCCCCGGCGGGGTGGTGGGTGCCTCCGGGGGGGCCATGCGCCCGGAGGACGGCGGGCGCAGCTGGCAGGCGGAGATCGGCCGACTGCCGAACGGGCGTGGa                                <  4:335728/135‑1 (MQ=255)
                                                                                                 aaGCCGTTGTTGGCCCTGAATTTCACCGCTGCCCGCAACGGCCGGGGCGTCGGTGCCTCCGGGGGGGCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGa                                <  5:209906/135‑1 (MQ=255)
                                                                                                    cggttgttGGCCCTCAATTTCACCGATGCCCGCACCGGCCTGGGCGGCGGGGCCTACGGGGGGGCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTg                             <  4:369733/133‑1 (MQ=255)
                                                                                                      gttgttGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGGCTACGGGGTGGCCATGCGCACCGACGACCGCGGGCGCACCAGGCAGTCGGAGAACGGCCGACTGCCCAACGGGCGTGGACTGCAc                          >  7:18909/1‑136 (MQ=255)
                                                                                                      gtggttgGCCCTTAATTCCCCCGTTGCCCCCACGGGCCGGGGGGTGGGGGCCCCCGGGGGGGCCTGGCGCCCCGCGGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCa                           <  3:176721/132‑1 (MQ=255)
                                                                                                            ggCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTAt                    >  7:266378/1‑136 (MQ=255)
                                                                                                                     tttCACCGATGCCCGCAACGGCCTGGTCCTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACCGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACAGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGgc            >  5:424263/1‑135 (MQ=255)
                                                                                                                           cGATGCCCGCAACGGCCTGGTCGTCGGTGGCGACCGGGTGGCCATGCGCCCCCCCGCCGGGGGGCGCCGCTGGCAGTCGGCGATCGGCCGCCTGCCGCACCGGCGGGGGCTGCACCAGTATGCCCGGGCGCGccc      >  4:249350/1‑134 (MQ=255)
                                                                                                                             aggcccgcccCGGCGTGGGCGGGGGGGCCTACGGGGGGGCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAg    <  2:293269/133‑1 (MQ=255)
                                                                                                                              ggCCCGCACCGCCGGGGTCGGGGGTCGCTCCGGGGGGGCCTTGCGACCCTACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGGGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc  <  8:405398/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                               
CAAGCGCGCCGCCGCCCTGGAGGTGCAGGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGC  >  ACB122/3625945‑3626206

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