Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,795,6500AC29.0% 150.7 / 11.7 62V39G (GTG→GGG) celFcryptic phospho‑beta‑glucosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (42/2);  new base C (0/18);  total (42/20)
Fisher's exact test for biased strand distribution p-value = 1.03e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCT  >  REL606/1795575‑1795723
                                                                           |                                                                         
aCGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATcca                                                                           >  1:1291116/1‑76 (MQ=255)
 cGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccac                                                                          <  1:230438/76‑1 (MQ=255)
    cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg                                                                       >  1:1139207/1‑76 (MQ=255)
     aGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGc                                                                      >  1:1037551/1‑76 (MQ=255)
     aGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGc                                                                      >  1:432213/1‑76 (MQ=255)
      gCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACTACCTTCGACATCCACCAGcc                                                                     >  1:1056336/1‑76 (MQ=255)
       cGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCa                                                                    >  1:839063/1‑76 (MQ=255)
        gTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCAt                                                                   >  2:261673/1‑76 (MQ=255)
          tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATaa                                                                 >  2:64277/1‑76 (MQ=255)
           atcaatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAAt                                                                >  2:1853935/1‑76 (MQ=255)
           atcaatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAAt                                                                >  1:975350/1‑76 (MQ=255)
           atcaatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAAt                                                                >  2:1813815/1‑76 (MQ=255)
            tcaatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAAtt                                                               >  2:1264272/1‑76 (MQ=255)
             caatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTc                                                              >  2:880158/1‑76 (MQ=255)
             caatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTc                                                              >  2:395081/1‑76 (MQ=255)
               atcatcCGTTGGCAGAGATCGAAAATAATAACCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTcac                                                            >  1:731084/1‑74 (MQ=255)
                tcatcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCt                                                           >  2:536686/1‑76 (MQ=255)
                 catcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTg                                                          >  1:12165/1‑76 (MQ=255)
                      gTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACATCCACCAGCCATAATTCGCTGACCgg                                                     <  2:1479473/76‑1 (MQ=255)
                      gTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCgg                                                     >  2:1544990/1‑76 (MQ=255)
                      gTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCgg                                                     >  1:924465/1‑76 (MQ=255)
                       ttGGCAGAGATCGAAAATAATATCCAGTTTCGATTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGc                                                    <  1:77836/76‑1 (MQ=255)
                         ggCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCaa                                                  >  1:1681918/1‑76 (MQ=255)
                            agagATCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTc                                               <  1:1615583/76‑1 (MQ=255)
                            agagATCAAAAATAATATCCAGTTTCGCTTTCCCCCCTGCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTc                                               <  1:985575/76‑1 (MQ=255)
                               gaTCGAAAATAATATCCATTTTCGCTTTCCCCCCTCCGACACCCACCAGCCATAATTCGCTGACCGGCAATTCGTg                                            <  2:1081912/76‑1 (MQ=255)
                                  caaaaaTAATATCCAGTTTCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                         <  1:1108439/74‑1 (MQ=255)
                                        taataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt                                   >  2:7039/1‑76 (MQ=255)
                                         aataTCCAGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta                                  <  1:167082/76‑1 (MQ=255)
                                         aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta                                  >  1:508372/1‑76 (MQ=255)
                                         aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta                                  >  2:730522/1‑76 (MQ=255)
                                         aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta                                  >  1:988561/1‑76 (MQ=255)
                                         aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta                                  >  2:153857/1‑76 (MQ=255)
                                          atatCCATTTTCTCTTTCCCCCCTCCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa                                 <  1:370527/76‑1 (MQ=255)
                                          atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa                                 >  1:859823/1‑76 (MQ=255)
                                          atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACAGGCAATTCGTGATAACGCTtaa                                 >  2:85492/1‑76 (MQ=255)
                                            atCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata                               >  1:121892/1‑76 (MQ=255)
                                             tCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa                              >  1:1213269/1‑76 (MQ=255)
                                              ccAGTTTCGCTTTACCACCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa                             <  2:947789/76‑1 (MQ=255)
                                               cAGTTTCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAAt                            <  2:119760/76‑1 (MQ=255)
                                                aGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATc                           >  1:505019/1‑76 (MQ=255)
                                                 gTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATcc                          >  2:593091/1‑76 (MQ=255)
                                                  tttCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt                         <  2:1706963/76‑1 (MQ=255)
                                                  tttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt                         >  1:1147883/1‑76 (MQ=255)
                                                  tttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt                         >  2:1499215/1‑76 (MQ=255)
                                                  tttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt                         >  1:452890/1‑76 (MQ=255)
                                                  ttccGCTTTCCCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt                         <  1:262584/73‑1 (MQ=255)
                                                   ttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt                        >  1:73417/1‑76 (MQ=255)
                                                     cGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                      <  1:1060276/76‑1 (MQ=255)
                                                     cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                      >  1:101454/1‑76 (MQ=255)
                                                     cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                      >  1:1438420/1‑76 (MQ=255)
                                                      gCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCa                     >  1:1576936/1‑76 (MQ=255)
                                                        tttCCCCCCTTCCACCCCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt                   <  2:707366/76‑1 (MQ=255)
                                                        tttACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt                   <  1:16279/76‑1 (MQ=255)
                                                         ttACCCCCTTCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa                  <  2:4637/76‑1 (MQ=255)
                                                          tACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTaa                 <  2:1880472/76‑1 (MQ=255)
                                                           cccaccTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAAc                <  1:1081652/75‑1 (MQ=255)
                                                            cccccTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACt               <  1:1603508/73‑1 (MQ=255)
                                                             ccccTCCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTc              <  2:605170/74‑1 (MQ=255)
                                                              accTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc             >  1:1842124/1‑76 (MQ=255)
                                                               ccTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCg            <  2:1057460/76‑1 (MQ=255)
                                                                 ttCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCggg          >  2:224231/1‑76 (MQ=255)
                                                                 ttCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCggg          >  1:478775/1‑76 (MQ=255)
                                                                 ttCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCggg          >  2:1841591/1‑76 (MQ=255)
                                                                 ttCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCggg          >  1:279997/1‑76 (MQ=255)
                                                                  tCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg         <  1:54571/76‑1 (MQ=255)
                                                                   cGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGt        >  1:1727185/1‑76 (MQ=255)
                                                                         ccaccaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgct  >  2:1761308/1‑76 (MQ=255)
                                                                           |                                                                         
ACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCT  >  REL606/1795575‑1795723

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