Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,626,117 A→C 50.3% T178P (ACC→CCC)  hcf136_3 → Ycf48‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,626,1170AC50.3% ‑1.2 / 11.9 25T178P (ACC→CCC) hcf136_3Ycf48‑like protein
Reads supporting (aligned to +/- strand):  ref base A (3/9);  new base C (5/8);  total (8/17)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCT  >  KT2440/3625989‑3626244
                                                                                                                                 |                                                                                                                               
cAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCCCCGATGCCCGCAACGGCCTGGTCCTCGGTGCCGACGGGGTGGCCATGCGCCCccacg                                                                                                                           >  5:175763/1‑135 (MQ=255)
 aGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGGCTACGGGGTGGCCATGCGCCccgccgc                                                                                                                          >  2:155119/1‑134 (MQ=255)
 aGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCGTCGGGGTGGCCCAGCGCACcgacga                                                                                                                          >  2:317559/1‑135 (MQ=255)
        cGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTACTAGGTGCCTACGGGGTGGCCATTCGCACCGACGACggcggg                                                                                                                   >  4:439609/1‑135 (MQ=255)
             cggggcggggccCCAAGGGCGGCGGGCCGACGGCGCGGACACGCCGTTGTGGCCCCTTAATTTCCCCGATCCCCCCACCGGCCTGGGCGGCGGGGCCTCCGGGGGGGCCATGCGCCCCGACGACGGCGGGCGCAGc                                                                                                             <  8:441046/133‑1 (MQ=255)
                                       ccGAGGGCGCGCACAAGCGGTTGTGGGCCCTCAATTGCCCCGAGGCCCGAACCGGCGGGGGGGGGGGGGCCTCGGGGGGGCCAATGCGCACCTACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                                   <  8:249667/136‑1 (MQ=255)
                                        ggccGGCGCGGACACGCGGTTGTGGGCCCTGAATTGCCCCGAGGCCCGCACCGGCCTGGGCGGCGGGGCCTCGGGGGGGGCCAGGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                                   <  2:44635/132‑1 (MQ=255)
                                        cgccGGCGCGCACACGCGGTTGTGGGCCCTTAATTCCCCCGAGCCCCGCACCGGCGGGGCGGGGGGTGCCTCGGGGGGGGCCAGGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                                   <  4:395579/132‑1 (MQ=255)
                                        cGAGGGCGCGCACAAGCCGTTGTTGCCCCTCAATTTCCCCCAGCCCCCCACCGGCCTGGGCGGCGGGGCCTCCGGGGGGCCCAGGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                                   <  3:184274/135‑1 (MQ=255)
                                              ggcgGACACGGGGTTGTTGCCCCTCAATTCCCCCGATCCCCGCACCGGCGGGGTGGGGGGTGCCTCGGGGGGGGCCAGCCGCCCGGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACgg                                                                             <  1:72596/134‑1 (MQ=255)
                                               gcgGACACGCGGTTGTGGGCCCTGAATTCCACCGAGGCCCCCAACGGCCGGGGGGGGGGTGCCTCCGGGGGGGCCTGGCGCCCCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGc                                                                           <  8:351514/136‑1 (MQ=255)
                                                   aCAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCGACGGGGTGGCCATGCGCCCCGCCGCCGGCGGGCGCAGCTGGCAGTCGGAGCTCGGCCGCCTGCCGCCCGGGCGTg                                                                        >  5:146012/1‑135 (MQ=255)
                                                    cAAGCGTTTGTTGGCCCTTAATTACCCCCAGGCCCGCACGGGGGGGGGCGTGGGTGCTTCGGGGGTGGCCATGCCCACCGTGGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGa                                                                      <  7:256152/136‑1 (MQ=255)
                                                      cgcggttgttGGCCTTTATTAGCACCTAGGCCCCAAAGGGCCTGTTCGTCGGGGCATACGGGGGGGCAAGCCTAACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGAc                                                                     <  6:433049/131‑1 (MQ=255)
                                                       gcggttgtgGGCCTTTATTTTCACCGTTCCCCGCACCGGGCGGGGGGGGGGTGCCTCGGGGGGGCCCATGCGCCCCGAGGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACt                                                                    <  4:310879/132‑1 (MQ=255)
                                                                    cTGATTTTCCCCCAGCCCCGCACCGGCGGGGTCGGGGGTGCCTAGGGGGGGCCCATGCCCACCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGccc                                                      <  7:139457/136‑1 (MQ=255)
                                                                      gAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCCACGCAGACGGCCGGCGCCGCTGGCAGTCGGAGATCGGCCGCCTGCCGAACGGGCGTGGACAGCACCAGTATGcccag                                                    >  8:56308/1‑134 (MQ=255)
                                                                                   gCCCGCAACGGCCTGGGGGTCGGTGCCTACGGGGTGCCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                                        <  4:308224/135‑1 (MQ=255)
                                                                                   ccccGCACCGGCCGGGGCGGCGGGGCCTCGGGGGGGGCCAGCCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                                        <  1:155119/134‑1 (MQ=255)
                                                                                       gCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCa                                   >  8:41220/1‑136 (MQ=255)
                                                                                              cggggTCGGCGGGGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTa                             <  3:426025/132‑1 (MQ=255)
                                                                                                      gcgGGGCCTCCGGGTGGGCCCGCCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCg                     <  6:175763/134‑1 (MQ=255)
                                                                                                                    ggggCCAGGCGCCCGGAGGACGGCGGGCGCAGCTGGCAGTCGGAGAGCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCctg       <  3:198723/133‑1 (MQ=255)
                                                                                                                      ggCCCTGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGGGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCctgct     <  3:116278/135‑1 (MQ=255)
                                                                                                                         cATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTgctgct  >  2:156356/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                               
CAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCT  >  KT2440/3625989‑3626244

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: