Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24406,039,8370CA30.5% 64.0 / 18.5 59A37S (GCC→TCC) puuD_3Gamma‑glutamyl‑gamma‑aminobutyrate hydrolasePuuD
Reads supporting (aligned to +/- strand):  ref base C (30/11);  new base A (0/18);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 5.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  KT2440/6039710‑6039968
                                                                                                                                |                                                                                                                                   
gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAgccg                                                                                                                              >  2:86202/1‑135 (MQ=255)
        gTAATGATGCGGCTCTATATTTGACGGCGAGCCGGTAAAAACAATGCGGGGCCGGCGAGCCAGCAGCGCTTGCGGGGCGCTGCCGCCATCCCGGGCCGGCAAGAGCAGGGGCAA‑GCCGGAAAAGCCGGCCGCCt                                                                                                                       <  6:233778/134‑1 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACtt                                                                                                             >  8:20018/1‑135 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACtt                                                                                                             >  3:244980/1‑135 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACtt                                                                                                             >  4:147252/1‑135 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACtt                                                                                                             >  4:54854/1‑135 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACtt                                                                                                             >  1:70599/1‑135 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACtt                                                                                                             >  7:192056/1‑135 (MQ=255)
                 cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCACGCCATAcat                                                                                                             >  2:133396/1‑133 (MQ=255)
                   gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAC‑GCCGGCAAAGCCGGCCGCCGCGACATACATGt                                                                                                           >  7:281677/1‑135 (MQ=255)
                   gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGCCATACTTGt                                                                                                           >  8:252749/1‑135 (MQ=255)
                   gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGt                                                                                                           >  3:42221/1‑135 (MQ=255)
                   gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGt                                                                                                           >  4:27151/1‑135 (MQ=255)
                   gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGt                                                                                                           >  5:218430/1‑135 (MQ=255)
                   gCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGt                                                                                                           >  2:161002/1‑135 (MQ=255)
                    cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTc                                                                                                          >  1:222108/1‑135 (MQ=255)
                    cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCGCGCCCTACGTgga                                                                                                          >  2:149398/1‑133 (MQ=255)
                                 ggCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGtgtg                                                                                              >  2:14829/1‑134 (MQ=255)
                                       ccGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgaggag                                                                                       >  7:271124/1‑135 (MQ=255)
                                        cGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGt                                                                                      >  1:181606/1‑135 (MQ=255)
                                                    tcccGTGCAGGCGAGCCAGCAGCGGTTGCGTGCCGCGGCCGCAATCAGGGGCCGGCAATATCTGGGGCAA‑GCCGCAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACct                                                                          <  2:179202/133‑1 (MQ=255)
                                                       cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGCCATACTTGTCGCCGCCCGTGTGCGAGGCGTTCTTCCCTACcagct                                                                       >  2:286729/1‑135 (MQ=255)
                                                       cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGCGGAGTTCTTGCCTACctgct                                                                       >  5:10875/1‑135 (MQ=255)
                                                       cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACctgc                                                                        >  7:64729/1‑134 (MQ=255)
                                                           ccggCGAGCCAGCAGCGCTTGCGGGGCGCGGCCGCCATCCCGGGCCGGCAAGATCAGTGGCAA‑GCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  7:69793/133‑1 (MQ=255)
                                                           cAGGGGCGCCTGCTGCGCTTGGGTGGGGCGGCCGCCAGCCCGGGCGGGAATGACCCGGGGCAC‑CCGGGCAACGCGGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  8:9277/135‑1 (MQ=255)
                                                            aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAC‑GCCGGCAAAGCCGGCCGCCGCGACCTACGTGTCGCCGCCCGTGTGCGAGGCGTTCTTCCCTACCTGCGGCCgg                                                                  >  6:17339/1‑135 (MQ=255)
                                                            aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGCCATACTTGTCGCCGCCCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  1:168476/1‑135 (MQ=255)
                                                                gAGCCCGCAGGGGTGGCGGGGCGCGGCCGCCATCCCGGGCCGGCAAGATCAGTGCCAC‑GCGGCAAAAGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAgg                                                              <  6:218430/135‑1 (MQ=255)
                                                                gAGCCAGCAGCGGTGGCGTGGCGCGGCCGCCAGCACGGGCCGGAAAGATCAGTGCCAC‑GCCGGCAACGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAgg                                                              <  6:139691/135‑1 (MQ=255)
                                                                    cagcagCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCt                                                          >  3:189551/1‑135 (MQ=255)
                                                                     agcagcGCTTGCGTGGCGCTGCCGCCATCACGGGCCGGCAAGATCTGTGGCAC‑GCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTg                                                         <  6:50929/135‑1 (MQ=255)
                                                                      gcagcGCTTGGGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAC‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa                                                        <  8:211146/135‑1 (MQ=255)
                                                                               gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGCGTTCTTCCCTACCTGCTGACGGCAGGCGCTGCCACCGATc                                                >  1:87545/1‑134 (MQ=255)
                                                                                   gtCGTTGCCTCCATCAGGGGCCGGCATTTTATGTGGAAA‑GCCTCAAATGCGGGCCGCCTCGACATACTTGTAGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAgg                                           <  2:40817/135‑1 (MQ=255)
                                                                                       cTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGCCGGCAGGCGCTGACACCGATCAAAGGGccc                                       >  5:7915/1‑135 (MQ=255)
                                                                                       cTGCCGCCATCACGGGCCGGCAAGATCACTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGCCCGTGGGCGAGGAGTTCTTCCCTACCTGCAGACGGCAGGCGCTGCCACCGATCAAACGGccc                                       >  1:89616/1‑135 (MQ=255)
                                                                                                 ctcgggcgggCAAGTTCTGTGCCAA‑GCGGGAAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  4:42221/133‑1 (MQ=255)
                                                                                                 ccggggCGGCCATGATCGGGGGCAC‑GCCGGAAATGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:20018/132‑1 (MQ=255)
                                                                                                 ccggggCCGCCAAGAGCCGGGCCAC‑GCCGCAAACGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:95149/132‑1 (MQ=255)
                                                                                                 ccggggCCGCAAACATCTGTGCCAC‑GCCGCAAACGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  1:186/132‑1 (MQ=255)
                                                                                                 cccGGGCCGGCAATATCCGTGCCAA‑GCCGGAAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:181606/133‑1 (MQ=255)
                                                                                                 cccGGGCCGGCAAGATCTGTGGCAC‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:168476/133‑1 (MQ=255)
                                                                                                 cccGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:232890/133‑1 (MQ=255)
                                                                                                 cccGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  1:38987/133‑1 (MQ=255)
                                                                                                 cccGGGCCGGAAATATCCG‑GGCAACGCCGCAAATGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:239774/133‑1 (MQ=255)
                                                                                                 cccGGGCCGGAAAGATCCGTGCCAC‑GCCGGAAAAGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:192056/133‑1 (MQ=255)
                                                                                                 cccGGGCCGGAAAGATCAGGGGCAC‑GCCGCAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:27151/133‑1 (MQ=255)
                                                                                                 cccGGGCCGCCACGAGCTGGGCCAA‑GCGGGAAACGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:42518/133‑1 (MQ=255)
                                                                                                 cACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:64729/135‑1 (MQ=255)
                                                                                                 cACGGGCCGGCAAGATCAGTGGCAA‑GCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:54854/135‑1 (MQ=255)
                                                                                                 cACGGGCCGGCAAGATCAGTGGCAA‑CCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:6471/135‑1 (MQ=255)
                                                                                                 cACGGGCCGGCAAGATCAGTGGAAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:217411/135‑1 (MQ=255)
                                                                                                 cACGGGCCGGCAAGATCAGTGCCAC‑GCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:150736/135‑1 (MQ=255)
                                                                                                         ggCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCAGACCGCCGGCGCTGACACCGATCAAATGGACCCCATTTGCCCTCATGGGc                     >  1:133312/1‑135 (MQ=255)
                                                                                                          gCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTATTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGc                     >  3:85998/1‑134 (MQ=255)
                                                                                                                   gTGGCAA‑GCGGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgtg           <  2:198118/135‑1 (MQ=255)
                                                                                                                    tGGCAA‑GCCGGCAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgtgt          <  7:210254/135‑1 (MQ=255)
                                                                                                                    gggCAT‑GCGGGAAACGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgtgt          <  8:40710/134‑1 (MQ=255)
                                                                                                                             ggccAACGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  <  1:161002/132‑1 (MQ=255)
                                                                                                                                |                                                                                                                                   
GCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAA‑GCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  KT2440/6039710‑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: ‑