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

TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG  >  REL606/1688916‑1689063
                                                                          |                                                                         
tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                          <  1:2056585/76‑1 (MQ=255)
  ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                        >  2:548896/1‑76 (MQ=255)
   gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                       <  1:163280/76‑1 (MQ=255)
    cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                      >  2:551223/1‑76 (MQ=255)
    cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                      >  2:1381669/1‑76 (MQ=255)
     gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                     >  1:1682507/1‑76 (MQ=255)
         aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                                 <  2:245216/76‑1 (MQ=255)
          aGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATc                                                                >  1:443150/1‑76 (MQ=255)
             ggTAAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTg                                                             <  2:1443352/76‑1 (MQ=255)
               taaaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCt                                                           <  2:854027/74‑1 (MQ=255)
               tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                           >  1:481260/1‑76 (MQ=255)
                gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg                                                          >  2:1455023/1‑76 (MQ=255)
                  aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATACCAACAAACCCTGCCGTATCTTGCTGGa                                                        <  2:91832/76‑1 (MQ=255)
                     ccTATATCAAAACTATTGGGCTTTAAAACAGCAAAGAAAAAAAAACCATCAAACCCTGCCGTATCTTGCTGGagca                                                     <  2:1609540/76‑1 (MQ=255)
                       tatataAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc                                                   <  1:2014622/76‑1 (MQ=255)
                       tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                   >  2:76306/1‑76 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca                                                  <  2:666569/76‑1 (MQ=255)
                        atatAAAAACGATTGGGCTTTATAACAGAAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAgcagca                                                  <  1:1270010/76‑1 (MQ=255)
                          atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATa                                                <  2:808515/76‑1 (MQ=255)
                          atCAAAACGATTGGGCTTAATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATa                                                <  2:1870691/76‑1 (MQ=255)
                            cAAAACGTTTGGGATTTATAACAGCAACGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                              <  2:1733001/76‑1 (MQ=255)
                            cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                              >  1:1799688/1‑76 (MQ=255)
                            aaaaaCGATTGGGTTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAAt                                              <  2:1842040/75‑1 (MQ=255)
                              aaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATgg                                            >  2:1360583/1‑76 (MQ=255)
                               aaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGc                                           <  1:2015453/76‑1 (MQ=255)
                                 cGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                         >  1:533981/1‑76 (MQ=255)
                                  gATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                        <  2:1348682/76‑1 (MQ=255)
                                   aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg                                       >  1:895346/1‑76 (MQ=255)
                                     tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                     >  1:280519/1‑76 (MQ=255)
                                       ggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc                                   >  2:280049/1‑76 (MQ=255)
                                        gCTTTAAAAAAACAAAACAAAAAATATAATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                  <  1:917482/76‑1 (MQ=255)
                                         cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                 >  2:1055024/1‑76 (MQ=255)
                                          tttATAACAGCAAAGCAGAAAATATCATCAAGACCTGCTGTATCTTGCTTGAGCAGCATAATAGCGCGGTTCCGGa                                >  2:83337/1‑76 (MQ=255)
                                          tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                                >  2:307594/1‑76 (MQ=255)
                                          tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                                >  2:1904019/1‑76 (MQ=255)
                                          tttATAACAGCAAAGAAAAAAATATCATAAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                                <  1:91629/76‑1 (MQ=255)
                                            tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg                              >  1:1130556/1‑76 (MQ=255)
                                             ataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGa                             >  2:1009421/1‑76 (MQ=255)
                                                    cAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                      >  2:1786587/1‑76 (MQ=255)
                                                     aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg                     >  2:374205/1‑76 (MQ=255)
                                                          aGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCtt                >  1:1613636/1‑76 (MQ=255)
                                                           gAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg               >  1:1314615/1‑76 (MQ=255)
                                                             aaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa             >  2:1727651/1‑76 (MQ=255)
                                                              aaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAg            <  1:2039969/76‑1 (MQ=255)
                                                                tatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc          >  1:1424421/1‑76 (MQ=255)
                                                                  tcatcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt        <  2:522717/76‑1 (MQ=255)
                                                                  tcatcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt        <  2:60421/76‑1 (MQ=255)
                                                                    ctcaATACCTGCCGTATATTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGCTCGTGCTGCGCTTGAAGCCCTGc      <  2:1934619/75‑1 (MQ=255)
                                                                    aacaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc      <  2:600394/74‑1 (MQ=255)
                                                                    aacaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc      <  2:1934627/74‑1 (MQ=255)
                                                                      caAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc    >  2:656369/1‑76 (MQ=255)
                                                                       aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg   <  2:520248/76‑1 (MQ=255)
                                                                       aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg   >  2:1884494/1‑76 (MQ=255)
                                                                        aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg  >  2:226224/1‑76 (MQ=255)
                                                                        aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgc  >  2:1296861/1‑75 (MQ=255)
                                                                          |                                                                         
TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG  >  REL606/1688916‑1689063

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