Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4640CG50.0% 3.0 / 23.9 20A22G (GCC→GGC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (3/7);  new base G (10/0);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 3.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720389‑720539
                                                                            |                                                                            
tCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgc                                                                              >  1:2654794/1‑76 (MQ=255)
  aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg                                                                            <  2:917135/76‑1 (MQ=255)
    gAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGTATTTACGTCTGCACGATAATCTCGCACTGGctgccgcc                                                                          <  2:6696229/76‑1 (MQ=255)
     ggTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCTCTGGctgccgcct                                                                        <  2:5202924/75‑1 (MQ=255)
     aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgccc                                                                         >  1:5530544/1‑75 (MQ=255)
     aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgccc                                                                         >  1:6580604/1‑75 (MQ=255)
     aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgccc                                                                         >  1:4787710/1‑75 (MQ=255)
     aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                        <  2:2820701/76‑1 (MQ=255)
     aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                        <  1:2071132/76‑1 (MQ=255)
      gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccc                                                                        >  1:5251492/1‑74 (MQ=255)
      gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                       <  1:2485821/76‑1 (MQ=255)
      gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                       <  2:3725001/76‑1 (MQ=255)
      gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                       <  1:2871223/76‑1 (MQ=255)
       ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgccccc                                                                       >  1:3898656/1‑73 (MQ=255)
       ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgc                                                                      >  1:3894024/1‑76 (MQ=255)
        tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgcc                                                                     >  2:5446144/1‑76 (MQ=255)
         ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccc                                                                    >  1:3664136/1‑75 (MQ=255)
          tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgccgc                                                                   >  2:3124286/1‑76 (MQ=255)
           aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccggc                                                                  >  1:4235772/1‑74 (MQ=255)
           aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccgcc                                                                  >  2:2797818/1‑75 (MQ=255)
            tGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCCCCGCac                                                                 >  2:1849788/1‑75 (MQ=255)
             gATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccgccca                                                                >  1:4220675/1‑73 (MQ=255)
             gATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAt                                                                >  1:3632373/1‑76 (MQ=255)
               ttACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTc                                                              >  1:6078683/1‑76 (MQ=255)
                tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCTGCCGCCATTcc                                                             >  2:1705253/1‑75 (MQ=255)
                tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCg                                                             >  1:1862012/1‑76 (MQ=255)
                 aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCACTCGt                                                            >  1:5851033/1‑76 (MQ=255)
                 aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCCCCGGCCACACTTTGt                                                            >  1:4811075/1‑76 (MQ=255)
                    cATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCGCCGCACATCCTCTg                                                         >  1:63727/1‑76 (MQ=255)
                     aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCAATCGTCtcc                                                        >  1:1423033/1‑74 (MQ=255)
                      tctgttctgtTTTCGCGAGGTTTTGCGTCGCCACGATACTCTCCCACGGGCTCCCCCCTGCCGCAATTCGTCTGCa                                                       <  2:927760/76‑1 (MQ=255)
                      tctggtctggTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCCATCGTCTggc                                                       >  1:2019691/1‑74 (MQ=255)
                       ctggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCTCCCGCCCTTCGTCTGCAc                                                      >  2:5270635/1‑76 (MQ=255)
                       ctggtctggtTCCGCCAGGTTTGCCGTCTGCACAATAACCCCCCCCTGGCTCCCCCCTGCCGCAATTCGTCTGCAc                                                      <  1:1964638/76‑1 (MQ=255)
                        tgttctgttGTCGCCAGGATTACCGTCTGCAAGATAATCTCCCCCGGGCTGCCCCCTGCCGCAATTCGTCTGCAcg                                                     <  2:3665459/73‑1 (MQ=255)
                         ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGCAcgc                                                    >  1:2371675/1‑76 (MQ=255)
                         ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCGCCCCCGCGCCACTGGTCTGCAcgc                                                    >  1:3027314/1‑76 (MQ=255)
                          gtctgttTCCGCCAGTTTTTACGTCCCCACGATAATCCCGCCCTGGCTGCCGCCTGCCGCAATTTGTCTGCAcgcg                                                   <  1:6495534/76‑1 (MQ=255)
                          gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCCTTCTTCCTCAcgcg                                                   >  2:2891536/1‑76 (MQ=255)
                            ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATTGTCTGCAcgcggg                                                 >  1:569529/1‑74 (MQ=255)
                             tggtTTGCGCGAGTATTTGCGTCGGCACAATAATCCCGCGCGGGCTCCCGCCTGCCGCAATTCGTCTGCACGCGTg                                                 <  2:592008/76‑1 (MQ=255)
                             tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCACTCGTCTTCCCACGTGc                                                >  2:6672304/1‑76 (MQ=255)
                              ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCATTAGTCTTCCCCCGTggg                                               >  2:1576451/1‑74 (MQ=255)
                               gtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCGTCTGCACGCgtggtg                                              >  1:6254595/1‑76 (MQ=255)
                                tttCGCCGGGTTTTCCGTCGGCACGATAATCCCGCCCTGGCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTgg                                             <  2:547809/76‑1 (MQ=255)
                                tttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCCTCTGCACGCCCGCTgg                                             >  1:420362/1‑76 (MQ=255)
                                  tCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCACTTGCCTTCCCCCGGGCTGGCg                                           >  2:1383929/1‑76 (MQ=255)
                                    cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAAATCGTCCGCACGCGCGCCGGCgg                                          >  1:2776453/1‑75 (MQ=255)
                                     gCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCGTCTGCACGCGGGCTGGCGtt                                         >  2:5183120/1‑76 (MQ=255)
                                      ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCAATTCGTCTGCACGCGTGCTGGCGttt                                        >  1:6416936/1‑75 (MQ=255)
                                       cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCGTCTGCCCGCGCGCTGGCGGtat                                       >  2:3214481/1‑76 (MQ=255)
                                        aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTTCACGCGTGCTGGCGGtata                                      >  1:1835660/1‑76 (MQ=255)
                                        aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACGCGGGCTGGCGGtatt                                      >  1:5771059/1‑75 (MQ=255)
                                          gATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACCCGTGGCGGCGGtttata                                    >  1:5549715/1‑76 (MQ=255)
                                            tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGGCCCCATTCGTCTGTACACGGGGTGGCGGTGTTTATc                                  >  1:5654214/1‑76 (MQ=255)
                                            tttACGTCGCCCCGATACTCCGGCGCTGGCCGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                  <  2:62696/76‑1 (MQ=255)
                                             ttACGTCTGCACGATAATCTCGCCCTGGCTGGGGCCTGCCCCAATTCATCTGCACGCGCGTTGGCGTTATATATCg                                 >  2:6793316/1‑76 (MQ=255)
                                             ttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAAATCGTCTGCCCCCGGGCGGGCGGTatatattc                                 >  1:3156607/1‑74 (MQ=255)
                                               gcctCTGCAAGATACCCTCCCGCGGCCCGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt                               <  1:6642888/73‑1 (MQ=255)
                                               ccGTCGGCAGGATACTCTCGCACGGGCTCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt                               <  2:4210526/75‑1 (MQ=255)
                                                cGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCAAATCGTCCGCACCCGTGCCGGCGGTATATATCGCTa                              >  1:1957616/1‑76 (MQ=255)
                                                 gTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGGGGTATATATTGCtcc                             >  1:1431049/1‑74 (MQ=255)
                                                    tccACGATACTCTCCCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                          <  1:5747765/74‑1 (MQ=255)
                                                     gCACGATAATCTCGCACTGGCTGCCGCCCGCCGCAAATCGTCTGCACGCGGGCTGGCGGTATATATTTCTAcacca                         >  1:5373357/1‑76 (MQ=255)
                                                        cGATAATCTCGCACTGGCTGGCGCCTGCCCCAATTCGTCTGCCCCCGGGCTGGGGTTTTATATCGCTAAACCACGc                      >  1:3305211/1‑76 (MQ=255)
                                                        cGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATTGCTACACCACGc                      >  1:5652461/1‑76 (MQ=255)
                                                          aTAGTCGCGCGCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCa                    <  1:3681137/76‑1 (MQ=255)
                                                          aTAATCTCGCACTGGCTGGGGCCCGCCGCCAATCGTCTGCACGCGCGCTGGCGTTATATATAGCTACACCAccccc                    >  2:6722013/1‑75 (MQ=255)
                                                           tactctcGCCCGGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg                   <  2:3338928/73‑1 (MQ=255)
                                                            aaTCTCGCACTGGCTGGCGCCGTCCCCCAATTGTCCGCCCGCCGGCCGGGGTTATATATCTCTAACCCACCCCCGt                  >  2:6878473/1‑76 (MQ=255)
                                                            aaTCTCGCACTGGCTGGCGCCCGCCCCCATTCCTCTGCCCGCCCGCCGGCCTTTTATATCGCTCCACCACGcccgc                  >  2:4425/1‑73 (MQ=255)
                                                                  ccACGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGa            <  2:4377546/75‑1 (MQ=255)
                                                                      tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGGCGCa        >  1:3745587/1‑76 (MQ=255)
                                                                      ctcctgccccctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCa        <  1:4425436/73‑1 (MQ=255)
                                                                       ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt       >  1:5323469/1‑76 (MQ=255)
                                                                       ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt       >  2:4023774/1‑76 (MQ=255)
                                                                         ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa     >  2:2706796/1‑76 (MQ=255)
                                                                         ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa     >  1:42203/1‑76 (MQ=255)
                                                                          tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:3773337/1‑76 (MQ=255)
                                                                          tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:5183215/1‑76 (MQ=255)
                                                                          tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:6378753/1‑76 (MQ=255)
                                                                          tgccacctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:4748944/1‑76 (MQ=255)
                                                                          cgccccctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    <  1:5912811/75‑1 (MQ=255)
                                                                          ccccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    <  2:4590405/74‑1 (MQ=255)
                                                                           gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGGCGCATAACa   >  1:496103/1‑76 (MQ=255)
                                                                           cccgcctgccgcAATTCGTCTGCACGCGTGCCGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   <  2:2873160/75‑1 (MQ=255)
                                                                            ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  2:4929728/1‑76 (MQ=255)
                                                                            |                                                                            
TCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720389‑720539

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: