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

TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG  >  REL606/1688916‑1689062
                                                                          |                                                                        
tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                         <  1:1808567/76‑1 (MQ=255)
 gggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc                                                                        >  1:1476033/1‑76 (MQ=255)
  ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                       <  2:183103/76‑1 (MQ=255)
   gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                      <  2:1612280/76‑1 (MQ=255)
   gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                      <  2:572809/76‑1 (MQ=255)
      ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                   >  1:422459/1‑76 (MQ=255)
         aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                                >  2:1011929/1‑76 (MQ=255)
           ggggTAAAACCCTATAAAAAAACGATTGGGTTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCt                                                              <  2:507064/76‑1 (MQ=255)
                 aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTgg                                                        <  2:1640482/76‑1 (MQ=255)
                    cccTATATAAAAACGATGGGGCTTTATAACAGAAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                     <  1:1565340/75‑1 (MQ=255)
                    aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                     <  2:1645817/76‑1 (MQ=255)
                    aCCTATAACAAAACGATTGGGCTTTATAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                     <  1:116655/76‑1 (MQ=255)
                     ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                    >  1:742691/1‑76 (MQ=255)
                       ttnatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                  >  2:1326173/4‑76 (MQ=255)
                       tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                  >  2:1148591/1‑76 (MQ=255)
                       tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                  >  2:1715583/1‑76 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca                                                 <  2:1565608/76‑1 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca                                                 <  1:1540024/76‑1 (MQ=255)
                          aaaaaacCGTTTGGGCTTTAAAACAGCAAAGCAAAAAATATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATa                                               <  2:844247/73‑1 (MQ=255)
                            cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                             >  2:128767/1‑76 (MQ=255)
                            cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                             >  1:1759496/1‑76 (MQ=255)
                            cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                             >  1:1455423/1‑76 (MQ=255)
                            cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                             >  1:1399269/1‑76 (MQ=255)
                                 cGATTGGGCTTTATAACAGCAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                        <  1:1457927/76‑1 (MQ=255)
                                 cGATGGGGCTTTAAAACAGCAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                        <  1:1067579/76‑1 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                     >  2:1714248/1‑76 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                     >  2:1482589/1‑76 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                     >  1:74431/1‑76 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                     >  1:21968/1‑76 (MQ=255)
                                     tGGGTTTTATAACAGCAAAGCAAAAATTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                    <  2:874003/76‑1 (MQ=255)
                                        gCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                 <  1:1445298/76‑1 (MQ=255)
                                         cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                >  1:1578547/1‑76 (MQ=255)
                                         cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                >  1:127882/1‑76 (MQ=255)
                                         cTTTAAACCAGCAAAGCAAAAAATAACATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                <  1:880249/76‑1 (MQ=255)
                                          tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                               >  1:412170/1‑76 (MQ=255)
                                           ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                              >  1:1545409/1‑76 (MQ=255)
                                           ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                              >  2:209944/1‑76 (MQ=255)
                                              aaaCAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt                           <  2:1484867/75‑1 (MQ=255)
                                                aCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCg                         >  1:551763/1‑76 (MQ=255)
                                                     aaaGAAGAAAATACCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg                    <  2:663368/76‑1 (MQ=255)
                                                        gCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTgcgc                 <  2:641065/76‑1 (MQ=255)
                                                            aaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa             >  1:1881343/1‑76 (MQ=255)
                                                            aaaaTACCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa             <  1:208329/76‑1 (MQ=255)
                                                              aaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAg           >  2:265367/1‑76 (MQ=255)
                                                               atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGc          <  1:1483524/76‑1 (MQ=255)
                                                               atatCAACAACCCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGTTCGTGCTGCGCTTGAAGc          <  1:1042771/76‑1 (MQ=255)
                                                                 atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc        >  1:780912/1‑76 (MQ=255)
                                                                   catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg      <  2:1799657/76‑1 (MQ=255)
                                                                       aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg  >  1:717750/1‑76 (MQ=255)
                                                                          |                                                                        
TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG  >  REL606/1688916‑1689062

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