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

CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCC  >  REL606/1688915‑1689061
                                                                           |                                                                       
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                         <  2:1491459/76‑1 (MQ=255)
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                         >  1:374904/1‑76 (MQ=255)
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                         <  2:56031/76‑1 (MQ=255)
  gggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc                                                                       >  1:380766/1‑76 (MQ=255)
   ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                      <  1:183378/76‑1 (MQ=255)
    gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                     <  1:582224/76‑1 (MQ=255)
    gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                     <  2:482835/76‑1 (MQ=255)
      gTTGACGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGcc                                                                   <  2:985706/76‑1 (MQ=255)
        tGAAGGGGTGAAACCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                 <  1:1592507/76‑1 (MQ=255)
        tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                 <  2:306470/76‑1 (MQ=255)
        tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                 <  2:984153/76‑1 (MQ=255)
        tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                 >  1:1393814/1‑76 (MQ=255)
            ggggTAAAACCCTATATCAAAACGATTGGTCTTTAAAACAGCAAAGCAAAAAATATCATCAAAGCCTGCCGTATCt                                                             <  2:897669/76‑1 (MQ=255)
                tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                         >  2:1124167/1‑76 (MQ=255)
                tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                         >  2:1009094/1‑76 (MQ=255)
                  aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTgg                                                       <  2:1861082/76‑1 (MQ=255)
                  aaaaCCGATAACAAAAGGATTGGGCTTTATAACAGCAAAGAAGAAAATATCATCAAACCCTGCCGTATCTTGCTgg                                                       <  2:151242/76‑1 (MQ=255)
                   aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa                                                      >  2:499592/1‑76 (MQ=255)
                      ccTATATAAAAACGATTGGGTTTTAAAACAGCAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGagca                                                   <  1:1337204/76‑1 (MQ=255)
                       cTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcag                                                  >  2:1516770/1‑76 (MQ=255)
                        tataTCAAAAGGATGGGGTTTAAAAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc                                                 <  2:1338814/76‑1 (MQ=255)
                        tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                 >  2:1388683/1‑76 (MQ=255)
                         atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                >  2:962284/1‑76 (MQ=255)
                         atatAAAAACGATTGGGCTTTATAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca                                                <  1:145354/76‑1 (MQ=255)
                              aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATg                                           <  1:1604829/76‑1 (MQ=255)
                               aaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATgg                                          <  2:1451642/76‑1 (MQ=255)
                                aaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGc                                         >  2:1947601/1‑76 (MQ=255)
                                 aCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                        >  1:127527/1‑76 (MQ=255)
                                  cGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                       <  2:1869634/76‑1 (MQ=255)
                                   gATTGGGCTTTAAAACAGAAAAGCAGAAAAAATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                      <  2:1984845/76‑1 (MQ=255)
                                     ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                    >  2:1297774/1‑76 (MQ=255)
                                     ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                    >  1:970551/1‑76 (MQ=255)
                                      tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                   >  2:1890278/1‑76 (MQ=255)
                                          ttttAAAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                               <  1:720317/75‑1 (MQ=255)
                                          cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                               >  2:881643/1‑76 (MQ=255)
                                            ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                             >  2:1491361/1‑76 (MQ=255)
                                              ataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGa                           >  1:53930/1‑76 (MQ=255)
                                              ataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGa                           >  1:1611652/1‑76 (MQ=255)
                                                aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc                         >  1:1448097/1‑76 (MQ=255)
                                                    gaaaaGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgc                     <  1:1879373/74‑1 (MQ=255)
                                                            naaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg             >  1:1812088/2‑76 (MQ=255)
                                                            aaaaaTATAAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg             <  1:1016388/75‑1 (MQ=255)
                                                             aaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa            >  2:609760/1‑76 (MQ=255)
                                                              aaaTATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa           <  2:1759940/76‑1 (MQ=255)
                                                                 tatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc        >  1:246193/1‑76 (MQ=255)
                                                                   acatcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt      <  1:1329987/75‑1 (MQ=255)
                                                                      tcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc   >  1:1303077/1‑76 (MQ=255)
                                                                      tcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc   >  2:1725839/1‑76 (MQ=255)
                                                                       caAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc  <  2:120258/76‑1 (MQ=255)
                                                                           |                                                                       
CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCC  >  REL606/1688915‑1689061

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: ‑