Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1226,038,6580CA29.8% 47.8 / 16.2 47A37S (GCC→TCC) puuD_3Gamma‑glutamyl‑gamma‑aminobutyrate hydrolasePuuD
Reads supporting (aligned to +/- strand):  ref base C (21/12);  new base A (0/14);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 4.41e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  ACB122/6038525‑6038789
                                                                                                                                      |                                                                                                                                    
cTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCa                                                                                                                                     >  2:164618/1‑135 (MQ=255)
      gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCTTGAAGGCGAGCCAGCACCGCTTGCGTGTCGCTGCCGCCAACAAGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑gccgg                                                                                                                             >  8:224079/1‑136 (MQ=255)
         ccATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑gccggccg                                                                                                                          >  7:20444/1‑136 (MQ=255)
                  tgatgCGGCGCGATATTTCACGGCGAGCCGGTAAAAACAATGCCGGGCAGGCGAGCCAGCAGCGCTTGCGTGGCGCGGCCGCCATCCCGGGCC‑GGAAAGATCAGGGGCACGCCGGCAAA‑GCCGGCCGCCTCGACa                                                                                                                  <  6:114379/135‑1 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACtt                                                                                                             >  6:115544/1‑135 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTg                                                                                                            >  3:58287/1‑136 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGAAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACtt                                                                                                             >  6:251830/1‑135 (MQ=255)
                            ctatatTTCACGGCCCGCCGGTAAAAACACCGCCGTGCCGGCGAGCCAGCCGCGCTGGCGGGGCGCGGCCGCCATCCCGGGCC‑GGCAAGATCAGGGGCAAGCCGGAAAA‑GCCGGCCGCCTCGACATACTTGTCGcc                                                                                                       <  3:187407/134‑1 (MQ=255)
                            cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGcc                                                                                                       >  1:190121/1‑136 (MQ=255)
                               tttttGACGGCGCGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGGCGCGGCCGCCATCACGGGCC‑GGCAAGATCAGTGCCAAGCCGGAAAA‑GCCGGCCGCCTCGACATACTTGTCGcc                                                                                                       <  8:20444/131‑1 (MQ=255)
                                      aGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACAGTGTGc                                                                                             >  8:226061/1‑136 (MQ=255)
                                            gCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgaggag                                                                                       >  8:153699/1‑136 (MQ=255)
                                                                gCAGGCGAGCCCGCCGCGCTGGCGTGGCGCTGCCGCCAGCCCGGGCC‑GGAAAGAGCCGGGCCACGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  4:61494/136‑1 (MQ=255)
                                                                ccAGGCGAGCCCGCAGCGCTTGCGTGGCGCGGCCGCCATCCCGGGCC‑GGAAAGACCTGTGCCAAGCCGGAAAC‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                   <  4:113581/135‑1 (MQ=255)
                                                                  aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  1:9355/1‑135 (MQ=255)
                                                                  aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACgg                                                                  >  3:161270/1‑135 (MQ=255)
                                                                  aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGc                                                                 >  2:151989/1‑136 (MQ=255)
                                                                  aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGGGTGCGAGGAGTTGTTTGCTACCTGCTGACgg                                                                  >  5:180633/1‑135 (MQ=255)
                                                                             caggGCTTGGGTGGGGCTGCCGCAATCCGGGGCC‑GGAAATATATGTGCCAAGCCGGAAAA‑GGCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa                                                        <  4:108131/134‑1 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGCCCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCa                                               >  2:202438/1‑135 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCa                                               >  1:18436/1‑135 (MQ=255)
                                                                                             cTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGCGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCg                                      >  1:240927/1‑136 (MQ=255)
                                                                                             cTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGCCACCGCTCAAAGGGAcc                                       >  5:72773/1‑135 (MQ=255)
                                                                                               gccgccATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCa                                    >  8:73422/1‑136 (MQ=255)
                                                                                               gccgccATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCa                                    >  8:76832/1‑136 (MQ=255)
                                                                                                ccgcaATAAGGGGCC‑GGAAAGATCTGTGGCACCCCGGAAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCAt                                   <  2:9355/136‑1 (MQ=255)
                                                                                                   cctcctcGGGCC‑CGCATCTTCCGTGCCAACCCCGCACAGGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCAttt                                 <  7:89426/130‑1 (MQ=255)
                                                                                                   ccATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGCCCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTg                                >  8:158362/1‑136 (MQ=255)
                                                                                                      tccggggCG‑GGAAAGATCAGTGCAATGCCGGAAAA‑GCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:190121/132‑1 (MQ=255)
                                                                                                      tcccGGGCC‑GTCAAGATCAGTGCCACGCCGGAAAA‑GCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  8:185763/133‑1 (MQ=255)
                                                                                                      tcccGGGCC‑GGCAACTCCTGGGGCAAGCCGCAAAA‑GCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  1:164618/133‑1 (MQ=255)
                                                                                                      tCAGGGCCC‑GGCATGATGAGTGGCAAGCCGGAAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  2:152400/136‑1 (MQ=255)
                                                                                                      tCACGGGCC‑GTCAATATCATTGGCAAGCCGGAAAA‑GCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:226061/136‑1 (MQ=255)
                                                                                                      tCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:76832/136‑1 (MQ=255)
                                                                                                      tCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:73422/136‑1 (MQ=255)
                                                                                                      tCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  4:102788/136‑1 (MQ=255)
                                                                                                      tCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  4:161270/136‑1 (MQ=255)
                                                                                                      ccACGGGCG‑GGCAACATCCGGCGCAAGCCGCAAAA‑CCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  7:231746/135‑1 (MQ=255)
                                                                                                       ccggggCG‑GCCATGATCTGTGCCACGCCGGAAAT‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:115544/132‑1 (MQ=255)
                                                                                                       cccgggcg‑ggCAAGAACCGTGGCAAGCCGGCAAC‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:72773/133‑1 (MQ=255)
                                                                                                       caggggCC‑GGAAAGATCAGTGGCAAGCCGGAAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:109087/132‑1 (MQ=255)
                                                                                                       cACGGGCC‑GGCATGATCAGGGGAAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  6:96103/135‑1 (MQ=255)
                                                                                                        ccGGGCG‑GCCAAGATCCGGGCCACCCCGGCAAC‑GCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTc                           <  1:198550/135‑1 (MQ=255)
                                                                                                        ccGGGCCGGGAAAGATCTGGGGAACGCCGCAAAC‑GCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                             <  5:202298/134‑1 (MQ=255)
                                                                                                             gcggtcAAGTTCAGTGCAAACCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATgg                       <  8:156887/131‑1 (MQ=255)
                                                                                                                 tcAAGATCTGGGGAACGCCGCCAAC‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCtt                   <  2:222757/135‑1 (MQ=255)
                                                                                                                                   cGGCAAA‑GCCGGCCGCCTCGACATACTTGTAGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACCCCGATCAAAGGGCCCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAc  >  5:139724/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCC‑GGCAAGATCAGTGGCAAGCCGGCAAA‑GCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAC  >  ACB122/6038525‑6038789

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