Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1436,043,1960CA31.1% 84.2 / 19.4 61A37S (GCC→TCC) puuD_3Gamma‑glutamyl‑gamma‑aminobutyrate hydrolasePuuD
Reads supporting (aligned to +/- strand):  ref base C (29/13);  new base A (0/19);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.24e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  ACB143/6043067‑6043327
                                                                                                                                 |                                                                                                                                   
gggCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCCGGCGAGCCAGCAGCGCTTGCGGGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgcc                                                                                                                               <  7:201190/136‑1 (MQ=255)
  gcgcCATTGTAATGATGCGGCTCGATATATGAAGGCGAGCCCGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCAGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgacgg                                                                                                                             >  2:43470/1‑136 (MQ=255)
       aTTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGgccgcc                                                                                                                        >  7:312812/1‑136 (MQ=255)
        ttGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCt                                                                                                                       >  6:133682/1‑136 (MQ=255)
          gTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCg                                                                                                                     >  5:210010/1‑136 (MQ=255)
                   cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTg                                                                                                            >  7:129192/1‑136 (MQ=255)
                   cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTg                                                                                                            >  2:82481/1‑136 (MQ=255)
                      cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCg                                                                                                         >  4:5687/1‑136 (MQ=255)
                      cTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACAGACTTGTCg                                                                                                         >  5:51966/1‑136 (MQ=255)
                        cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGcc                                                                                                       >  7:108472/1‑136 (MQ=255)
                        cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGcc                                                                                                       >  3:94648/1‑136 (MQ=255)
                        cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGcc                                                                                                       >  1:312/1‑136 (MQ=255)
                                  aGGCGAGCCGTTAAAAACAATGCCTTGCAGGCGAGCCAGCCGCGCTTGCGTGTCGCGGCCGCCATCACGGGCCGGAAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGc                                                                                             <  4:265551/136‑1 (MQ=255)
                                     cgcgCCGGTAAAAACACTGCCGGCCCGGCGAGCCAGCAGCGCTGGCGGGTCGCGGCCGCCATCACGGGCGGGCAAGATCAGGGGAACGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgag                                                                                          <  5:300400/133‑1 (MQ=255)
                                     cgcgCCGGAAAAAACAAGCCCGGGCCGGCGAGCCAGCCGCGGTGGCGGGCCGCGCCCGCCATCACGGGCCGCCAAGATCAGGGCCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgag                                                                                          <  5:68927/133‑1 (MQ=255)
                                      gAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCGCGACATACTTGTCTCCGACCGTGTGCgag                                                                                          >  3:265159/1‑135 (MQ=255)
                                        gccgtaaaaaaCACTGCCGTGCTGGCGAGCCAGCCGCGCTTGCGGGGCGCTGCCGCCATCAGGGGCCGCCAAGATCTGTGCCACGCCGTAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgaggag                                                                                       <  3:72204/130‑1 (MQ=255)
                                          cGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGtt                                                                                     >  1:310629/1‑136 (MQ=255)
                                          cGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGtt                                                                                     >  5:32986/1‑136 (MQ=255)
                                          cGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGt                                                                                      >  2:174132/1‑135 (MQ=255)
                                                         cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACctgct                                                                       >  4:104502/1‑135 (MQ=255)
                                                            gCAGGCGAGCCAGCAGCGCTTGCGTGTCGCGGCCGCCATCACGGGCCGGCAATATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  1:69659/136‑1 (MQ=255)
                                                             cAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  6:138517/135‑1 (MQ=255)
                                                              aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  1:270382/1‑135 (MQ=255)
                                                              aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  7:133227/1‑135 (MQ=255)
                                                              aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  2:60154/1‑135 (MQ=255)
                                                                                 gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATTCTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCa                                               >  3:308652/1‑135 (MQ=255)
                                                                                 gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATcg                                               >  1:2775/1‑134 (MQ=255)
                                                                                    tgtgGCTGCCGCCATCACGGGCCGGCAATATCAGGGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAgg                                           <  3:48653/136‑1 (MQ=255)
                                                                                      tCGCTGCCGCAATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATa                                                                                                                >  4:133900/1‑65 (MQ=255)
                                                                                      tCGCTGCCGCAATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATa                                                                                                                <  3:133900/65‑1 (MQ=255)
                                                                                         cTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCg                                      >  5:115659/1‑136 (MQ=255)
                                                                                              gccATCACGGGCCGGCAAGATCAGGGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCAttt                                 <  5:345992/136‑1 (MQ=255)
                                                                                                  tcccgggcgggAAAGATCTGGGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:200878/133‑1 (MQ=255)
                                                                                                  tcccGGGCCGGCAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:124672/133‑1 (MQ=255)
                                                                                                  tCAGGGGCGGGCAAGTTCAGTGGCAAGCGGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:42102/136‑1 (MQ=255)
                                                                                                  tCAGGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  1:269518/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGTCAAGATCTGGGGCAAGCGGTAAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:310629/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:324102/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCACGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:257807/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACCTACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:98093/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:92038/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:26921/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:191661/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:204172/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  4:290975/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:104502/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:120195/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:172789/136‑1 (MQ=255)
                                                                                                  tCACGGGCCGGAAATATCAGTGGAAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGGCGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:119671/136‑1 (MQ=255)
                                                                                                  cccggggCCGCAAAGATCTGGGGAAAGCCGCCAAAGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGGTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  3:348898/132‑1 (MQ=255)
                                                                                                    aggggCCGGAAATATCAGTGGAACGCCGGAAATGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCt                            <  3:50641/133‑1 (MQ=255)
                                                                                                         ccGGCAATATCAGTGGCATGCCGCAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATggg                      <  4:302773/136‑1 (MQ=255)
                                                                                                            gCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCtt                   <  4:94648/136‑1 (MQ=255)
                                                                                                               aGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTaca                >  3:143053/1‑136 (MQ=255)
                                                                                                                    aGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgtg           >  4:137868/1‑136 (MQ=255)
                                                                                                                         cAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAg      >  6:276923/1‑136 (MQ=255)
                                                                                                                          aaGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAgg     <  4:332953/136‑1 (MQ=255)
                                                                                                                             ccGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  >  2:45667/1‑136 (MQ=255)
                                                                                                                              cGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  >  6:131973/1‑135 (MQ=255)
                                                                                                                              cGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  >  5:279635/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                   
GGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  ACB143/6043067‑6043327

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