Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24406,039,8370CA33.3% 35.9 / 21.9 48A37S (GCC→TCC) puuD_3Gamma‑glutamyl‑gamma‑aminobutyrate hydrolasePuuD
Reads supporting (aligned to +/- strand):  ref base C (22/10);  new base A (0/16);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 4.36e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  KT2440/6039708‑6039968
                                                                                                                                 |                                                                                                                                    
gggCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgc                                                                                                                                 >  2:243622/1‑135 (MQ=255)
  gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccg                                                                                                                               >  4:126316/1‑135 (MQ=255)
   cgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccgg                                                                                                                              >  6:248115/1‑135 (MQ=255)
     ccATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccgg‑cc                                                                                                                           >  5:6254/1‑135 (MQ=255)
       aTTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGg‑ccgc                                                                                                                         >  6:255855/1‑135 (MQ=255)
        ttGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGg‑ccgcc                                                                                                                        >  1:147503/1‑135 (MQ=255)
                   cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACtt                                                                                                             >  8:51895/1‑135 (MQ=255)
                   cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACtt                                                                                                             >  8:236420/1‑135 (MQ=255)
                   cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACtt                                                                                                             >  7:78891/1‑135 (MQ=255)
                   cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACtt                                                                                                             >  4:108478/1‑135 (MQ=255)
                      cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTc                                                                                                          >  6:109127/1‑135 (MQ=255)
                      cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTc                                                                                                          >  3:81879/1‑135 (MQ=255)
                        cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCCAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGc                                                                                                        >  2:264104/1‑135 (MQ=255)
                          ttatTTGACGGCGCGGTGGTAAAAACATTGCCGTGCAGGGTACCCTGCAGCGCTGGCGTGTGGCTGACGCCATCACGGGCCGGAAAGATAAGGGTCAAGCCGTAAAAGCCGG‑CCGCCTCGACATACTTGTCGcc                                                                                                       <  8:268237/133‑1 (MQ=255)
                                      gAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCCGCGCTTGCGTGGCGCGGCCGCCATCACGGGCCGCCAAGATCAGTGGCAAGCCGGAAACGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCgag                                                                                          <  6:220503/135‑1 (MQ=255)
                                      gAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CAGCCTCGACATACTTGTCGCCGCCCGTGTGCgag                                                                                          >  6:207175/1‑135 (MQ=255)
                                          cGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGt                                                                                      >  2:168579/1‑135 (MQ=255)
                                                         cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACctgct                                                                       >  3:1763/1‑135 (MQ=255)
                                                         cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCAGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACctgc                                                                        >  8:236160/1‑134 (MQ=255)
                                                             cAGGCGTGCCAGCAGCGGTTGGGTGTCGCGGCCGCCATCCCGGGCGGGCAAGATCAGGGGCAACCCGGCAAAGCGGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  4:60605/135‑1 (MQ=255)
                                                             cAGGCGAGCCAGCAGCGCTTGCGTGGCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  8:44008/135‑1 (MQ=255)
                                                              aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  1:160309/1‑135 (MQ=255)
                                                              aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   >  5:63472/1‑134 (MQ=255)
                                                                        gcagcGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGCCCTACTTGTCGCCGACCGTGGGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa                                                        >  6:143494/1‑135 (MQ=255)
                                                                               gtGCGGGGCGCGGCCGCCATCCGGGGCGGCCAAGATCCGTGGCACGCCGGAAACGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGAt                                                 <  5:98254/134‑1 (MQ=255)
                                                                                                   ccggggCCGGCATGATCAGGGGCATCCCGGAAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:63472/132‑1 (MQ=255)
                                                                                                   ccggggCCGGCAAGATCTGGGCCAAGCCGCAAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGCCGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:252016/132‑1 (MQ=255)
                                                                                                   cccgggcgggCAAGTTCAGTGGCAAGGCGGAAAAGCGGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  4:117303/133‑1 (MQ=255)
                                                                                                   cccGGGGCGGCAAGATCTGTCTCAAGCCGTCAAAGCCCG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:83390/133‑1 (MQ=255)
                                                                                                   cccGGGCCGGAAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:27159/133‑1 (MQ=255)
                                                                                                   cccGGGCCGCCAAGACCTGTGCCACCCCGGAAACGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:88077/133‑1 (MQ=255)
                                                                                                   cccGGGCCGCAAACATCAGGGCCAAGCCGGAAACGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:23721/133‑1 (MQ=255)
                                                                                                   caggggCCGGAAAGAACAGTGCACAGCGGGAAATCCCGG‑ACGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:34948/132‑1 (MQ=255)
                                                                                                   cACGGGCGGGCAAGATCAGTGCAAAGCCGCCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:29947/135‑1 (MQ=255)
                                                                                                   cACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:236420/135‑1 (MQ=255)
                                                                                                   cACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:116281/135‑1 (MQ=255)
                                                                                                   cACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  4:63972/135‑1 (MQ=255)
                                                                                                   cACGGGCCGGCAAGATCAGGGGCAAGCCGGAAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:108162/135‑1 (MQ=255)
                                                                                                   cACGGGCCGGAAAGATCTGGGCCACGGCGGAAACCCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:250842/135‑1 (MQ=255)
                                                                                                     cGGGCCGGCAAGAGCTGTGGCACGCCGGAAACGCGGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTc                           <  3:66732/135‑1 (MQ=255)
                                                                                                     cGGGCCGGCAAGAGCAGTGGCACGCCGGAAACGCGGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTc                           <  5:119585/135‑1 (MQ=255)
                                                                                                                     ggggCACGCCGGAAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgt            <  8:186910/132‑1 (MQ=255)
                                                                                                                     ggggCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgt            <  7:284637/132‑1 (MQ=255)
                                                                                                                      gggCAAGCCGGAAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgtgt          <  4:181663/134‑1 (MQ=255)
                                                                                                                         cAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCa       >  8:209428/1‑135 (MQ=255)
                                                                                                                              cGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTa   >  4:174985/1‑134 (MQ=255)
                                                                                                                              cGGAAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  <  7:14710/135‑1 (MQ=255)
                                                                                                                               ggCCAACCCGGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  <  5:255855/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
GGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGG‑CCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  KT2440/6039708‑6039968

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