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

CCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGC  >  REL606/1795578‑1795725
                                                                        |                                                                           
ccAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacca                                                                          <  2:1165945/76‑1 (MQ=255)
 cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg                                                                         >  2:824559/1‑76 (MQ=255)
     gTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCAt                                                                     <  1:889442/76‑1 (MQ=255)
       tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATaa                                                                   >  1:928577/1‑76 (MQ=255)
       tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATaa                                                                   >  1:1449918/1‑76 (MQ=255)
          caatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTc                                                                >  1:1669607/1‑76 (MQ=255)
          caatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTc                                                                >  1:1834548/1‑76 (MQ=255)
            atcatcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGc                                                              >  2:1055852/1‑76 (MQ=255)
            atcatcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGc                                                              >  2:821374/1‑76 (MQ=255)
             tcatcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCt                                                             >  2:94394/1‑76 (MQ=255)
             tcatcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCt                                                             >  1:818367/1‑76 (MQ=255)
              cagccGTTGGCAGAGATCGAAAATAATACCCAGTTTCGCTTCACCACCTTCGACACCCCCCAGCCATAATTCGCTg                                                            <  2:698872/73‑1 (MQ=255)
                   gTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCgg                                                       >  2:1986374/1‑76 (MQ=255)
                    ttGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGc                                                      >  2:839356/1‑76 (MQ=255)
                     tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCa                                                     <  2:592443/76‑1 (MQ=255)
                     tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                                                     >  2:502993/1‑76 (MQ=255)
                     tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                                                     >  1:859694/1‑76 (MQ=255)
                       gCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAt                                                   >  2:1021737/1‑76 (MQ=255)
                        cAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAtt                                                  >  1:672435/1‑76 (MQ=255)
                        cAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAtt                                                  >  2:1809468/1‑76 (MQ=255)
                        cAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAtt                                                  >  2:1211761/1‑76 (MQ=255)
                        cAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAtt                                                  >  1:1265047/1‑76 (MQ=255)
                            gaTCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTg                                              <  1:1358801/76‑1 (MQ=255)
                            gaTCAAATATAATACCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTg                                              <  2:813741/76‑1 (MQ=255)
                             aTCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGa                                             >  2:1966848/1‑76 (MQ=255)
                               caaaaaaaATATCCGGTTTCGCTTTTCCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAGTCGTGATa                                           <  2:90496/74‑1 (MQ=255)
                               cGAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                           <  1:26086/76‑1 (MQ=255)
                               cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                           >  2:982917/1‑76 (MQ=255)
                               cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa                                           >  1:737878/1‑76 (MQ=255)
                                gAAAATAATATCCTGTTTCGCTTTCCCCCCTGAGCCAACCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATaa                                          <  2:1408463/76‑1 (MQ=255)
                                gAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATaa                                          <  1:1677775/76‑1 (MQ=255)
                                  aaaTAATATCCAGTTTCCTTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg                                        <  1:188886/76‑1 (MQ=255)
                                  aaaTAATACCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg                                        <  2:1512860/76‑1 (MQ=255)
                                   aataatATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGc                                       >  1:1387253/1‑76 (MQ=255)
                                     taataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt                                     >  1:713880/1‑76 (MQ=255)
                                     taataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt                                     >  1:1952628/1‑76 (MQ=255)
                                      aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta                                    >  2:1156569/1‑76 (MQ=255)
                                       atatCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa                                   <  2:1072157/76‑1 (MQ=255)
                                        tatCCAGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat                                  <  2:149011/76‑1 (MQ=255)
                                        tatCCAGTTTCGCTTTCCCACCTGCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat                                  <  1:48360/76‑1 (MQ=255)
                                        tatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat                                  >  2:1691512/1‑76 (MQ=255)
                                        tatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat                                  >  2:960692/1‑76 (MQ=255)
                                        tatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat                                  >  1:303561/1‑76 (MQ=255)
                                        tacccATTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat                                  <  1:1048159/73‑1 (MQ=255)
                                         atCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata                                 >  1:1693242/1‑76 (MQ=255)
                                          cccAGTTCCCCTTTCCCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa                                <  2:1189765/75‑1 (MQ=255)
                                            cAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAAt                              >  2:1740428/1‑76 (MQ=255)
                                             aGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATc                             <  1:295347/76‑1 (MQ=255)
                                                ttCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt                          <  2:1238484/76‑1 (MQ=255)
                                                ttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt                          >  2:468112/1‑76 (MQ=255)
                                                  cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc                        >  2:319352/1‑76 (MQ=255)
                                                   ccTTTCCCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCa                       <  2:1025157/75‑1 (MQ=255)
                                                     tttCCCCCCTCCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAgn                     <  1:622679/76‑2 (MQ=255)
                                                      ttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa                    >  2:841877/1‑76 (MQ=255)
                                                       tACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTaa                   >  1:977956/1‑76 (MQ=255)
                                                       tACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTaa                   >  1:409953/1‑76 (MQ=255)
                                                         cccccTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACt                 <  2:70796/73‑1 (MQ=255)
                                                         cccccTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACt                 <  2:1243380/73‑1 (MQ=255)
                                                         ccaccTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACt                 >  1:1813050/1‑76 (MQ=255)
                                                          caccTTCGACACCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTc                <  1:1872753/76‑1 (MQ=255)
                                                           accTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc               >  1:769794/1‑76 (MQ=255)
                                                           accTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc               >  2:365515/1‑76 (MQ=255)
                                                           accTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTAGTGATAACGCTTAATAAATCCTTCCAGTAACTcc               >  1:1169400/1‑76 (MQ=255)
                                                             cTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgg             >  1:965965/1‑76 (MQ=255)
                                                               tCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg           >  2:1600761/1‑76 (MQ=255)
                                                               tCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg           >  1:1586485/1‑76 (MQ=255)
                                                               tCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg           >  2:1084864/1‑76 (MQ=255)
                                                                caacacCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGt          <  1:632095/74‑1 (MQ=255)
                                                                 gACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGta         >  1:857585/1‑76 (MQ=255)
                                                                     tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc     >  1:2015460/1‑76 (MQ=255)
                                                                     cccaccaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc     <  1:354855/75‑1 (MQ=255)
                                                                        accaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgctgc  >  1:1032528/1‑76 (MQ=255)
                                                                        |                                                                           
CCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGC  >  REL606/1795578‑1795725

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