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

CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG  >  REL606/1688915‑1689062
                                                                           |                                                                        
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                          >  1:1701984/1‑76 (MQ=255)
cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa                                                                          <  1:845294/76‑1 (MQ=255)
 tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATCGGGCTTTATAATAGCANAGCAGAAAATATCATCAAAAc                                                                         >  2:304653/1‑76 (MQ=255)
   ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                       <  1:1556609/76‑1 (MQ=255)
       ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                   >  1:655372/1‑76 (MQ=255)
        tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                  >  2:803236/1‑76 (MQ=255)
        tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                  <  2:515841/76‑1 (MQ=255)
          aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                                >  2:262369/1‑76 (MQ=255)
             gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt                                                             >  1:564930/1‑76 (MQ=255)
              ggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTg                                                            >  1:1422683/1‑76 (MQ=255)
                tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                          >  2:94357/1‑76 (MQ=255)
                tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                          >  2:1668898/1‑76 (MQ=255)
                tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                          >  2:1272333/1‑76 (MQ=255)
                 gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg                                                         >  1:610607/1‑76 (MQ=255)
                 gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg                                                         >  2:887505/1‑76 (MQ=255)
                  aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg                                                        >  2:81940/1‑76 (MQ=255)
                  aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg                                                        <  2:117225/76‑1 (MQ=255)
                     aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                     <  2:148960/76‑1 (MQ=255)
                       cTATATCAAAACTATTGGGCTTTATGACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcag                                                   <  2:426889/76‑1 (MQ=255)
                       cTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcag                                                   >  1:268706/1‑76 (MQ=255)
                        tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                  >  1:387940/1‑76 (MQ=255)
                        tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc                                                  >  1:728818/1‑76 (MQ=255)
                          tataAAAACGTTTGGGTTTTAAAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                <  1:81629/76‑1 (MQ=255)
                          tatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                <  1:1734614/76‑1 (MQ=255)
                           atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATa                                               <  2:76669/76‑1 (MQ=255)
                             cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                             >  1:389832/1‑76 (MQ=255)
                             cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                             >  2:37826/1‑76 (MQ=255)
                              aaaaCGATTGGGTTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAAta                                            <  2:796091/76‑2 (MQ=255)
                              aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATg                                            >  2:367736/1‑76 (MQ=255)
                              aaaaCGATGGGGTTTAATAACAACAAACCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATg                                            <  2:892586/76‑1 (MQ=255)
                               aaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATgg                                           >  2:1078059/1‑76 (MQ=255)
                                         gCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGTTTCCg                                 <  2:1670079/76‑1 (MQ=255)
                                         gCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                 >  1:630638/1‑76 (MQ=255)
                                           tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                               >  1:230459/1‑76 (MQ=255)
                                            ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                              >  2:411199/1‑76 (MQ=255)
                                                aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc                          >  2:482417/1‑76 (MQ=255)
                                                 aCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCg                         >  1:1833708/1‑76 (MQ=255)
                                                  cAGCAAAGCAAAAAATATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGt                        <  2:831679/76‑1 (MQ=255)
                                                      aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg                    >  2:1130910/1‑76 (MQ=255)
                                                      aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg                    >  2:495537/1‑76 (MQ=255)
                                                          aaGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCt                <  2:32853/75‑1 (MQ=255)
                                                           aaaaaaTATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCtt               <  2:1475682/74‑1 (MQ=255)
                                                             aaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa             >  1:1651856/1‑76 (MQ=255)
                                                             aaaaTATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa             <  2:802288/76‑1 (MQ=255)
                                                                 tatCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc         <  1:359263/76‑1 (MQ=255)
                                                                   tcatcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt       >  1:296698/1‑76 (MQ=255)
                                                                    catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg      <  2:711705/76‑1 (MQ=255)
                                                                    catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg      <  2:268706/76‑1 (MQ=255)
                                                                     atcaAACCCTGCCGTATCTTGCTGGAGCAACATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc     <  2:309915/76‑1 (MQ=255)
                                                                        aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg  <  1:34454/76‑1 (MQ=255)
                                                                        aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg  >  2:360719/1‑76 (MQ=255)
                                                                        aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg  >  1:706000/1‑76 (MQ=255)
                                                                           |                                                                        
CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG  >  REL606/1688915‑1689062

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