Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,795,6500AC20.0% 125.7 / 10.3 50V39G (GTG→GGG) celFcryptic phospho‑beta‑glucosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (33/7);  new base C (0/10);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 1.89e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGC  >  REL606/1795577‑1795722
                                                                         |                                                                        
gCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacc                                                                        >  1:641794/1‑76 (MQ=255)
  cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg                                                                      >  2:1105422/1‑76 (MQ=255)
  cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGCTTCGCTTTACCACCTTCGACATCCACCAg                                                                      >  1:1128765/1‑76 (MQ=255)
     cGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCa                                                                   >  1:808994/1‑76 (MQ=255)
      gTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCAt                                                                  <  2:438752/76‑1 (MQ=255)
       ttATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATa                                                                 >  2:843420/1‑76 (MQ=255)
           caatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTc                                                             >  2:1426989/1‑76 (MQ=255)
                 tcCGTTGGCAGAGATCGAAAATAATATCCATTTTCGCTTTACCCCCTTCGACATCCACCAGCCATAATTCGCTGAc                                                       <  2:590678/76‑1 (MQ=255)
                   cGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCg                                                     >  1:1511616/1‑76 (MQ=255)
                     ttGGCAGAGATCGAAAATAATATCCATTTTCGCTTTACCACCTTCGACACCCACCAGCCATAATTCGCTGACCGGc                                                   <  1:1268154/76‑1 (MQ=255)
                      tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                                                  >  2:99099/1‑76 (MQ=255)
                      tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                                                  >  2:422341/1‑76 (MQ=255)
                      tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                                                  >  1:367995/1‑76 (MQ=255)
                          agagATCAAAAATAATATCCTGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTc                                              <  1:60666/76‑1 (MQ=255)
                           gagaTCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTcn                                             <  1:600270/76‑2 (MQ=255)
                                cGAAAATAATATCCATTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                        <  1:558486/76‑1 (MQ=255)
                                cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCACTTCGTGATa                                        >  2:1731489/1‑76 (MQ=255)
                                cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                        >  1:568324/1‑76 (MQ=255)
                                cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                        >  2:1366407/1‑76 (MQ=255)
                                 gAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa                                       >  2:505763/1‑76 (MQ=255)
                                 gAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa                                       >  1:1109430/1‑76 (MQ=255)
                                 gAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa                                       >  2:1713547/1‑76 (MQ=255)
                                  aaaaTAATATCCAGTTCCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAAc                                      <  1:1191839/76‑1 (MQ=255)
                                     ataataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCt                                   >  2:1033268/1‑76 (MQ=255)
                                      aaataTCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt                                  <  1:381126/75‑1 (MQ=255)
                                        atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa                                >  2:758978/1‑76 (MQ=255)
                                          atCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCttatta                              >  2:710910/1‑76 (MQ=255)
                                           tCCAGTTTCGCTTTACCACCTTCGACCCCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa                             <  1:1742217/76‑1 (MQ=255)
                                           tCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa                             >  2:1196867/1‑76 (MQ=255)
                                            ccAGTTTCGCTTTACCCCCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa                            <  2:1315583/76‑1 (MQ=255)
                                            ccAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa                            >  2:1338380/1‑76 (MQ=255)
                                            ccAGTTCCGCTTCACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa                            <  2:1687381/76‑1 (MQ=255)
                                             cAGTTCCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAAt                           <  1:934040/76‑1 (MQ=255)
                                                tttCGCTTTACCCCCTCCGACACCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt                        <  1:748044/76‑1 (MQ=255)
                                                 ttCGCTTTCCCCCCTCCACCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt                       <  1:678972/76‑1 (MQ=255)
                                                 ttCGCTTTCCACCCTTCGACACCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt                       <  2:1262477/76‑1 (MQ=255)
                                                   cGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                     <  1:1253762/76‑1 (MQ=255)
                                                   cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                     >  1:443390/1‑76 (MQ=255)
                                                   cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                     >  2:1554236/1‑76 (MQ=255)
                                                    gCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCa                    >  2:1248426/1‑76 (MQ=255)
                                                    gCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCa                    >  1:1165571/1‑76 (MQ=255)
                                                     cTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAg                   >  2:1250570/1‑76 (MQ=255)
                                                      tttCCCCCCTTCGCCATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt                  <  1:1109863/76‑1 (MQ=255)
                                                      tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt                  >  1:86057/1‑76 (MQ=255)
                                                       ttACCCCCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa                 <  1:658869/76‑1 (MQ=255)
                                                       ttACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa                 <  2:1611503/76‑1 (MQ=255)
                                                        tACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTaa                >  1:877021/1‑76 (MQ=255)
                                                         accaccTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAAc               >  1:235754/1‑76 (MQ=255)
                                                          cccccTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACt              <  1:1380632/73‑1 (MQ=255)
                                                           ccccTTCAACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTc             <  1:730336/74‑1 (MQ=255)
                                                                tCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg        >  1:465946/1‑76 (MQ=255)
                                                                tCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg        >  1:658341/1‑76 (MQ=255)
                                                                   ccttCCGCCAGCCGTAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtat     <  2:1162731/73‑1 (MQ=255)
                                                                   acacCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtat     <  2:1060249/76‑1 (MQ=255)
                                                                    cATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtata    >  2:1304320/1‑76 (MQ=255)
                                                                     ccccaccaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAg   <  2:1384752/74‑1 (MQ=255)
                                                                      tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc  >  2:582109/1‑76 (MQ=255)
                                                                      tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc  >  2:441434/1‑76 (MQ=255)
                                                                      tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAAAGCTTAATAAATCCTTCCAGTAACTCCGGGGGATAgc  >  1:715987/1‑76 (MQ=255)
                                                                         |                                                                        
GCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGC  >  REL606/1795577‑1795722

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