Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,688,9900AC30.3% 178.7 / 14.6 76T92P (ACC→CCC) nthDNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)
Reads supporting (aligned to +/- strand):  ref base A (42/11);  new base C (0/23);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.68e-11
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
                                                                          |                                                                        
tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                          <  1:378664/75‑1 (MQ=255)
 gggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                         <  1:202991/75‑1 (MQ=255)
  ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc                                                                        <  2:634448/75‑1 (MQ=255)
  ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc                                                                        >  1:2014717/1‑75 (MQ=255)
   gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                       <  2:317631/75‑1 (MQ=255)
    cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                      <  1:1460544/75‑1 (MQ=255)
     gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                     >  1:319389/1‑75 (MQ=255)
     gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                     >  1:1715750/1‑75 (MQ=255)
      ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                    >  1:275653/1‑75 (MQ=255)
      ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                    >  1:1392497/1‑75 (MQ=255)
       tGAAGGGGTGAAGACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                   >  1:256891/1‑75 (MQ=255)
       tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                   <  1:1351429/75‑1 (MQ=255)
        gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                  >  1:5800/1‑75 (MQ=255)
        gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                  <  2:428359/75‑1 (MQ=255)
        gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                  >  1:345716/1‑75 (MQ=255)
        gAAGGGGTGAAAACCTATATCAAAACGACTGGGCTTTATAACAGCAAAGCAGAAAATATCATCCAAACCTGCCGt                                                                  >  1:5819/1‑75 (MQ=255)
         aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa                                                                 >  1:2020697/1‑75 (MQ=255)
           ggagTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATc                                                               >  2:370294/4‑75 (MQ=255)
            gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt                                                              <  2:498215/75‑1 (MQ=255)
            gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCt                                                              <  2:289651/75‑1 (MQ=255)
            gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGAAAAGAAGAAAATATCATCAAACCCTGCCGTATCt                                                              <  2:35394/75‑1 (MQ=255)
               tGAAACCCTATACCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAACCCCTGCCGTATCTTGc                                                           <  2:600813/75‑1 (MQ=255)
                gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCt                                                          <  2:1369528/75‑1 (MQ=255)
                gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                          <  2:1949566/75‑1 (MQ=255)
                 aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg                                                         >  1:771888/1‑75 (MQ=255)
                  aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg                                                        >  1:211880/1‑75 (MQ=255)
                  aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg                                                        >  2:93441/1‑75 (MQ=255)
                   aaCCTATATCAAAACGATTGGGTTTTATAACAGAAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGa                                                       <  2:2050974/75‑1 (MQ=255)
                   aaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa                                                       <  2:693607/75‑1 (MQ=255)
                   aaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa                                                       >  2:1333689/1‑75 (MQ=255)
                    aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGag                                                      >  2:254859/1‑75 (MQ=255)
                    aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGag                                                      >  2:1469962/1‑75 (MQ=255)
                    aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGag                                                      >  2:1215103/1‑75 (MQ=255)
                      cTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                    >  2:2038514/1‑75 (MQ=255)
                       tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcag                                                   >  2:230103/1‑75 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc                                                  <  2:577762/75‑1 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc                                                  <  2:1392497/75‑1 (MQ=255)
                          atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAt                                                >  1:1892063/1‑75 (MQ=255)
                          atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAt                                                >  1:61294/1‑75 (MQ=255)
                          atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                <  1:27211/75‑1 (MQ=255)
                          aacAAAACGATTGGGCTTTATAACAGCAAAGAAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                <  1:1566047/73‑1 (MQ=255)
                           tCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa                                               >  1:1606895/1‑75 (MQ=255)
                           tCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa                                               >  2:1251225/1‑75 (MQ=255)
                           acAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa                                               <  1:1391201/74‑1 (MQ=255)
                             aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAAt                                             <  2:526632/75‑1 (MQ=255)
                              aaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATg                                            >  2:2097146/1‑75 (MQ=255)
                              aaaCGATTGGGCTTTATAACAGAAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATg                                            <  1:559573/75‑1 (MQ=255)
                               aaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATgg                                           >  2:1667362/1‑75 (MQ=255)
                                aCGTTTGGGTTTTAAAACAGAAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGc                                          <  1:1726430/75‑1 (MQ=255)
                                aCGATTGGGTTTTAAAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGc                                          <  2:1696737/75‑1 (MQ=255)
                                  gtttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                        <  2:944711/73‑1 (MQ=255)
                                  gATGGGGTTTTAAAACAGAAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                        <  1:1299057/75‑1 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt                                                             <  2:568119/52‑1 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt                                                             >  1:568119/1‑52 (MQ=255)
                                      gggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                    >  1:644898/1‑75 (MQ=255)
                                       ggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGATc                                   >  1:1485000/1‑75 (MQ=255)
                                         cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                 >  2:1256075/1‑75 (MQ=255)
                                         cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                 >  1:943587/1‑75 (MQ=255)
                                         cTTTAAAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                 <  1:317675/75‑1 (MQ=255)
                                          tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg                                >  1:400338/1‑75 (MQ=255)
                                           ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                               >  2:1753669/1‑75 (MQ=255)
                                            tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                              >  1:1080603/1‑75 (MQ=255)
                                            tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                              >  1:1080585/1‑75 (MQ=255)
                                               aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt                           >  2:75082/1‑75 (MQ=255)
                                                     aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                     >  1:332939/1‑75 (MQ=255)
                                                     aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                     >  2:1672252/1‑75 (MQ=255)
                                                     aaaGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                     <  1:1465865/75‑1 (MQ=255)
                                                     aaaGCAAAAAATACCACAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                     <  2:633916/75‑1 (MQ=255)
                                                        gCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTgcg                  >  1:1274829/1‑75 (MQ=255)
                                                          aGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCt                >  1:827131/1‑75 (MQ=255)
                                                            aaaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg              <  1:2125423/75‑1 (MQ=255)
                                                            aaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg              >  1:455032/1‑75 (MQ=255)
                                                             aaaTATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa             <  2:1295177/75‑1 (MQ=255)
                                                             aaaTATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa             <  2:1522481/75‑1 (MQ=255)
                                                                 atcatcAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc         <  2:1984503/75‑1 (MQ=255)
                                                                  tcataaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc        <  1:523782/75‑1 (MQ=255)
                                                                   catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt       <  1:1820136/75‑1 (MQ=255)
                                                                   catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt       >  1:301633/1‑75 (MQ=255)
                                                                      caAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc    <  2:1099181/75‑1 (MQ=255)
                                                                        aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg  >  2:237425/1‑75 (MQ=255)
                                                                          |                                                                        
TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG  >  REL606/1688916‑1689062

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