Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,866,8860TG30.2% 56.4 / 14.1 43N188T (AAC→ACC) aqpZAquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (7/23);  new base G (13/0);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 2.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGC  >  KT2440/4866750‑4867016
                                                                                                                                        |                                                                                                                                  
ttAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGt                                                                                                                                    <  1:239691/137‑1 (MQ=255)
   aGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGa                                                                                                                                 <  1:93512/137‑1 (MQ=255)
   aGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGa                                                                                                                                 <  3:76718/137‑1 (MQ=255)
        cTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACagag                                                                                                                            <  2:98880/137‑1 (MQ=255)
           tCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTa                                                                                                                         <  6:177678/137‑1 (MQ=255)
           tCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCGACCCGCCGGGGTGCCacagtgt                                                                                                                         >  6:26189/1‑134 (MQ=255)
                aGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGt                                                                                                                    <  5:5747/137‑1 (MQ=255)
                      cGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCgg                                                                                                              <  8:203834/137‑1 (MQ=255)
                      cGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCgg                                                                                                              <  8:200427/137‑1 (MQ=255)
                       gATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCggg                                                                                                             <  5:351381/137‑1 (MQ=255)
                                cgccgcCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATc                                                                                                    <  1:218101/137‑1 (MQ=255)
                                cgccgcCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGGTACGCCCGGGGGGGGCCGAGGGGTTGGTGACGGGGGTCgaggtg                                                                                                    >  7:308704/1‑134 (MQ=255)
                                 gccgccGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGGCCCCGGCGTCGCGGGGGGGGGGGAGAGGGTTTGGGGAACGGGATGGGGATCa                                                                                                   >  3:325974/1‑137 (MQ=255)
                                 gccgccGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGGTCCGCGCGGGGTTGACAGAGGTATTGGTGACCGGGATCGAGAACa                                                                                                   >  2:218101/1‑137 (MQ=255)
                                 gccgccGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGGTGTCAGAGGTATTGGGGACCGGGTTCGAGATCa                                                                                                   >  2:46798/1‑137 (MQ=255)
                                     ccgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                               <  4:160181/137‑1 (MQ=255)
                                     ccgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                               <  4:17285/137‑1 (MQ=255)
                                     ccgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGGCACCGGGATCGAGATCAGGTg                                                                                               <  4:189570/137‑1 (MQ=255)
                                      cgATAGCAGCACCGATCTGCGGCGCCACCCAGAACAGCCTCAGCTGCTGCAGGGCCCAGCGACGGTGGAAGAGGGCCGGCCCCGTGCTACGCGCGGGGTTGAGAGTGGTATTGGTCACCGGGATCTAGATCAGGTgg                                                                                              <  5:182052/137‑1 (MQ=255)
                                      cgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGGCCCCGGGGACCGGCGGGGGTGGGAGCGGGTTTGGTGGCCGGGGTCGAGCTGAGgggg                                                                                              >  2:147635/1‑137 (MQ=255)
                                                 ccGCTCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTTCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcag                                                                                   <  8:194798/137‑1 (MQ=255)
                                                 ccGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTAGGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcag                                                                                   <  3:130575/137‑1 (MQ=255)
                                                 ccGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagggtcag                                                                                   <  6:174406/137‑1 (MQ=255)
                                                     tCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGGTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGTTCAGGGTCAGGGcc                                                                               >  8:340110/1‑137 (MQ=255)
                                                       aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGGCCCCGCCGCCGCCCGGGGTGGTGAGAGGGTTTGGTGGCCCGGGGAGCGATCACGTGGGtaaggggcaggggccc                                                                             >  6:45321/1‑133 (MQ=255)
                                                       aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGGGCTACGCGCGGGGGTGTGAGAGGGATTGGTGCCCGGGCTCGGCATCAGGTGGGtgagggtccggggccg                                                                             >  4:217561/1‑133 (MQ=255)
                                                       aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGGGCGAACCGCGGGGGTGGGAGAGGTATTTGTGGACGGGCTCGCGTTGATGTGGGTGAGGGGGAGGGcccc                                                                             >  2:152573/1‑135 (MQ=255)
                                                       aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCTTCGCGCGGGGGTGTCCGCGGTTTTGGTGCCCGGGCGAGAGGTCATGGGGGtcaggggcaggggccg                                                                             >  8:5884/1‑133 (MQ=255)
                                                       aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCCGCGACGCGCGGGGTGGTCACAGGTATTGGGGGCCGGGATCGGGATGACGTGGATGAGGGGCAGGCCCAg                                                                             >  3:114393/1‑137 (MQ=255)
                                                                 cccAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCcgatggcg                                                                   <  6:201537/137‑1 (MQ=255)
                                                                               aGCTGCTGCAGGGCCCAGCCACCGCCGACCAGGGCCGGCCCCGTGCGACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAAcc                                                     <  2:274753/137‑1 (MQ=255)
                                                                                   gctgcAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCTGCTACGCGCGGGGGTGACAGAGGGATTGGGCCCCGGGGGAGCGATCACGGGGAGCCGGGTGAGGGACAGGCCGCAGGGGCTGGGTGGAAAACCAGcc                                                 >  8:341725/1‑137 (MQ=255)
                                                                                      gcAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTTTTTGTCACCGGGGTAGAGATCAGGGGGCTGAGGGTCAGGGCCCCGCCGCTGGCGATTGGTGCAAAACCAGcccgg                                              >  8:258870/1‑134 (MQ=255)
                                                                                      gcAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGGCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGt                                              <  5:13205/137‑1 (MQ=255)
                                                                                       cAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGGGGGACGCGCGGGGGTGACAGAGGTTTTGGTCACCGGGAGCGCGTGCAGGTGGGTGAGGGGCCGGGCCCGGCCGCTGGGGGTGGGTGGAAAACCTGCCGGTg                                             >  4:267740/1‑137 (MQ=255)
                                                                                        aGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCTGCTACGCGCGGGGGTGGCAGAGGGATTTGGGACCCGGGTGGAGATCAGGTTGGTCCGGGGGAGGGCCAGGCCGATTGCGATGGGTGCAAAACCTTCCCGTgc                                            >  6:69403/1‑137 (MQ=255)
                                                                                               aGCCACCGACGAACAGGGCCGGCCCCGGGCGTAGCGCGGGGGTGGCACAGGGATTTGGGCCCGGGGTCGCGATGAGGGGGATGAGGGGCAGGGCCAGGCCGCTGGCGGTGGGGGCAAAAACTGCCGGGGCGCCGGCg                                     >  2:21033/1‑137 (MQ=255)
                                                                                                 ccACCGACGAACAGGGCCGGCCCCGTGGTTCGCGCGGGGTTTACACAGGTATTGGTGAACGGGATCGAGATGAGGGGGGTACGGGGCACGGCCCCGCCGCTGGCGCTGGGGGCTAAAACGGCCGGGGCGCGGGCggg                                   >  2:103482/1‑135 (MQ=255)
                                                                                                     cgacgaACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCgtcggtc                               <  8:347302/137‑1 (MQ=255)
                                                                                                               ggccggccCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGTGATCGAGATCAGGTGGATCAGGGTCAGGGCCACGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatga                     <  4:325974/137‑1 (MQ=255)
                                                                                                                  cggccCCGAGCTACGCGCGGGGTTGACACAGGTATTGGTACACGGGGTGCAGATGAGGTGGATAAGGGTCACGGCCACGCCGAAGGCGATGGGTGCAAAAACTTCCGGTGCGCGGGCGGCGGTCGCCCCCAtgaaga                  >  8:278554/1‑137 (MQ=255)
                                                                                                                    gccCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGacc                <  1:209361/137‑1 (MQ=255)
                                                                                                                                gcgctGGGTTGTCAGTGGTATTGGGCTCCGGGATCGAGATCAGGTTGATCAGGGGCAGGTCCAGGCCGATGGCTTTGTGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATc    <  6:157958/137‑1 (MQ=255)
                                                                                                                                  gcgGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGc  >  3:18327/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                  
TTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGTCAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGC  >  KT2440/4866750‑4867016

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