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

TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCG  >  REL606/1688916‑1689065
                                                                          |                                                                           
tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                            >  1:1008088/1‑76 (MQ=255)
tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                            >  1:1710829/1‑76 (MQ=255)
 gggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc                                                                           <  2:616226/76‑1 (MQ=255)
  ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                          <  1:1236257/76‑1 (MQ=255)
  ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                          <  1:29630/76‑1 (MQ=255)
   gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                         <  1:37960/76‑1 (MQ=255)
     ggtGAAGGTGTAAAAACCGATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                       <  2:529680/74‑1 (MQ=255)
     gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                       <  2:529683/76‑1 (MQ=255)
     gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                       <  2:273553/76‑1 (MQ=255)
        gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa                                                                    >  2:772891/1‑76 (MQ=255)
           ggggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt                                                                 >  1:1206152/1‑76 (MQ=255)
              gTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGc                                                              >  1:1662796/1‑76 (MQ=255)
                 aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg                                                           >  2:1808003/1‑76 (MQ=255)
                   aaCCTATATCAAAACGATGGGGCTTTATAACAGCAAAGCAGAAAATATCACCAAACCCTGCCGTATCTTGCTGGag                                                         <  2:1258625/76‑1 (MQ=255)
                    aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                        <  1:87818/76‑1 (MQ=255)
                     ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                       >  2:304384/1‑76 (MQ=255)
                     ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                       >  1:457860/1‑76 (MQ=255)
                      cTATATCAAAACGATGGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcag                                                      <  2:941169/76‑1 (MQ=255)
                       tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                     >  1:1489910/1‑76 (MQ=255)
                       tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                     >  2:1037381/1‑76 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                    >  2:1419663/1‑76 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                    >  2:1975994/1‑76 (MQ=255)
                        atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                    >  1:1929451/1‑76 (MQ=255)
                         tataAAAACGATGGGGTTTTAAAACAGCAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                   <  2:779162/76‑1 (MQ=255)
                         tatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAt                                                   >  2:1041587/1‑76 (MQ=255)
                           acAAAACGATGGGGCTTTAAAACAGAAAGCCAGAAAATATCACAAAACCCTGCCGTATCTTGCTGGAGCAGCATaa                                                 <  1:1721424/75‑1 (MQ=255)
                             aaaaCGATGGGGCTTGAAAACAGAAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATg                                               <  2:1206060/76‑1 (MQ=255)
                                aCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                            >  1:715771/1‑76 (MQ=255)
                                aCGATTGGGCTGTATAACAGCAAAGCAGAAAATATCACCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                            >  2:1974210/1‑76 (MQ=255)
                                  gATTGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                          <  1:814488/76‑1 (MQ=255)
                                    ttGGGCTTTATAACAGCAACGCAGAAAATATCATCAAAACCTGCCGAATCTTGCTGCAGCAGCATAATGGCGAggc                                        >  2:786150/1‑75 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                        >  2:786161/1‑76 (MQ=255)
                                    ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                        >  1:1522788/1‑76 (MQ=255)
                                     tGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                       <  2:1553025/76‑1 (MQ=255)
                                       ggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc                                     >  1:1848580/1‑76 (MQ=255)
                                          tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                                  >  2:243739/1‑76 (MQ=255)
                                           ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                                 >  1:659447/1‑76 (MQ=255)
                                           ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTCGAGCAGCATAATGGCGAGGTTCCGGaa                                 >  2:290305/1‑76 (MQ=255)
                                            nataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg                                >  1:1445315/2‑76 (MQ=255)
                                              taACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt                              >  1:730289/1‑76 (MQ=255)
                                               aaCAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc                             <  1:251664/76‑1 (MQ=255)
                                                  aaaaaaGAAAAAAATATCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtg                          <  1:826145/73‑1 (MQ=255)
                                                  aGCAAAGAAAAAAATATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtg                          <  2:910054/76‑1 (MQ=255)
                                                   gCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgc                         >  2:1144669/1‑76 (MQ=255)
                                                    cAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTCCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                        >  1:588853/1‑76 (MQ=255)
                                                     aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg                       >  2:1143779/1‑76 (MQ=255)
                                                     aaaGAAAAAAATATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg                       <  1:800937/76‑1 (MQ=255)
                                                            aaaaTATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa                <  2:53425/76‑1 (MQ=255)
                                                             aaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa               >  1:1196003/1‑76 (MQ=255)
                                                             aaaTATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa               <  1:723405/76‑1 (MQ=255)
                                                                  tcataaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt          <  1:1205544/76‑1 (MQ=255)
                                                                   catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg         >  1:861526/1‑76 (MQ=255)
                                                                   caccaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg         <  2:27170/73‑1 (MQ=255)
                                                                    atcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc        >  2:1131331/1‑76 (MQ=255)
                                                                          aCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCg  >  1:125275/1‑76 (MQ=255)
                                                                          |                                                                           
TGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCG  >  REL606/1688916‑1689065

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