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

GGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCC  >  REL606/1688918‑1689061
                                                                        |                                                                       
ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                      <  1:170914/76‑1 (MQ=255)
 gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                     >  1:1014960/1‑76 (MQ=255)
 gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                     <  1:1123986/76‑1 (MQ=255)
  cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                    >  2:254512/1‑76 (MQ=255)
    ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg                                                                  <  1:125669/76‑1 (MQ=255)
      gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa                                                                <  2:139956/76‑1 (MQ=255)
       aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                               <  1:1271203/76‑1 (MQ=255)
       aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                               <  2:471127/76‑1 (MQ=255)
          gggTGAAAACCTATATAAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCtt                                                            <  2:801758/76‑1 (MQ=255)
           ggTGAAACCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTg                                                           <  2:419119/76‑1 (MQ=255)
           ggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTg                                                           <  1:57192/76‑1 (MQ=255)
           ggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTg                                                           <  1:1316834/76‑1 (MQ=255)
             tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt                                                         >  1:1332797/1‑76 (MQ=255)
              gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg                                                        >  2:528004/1‑76 (MQ=255)
               aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg                                                       >  1:979238/1‑76 (MQ=255)
                aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa                                                      >  1:436879/1‑76 (MQ=255)
                  acatataTCAAAACGATTGGGCTTTATAACAGCAAAACAAAAAATATCATCAAACCCTGCCGTATCTTTCTGGagc                                                    <  2:187238/73‑1 (MQ=255)
                  aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                    <  1:1319880/76‑1 (MQ=255)
                   ccTATATCAAAACGATTGGGTTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca                                                   <  2:574718/76‑1 (MQ=255)
                   ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca                                                   >  1:329675/1‑76 (MQ=255)
                     tatataAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc                                                 <  2:1145079/76‑1 (MQ=255)
                     tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc                                                 <  1:659011/76‑1 (MQ=255)
                      atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                >  1:956355/1‑76 (MQ=255)
                      atatAAAAACGATTGGGCTTTATAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca                                                <  1:444989/76‑1 (MQ=255)
                       tatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                               <  1:528358/76‑1 (MQ=255)
                        atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa                                              >  1:116193/1‑76 (MQ=255)
                           aaaaCGATGGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATg                                           <  1:696943/76‑1 (MQ=255)
                            aaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATgg                                          >  1:200615/1‑76 (MQ=255)
                               cGATTGGGTTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                       <  2:1148722/76‑1 (MQ=255)
                                gATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg                                      >  1:101970/1‑76 (MQ=255)
                                 aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg                                     >  2:786845/1‑76 (MQ=255)
                                  ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                    >  2:914858/1‑76 (MQ=255)
                                    gggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTc                                  >  2:1267156/1‑76 (MQ=255)
                                     ggCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc                                 <  2:65677/76‑1 (MQ=255)
                                     ggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc                                 >  2:73438/1‑76 (MQ=255)
                                      gCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                <  2:362196/76‑1 (MQ=255)
                                      gCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg                                <  2:519723/76‑1 (MQ=255)
                                         ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa                             >  2:530052/1‑76 (MQ=255)
                                          tttaACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAGTGGCGAGGTTCCGGAAg                            >  1:469792/3‑76 (MQ=255)
                                          tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg                            >  1:76653/1‑76 (MQ=255)
                                            aaaCAGAAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt                          <  1:985565/75‑1 (MQ=255)
                                             aaCAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc                         <  1:1051021/76‑1 (MQ=255)
                                                 gCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgc                     <  1:361696/76‑1 (MQ=255)
                                                  cAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgat                    >  2:1204221/1‑74 (MQ=255)
                                                    aaGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctgc                  <  2:341637/76‑1 (MQ=255)
                                                           aaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa           <  1:1129271/76‑1 (MQ=255)
                                                           aaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa           <  2:4615/76‑1 (MQ=255)
                                                            aaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAg          >  1:1052881/1‑76 (MQ=255)
                                                             atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGc         >  1:1311258/1‑76 (MQ=255)
                                                             atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGc         >  2:356208/1‑76 (MQ=255)
                                                                 catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg     <  1:650321/76‑1 (MQ=255)
                                                                   tcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc   <  2:639843/76‑1 (MQ=255)
                                                                    caAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc  >  2:330077/1‑76 (MQ=255)
                                                                        |                                                                       
GGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCC  >  REL606/1688918‑1689061

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