Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,688,9900AC38.5% 79.4 / 23.7 52T92P (ACC→CCC) nthDNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)
Reads supporting (aligned to +/- strand):  ref base A (23/9);  new base C (0/20);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 8.15e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCG  >  REL606/1688915‑1689065
                                                                           |                                                                           
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                             <  1:1161814/76‑1 (MQ=255)
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                             >  1:1639992/1‑76 (MQ=255)
 tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                            <  1:1287992/76‑1 (MQ=255)
    gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                         <  1:259503/76‑1 (MQ=255)
    gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                         <  1:1697219/76‑1 (MQ=255)
     cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                        >  1:1693326/1‑76 (MQ=255)
      gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                       >  2:551144/1‑76 (MQ=255)
       ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                      <  2:1531125/76‑1 (MQ=255)
       ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                      <  1:1470673/76‑1 (MQ=255)
            ggggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt                                                                 >  1:779078/1‑76 (MQ=255)
            ggggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt                                                                 >  1:947203/1‑76 (MQ=255)
               gTGAAAACCTATATAAAAACGATTGGGCTTTATAACAGCAAAGAAGAAAATATCATCAAACCCTGCCGTATCTTGc                                                              <  1:1092865/76‑1 (MQ=255)
                  aaacCCTATATCAAAACGATTGGGCTTTATAGCAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTgg                                                           <  2:2087225/76‑1 (MQ=255)
                  aaaaCCTATATCAAAACGATTGGGTTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTgg                                                           <  2:1170773/76‑1 (MQ=255)
                   aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa                                                          >  1:745502/1‑76 (MQ=255)
                     aCCTATATCAAAACAATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                        <  1:1534797/76‑1 (MQ=255)
                     aCCTATATAAAAACGATTGGGCTTTAAAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                        <  1:1689260/76‑1 (MQ=255)
                      ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca                                                       <  2:1155438/76‑1 (MQ=255)
                      ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                       >  2:372450/1‑76 (MQ=255)
                      ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                       <  2:1520840/76‑1 (MQ=255)
                         atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                    >  1:1024381/1‑76 (MQ=255)
                         atatCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca                                                    <  1:469600/76‑1 (MQ=255)
                          tatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAt                                                   >  2:817702/1‑76 (MQ=255)
                          tatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAt                                                   >  2:1375466/1‑76 (MQ=255)
                          tatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAt                                                   >  2:292515/1‑76 (MQ=255)
                          aaacAAAACGATTGGGCTTTATAACAGCAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                   <  1:1263185/73‑1 (MQ=255)
                            acAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATaa                                                 <  2:756739/75‑1 (MQ=255)
                              aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATg                                               <  2:779078/76‑1 (MQ=255)
                              aaaaCGATTGGGCTTTATAACAGAAAAGAAGAAAATATCATAAAAACCTGCCGTATCTTGCTGGAGCAGCATAATg                                               <  2:947203/76‑1 (MQ=255)
                               aaaCGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATgg                                              <  1:347570/76‑1 (MQ=255)
                                aaCGATTGGGTTTTAAAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGc                                             <  2:1239745/76‑1 (MQ=255)
                                 aCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                            >  1:884379/1‑76 (MQ=255)
                                 aCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                            >  1:1388830/1‑76 (MQ=255)
                                  cGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                           >  2:1669900/1‑76 (MQ=255)
                                   gATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                          >  2:1894656/1‑76 (MQ=255)
                                   gATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                          <  1:2015086/76‑1 (MQ=255)
                                   gATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                          >  2:501720/1‑76 (MQ=255)
                                      tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                       >  2:564478/1‑76 (MQ=255)
                                      tGGGCTTTAAAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                       <  1:1526096/76‑1 (MQ=255)
                                        ggCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc                                     <  2:138883/76‑1 (MQ=255)
                                          cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                   >  1:1057876/1‑76 (MQ=255)
                                          cTTTATAACAACAAAGCAGAAAATATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                   <  2:648408/76‑1 (MQ=255)
                                          cTTTAAAAAAGCAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                   <  1:232409/76‑1 (MQ=255)
                                           tttATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                                  <  1:1291238/76‑1 (MQ=255)
                                            ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                                 >  2:1000944/1‑76 (MQ=255)
                                            taataACAGCAAAGAAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                                 <  2:231336/74‑1 (MQ=255)
                                                  cAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGt                           <  2:869827/76‑1 (MQ=255)
                                                  cAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGt                           >  2:1348645/1‑76 (MQ=255)
                                                  cAGCAAAGAAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGt                           <  2:1480805/76‑1 (MQ=255)
                                                   aacAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtg                          <  2:832505/74‑1 (MQ=255)
                                                            aaaaaTATCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg                 <  1:556574/75‑1 (MQ=255)
                                                                 tatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc            >  1:1247812/1‑76 (MQ=255)
                                                                  atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc           >  2:1442884/1‑76 (MQ=255)
                                                                   tcatcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt          >  1:965787/1‑76 (MQ=255)
                                                                      tcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc       >  1:973793/1‑76 (MQ=255)
                                                                       caAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc      <  2:884379/76‑1 (MQ=255)
                                                                           aCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCg  >  2:713572/1‑76 (MQ=255)
                                                                           |                                                                           
CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCG  >  REL606/1688915‑1689065

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