Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4640CG50.0% 0.2 / 22.1 20A22G (GCC→GGC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base G (10/0);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
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.

Only 100 of 106 total aligned reads displayed.

AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720391‑720539
                                                                         |                                                                             
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                             >  2:8361971/1‑76 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                             >  2:8320164/1‑76 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg                                                                             <  2:6768984/76‑1 (MQ=255)
  gAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcc                                                                           <  2:2673380/76‑1 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                         <  2:4372579/76‑1 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                         <  1:2689819/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                        <  1:4442986/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                        <  1:2297103/76‑1 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgc                                                                      >  1:1054547/1‑76 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgc                                                                      >  2:1067235/1‑76 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgc                                                                      >  1:181170/1‑76 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctgc                                                                      <  1:5089141/76‑1 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgcc                                                                     >  1:60837/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctgcc                                                                     >  1:2825572/1‑76 (MQ=255)
       ttATGATTCCCCATCTGGTCTGGTTCCGCCAGAATTTACGTCGGCAGGATAATCTCGCACTGGctgccccctgccg                                                                    <  2:6710850/76‑1 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgccc                                                                    >  2:3581869/1‑75 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggccgccc                                                                    >  1:4978749/1‑75 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccc                                                                    >  1:1965067/1‑75 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccc                                                                    >  1:3709454/1‑75 (MQ=255)
        tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgcccc                                                                   >  2:1950501/1‑74 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggccgccgcc                                                                  >  2:1795141/1‑75 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgccggccgcc                                                                  >  1:3870114/1‑75 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccccc                                                                  >  1:5760397/1‑73 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCGCCGCa                                                                  >  1:8332936/1‑76 (MQ=255)
          tGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccgcca                                                                 >  1:959287/1‑74 (MQ=255)
          tGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGACGCac                                                                 >  1:3460797/1‑75 (MQ=255)
           gATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccgccac                                                                >  1:625903/1‑73 (MQ=255)
            aTTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCGGCCGCCAtt                                                               >  1:5664496/1‑76 (MQ=255)
            aTTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAtt                                                               >  1:5495811/1‑76 (MQ=255)
             ttACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTc                                                              >  1:4255015/1‑76 (MQ=255)
              tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCCGCAAtcg                                                              >  2:6204402/1‑74 (MQ=255)
              tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCg                                                             >  1:457663/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCt                                                          >  1:3953316/1‑76 (MQ=255)
                   aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGGCGCCAATCGTCttc                                                        >  2:7039889/1‑74 (MQ=255)
                     ctggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCGTCTGCAc                                                      >  1:321780/1‑76 (MQ=255)
                      tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCCCCCATCGTCTGCAcc                                                     >  2:3765643/1‑75 (MQ=255)
                      tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCCATCCTCTGCAcg                                                     >  1:6206513/1‑76 (MQ=255)
                        gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCCCTCGTCTTCAcgcc                                                   >  2:8350115/1‑75 (MQ=255)
                            ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCCCCGCCCATCGTCTTCACCCGTggt                                               >  1:3988015/1‑74 (MQ=255)
                            ggtTCCGCCAGGTTTTCCTTCAGCAGAATAATCCCGCACTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCt                                               <  2:7610845/76‑1 (MQ=255)
                             gtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCAATTCGTCTGCACGCGTggcg                                              >  1:8090466/1‑73 (MQ=255)
                             gtTTCGCAGGGATTTGCGTCTGCAGAATATTCTCGCTCTGCCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTg                                              <  2:6002955/76‑1 (MQ=255)
                              tttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCCATCGTCTTCTCGCGCGGTgg                                             >  2:7644549/1‑76 (MQ=255)
                              ttgcGCCAGTATTTAGGTCTGCAGAATACCCCCGCCTTGCCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTgg                                             <  2:2244740/73‑1 (MQ=255)
                              ttgcGCCAGGTTTTACGTCGGCGCGATAATCCCGCCCGGCCTCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTgg                                             <  1:3934289/73‑1 (MQ=255)
                                 cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCACTCCTCTGCCCGCGGGCTggggg                                          >  2:7150961/1‑73 (MQ=255)
                                 cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTTCCCCCCGGCTGGCgg                                          >  1:6555278/1‑75 (MQ=255)
                                 cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGt                                          >  1:1798654/1‑76 (MQ=255)
                                   ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCGTCTTCACGCGGGGTGGCggtt                                        >  1:813186/1‑73 (MQ=255)
                                    cAGGATTTCCTTCTGCACGATAATCCCGCCCGGGCGGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtat                                       <  1:910335/76‑1 (MQ=255)
                                    cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACGCGCGCTGGCGCtat                                       >  1:1213439/1‑76 (MQ=255)
                                     aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGGAATTCGTCTGCACGCGTGCTGGCGGtata                                      >  1:1966180/1‑76 (MQ=255)
                                      gtttttACGTCTGCAAGATACTCTCGCCGTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatat                                     <  2:271451/73‑1 (MQ=255)
                                      ggttttCCGTCTGCAGGATAATCCCGCACTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatat                                     <  2:4659464/73‑1 (MQ=255)
                                       tttttCCGTCTGCACAATAATCTCGCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatata                                    <  1:6713449/74‑1 (MQ=255)
                                       taTTTCCGTCTGCACAATATCCCCCCCCTGCCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatata                                    <  1:4621101/75‑1 (MQ=255)
                                        ttttGCGCCTCCACAATAATCTCGCCCGGGCTCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTatatat                                   <  1:4066959/75‑1 (MQ=255)
                                         tttCCGTCTGCCAGATACTCCCGCCCTGCCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                  <  2:6478063/76‑1 (MQ=255)
                                         tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGGCAATCCTCTGCACGCGCGCTGGCGTTatatata                                  >  1:1038078/1‑75 (MQ=255)
                                          ttACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCAATCGTCTTCCCGCGCGCTGGGGGTatttatat                                 >  2:2650637/1‑74 (MQ=255)
                                           tctgTCTGAACAATAATCGCGCCCTGCCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGc                                <  1:6604242/73‑1 (MQ=255)
                                            tcGTCCGCACGATAATCTCGCCGGGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt                               <  1:5633250/75‑1 (MQ=255)
                                            aCGTCTGCACGATAATCTCGCACTGGCTGGGGCCGCGCCCCAATTTTCTCCACCCGCGCGGGCGTTGTATATTGCt                               >  2:3660422/1‑76 (MQ=255)
                                             cGTCTGCCCGATACTCCCGCCCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTa                              <  2:3636592/76‑1 (MQ=255)
                                               tCTGCACGATAATCTCGCACTGGCTGCCGCCCGCCGCCATTCGTCTGCACCCCCGCTGGGGGTTTATATCGCTaca                            >  1:5586474/1‑76 (MQ=255)
                                                cTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACGCGTGCTGGCGGTATATATCGCTacac                           >  1:1497442/1‑76 (MQ=255)
                                                cTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCAATCGTCTGCACCCGGGCTGGCGGTTTATATCGCTacac                           >  2:7923937/1‑76 (MQ=255)
                                                cTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCCATTCGTCTGCACCCGTGGTGGGGGTATATATCGCTacac                           >  1:3783361/1‑76 (MQ=255)
                                                cTGCACGATAATCTCGCACTGGCTGCCGCCTCCCGCCATTCCCCTGCCCCCGTGCTGGGGTTTTATATTGCCacac                           >  1:4522078/1‑76 (MQ=255)
                                                 tccACGATACCCGCGCTCTCGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                          <  1:4419467/74‑1 (MQ=255)
                                                 tGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCCTTCGTCTGCACGCGGGCTGGCGGTATATATTTCTAcacc                          >  2:4318216/1‑76 (MQ=255)
                                                 ggCACGATACCCCCGCCCTGCCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                          <  2:6340055/75‑1 (MQ=255)
                                                 gccacaATAATCCCGCACTGCCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                          <  1:6899750/74‑1 (MQ=255)
                                                 gccACGATATCCTCGCCGGGCCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                          <  2:791374/74‑1 (MQ=255)
                                                  gCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacca                         >  1:1434100/1‑76 (MQ=255)
                                                  gCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTTCCCCCGTGCTGGCGGTTTATTTCGCTAcacca                         >  1:1199857/1‑76 (MQ=255)
                                                   cACTATAATCCCGCCCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcaccac                        <  1:3829329/76‑1 (MQ=255)
                                                    aCGATAATCTCGCCCTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg                       <  2:2567412/76‑1 (MQ=255)
                                                    aCGATAATCTCGCACTGGCTGGCGCCCGCCCCCAATCGTCTTCTCCCCCGCTGGCGGTTTATATTTCTTcaccacc                       >  2:3348617/1‑75 (MQ=255)
                                                     caaTAGTCTCGCCGTGGCGGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc                      <  2:4929666/74‑1 (MQ=255)
                                                     cGATAATCTCGCACTGGCTGGCGCCTGCCGCAATTCGTCTGCACGCGTGCCGGCGTTATATATTGCTAcaccaacc                      >  1:707895/1‑73 (MQ=255)
                                                      gATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCGTCTGCCCCCGCGCTTGGGGTTTATATCGCTACACCAcccc                     >  2:4372231/1‑76 (MQ=255)
                                                       aTAATCTCGCACTGGCTGGCCCCCGCCCCCATTTCTCTGCACGCCGGCCGGGGGTATATATTGCTTCACCACAccc                    >  2:6661725/1‑75 (MQ=255)
                                                        tAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTTCACCCGGGCTGGCGGTTTATATCGCTTCACCCCCCCAg                   >  2:2871819/1‑76 (MQ=255)
                                                         aaTCTCGCACTGGCTGGCGCCCGCCCCCATTTGTCTGCACCCGTGCTTGGGGTTTATATCGCTACACCCCCCCaca                  >  2:8152880/1‑74 (MQ=255)
                                                          aTCTCGCCCGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTg                 <  2:8320882/76‑1 (MQ=255)
                                                           ccccgcgCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTgg                <  2:5227659/73‑1 (MQ=255)
                                                                  cTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCCGTGGGCGACGc         >  1:3169068/1‑76 (MQ=255)
                                                                  cTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc         >  2:110568/1‑76 (MQ=255)
                                                                  cTGGCTGACCCCTGCCGCAATTCGTCAGCACGCCTGCTGGCGTTATATATCGCTGAAACACGCAAGCGGGCGACGc         >  2:2030324/1‑76 (MQ=255)
                                                                   tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACCCCACGCCAGGGGGCGACGCa        >  2:2857306/1‑76 (MQ=255)
                                                                    ggCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt       <  2:6202407/76‑1 (MQ=255)
                                                                    ggCCGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt       <  2:4123238/76‑1 (MQ=255)
                                                                     gCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATa      >  1:1632217/1‑76 (MQ=255)
                                                                       tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:6025175/1‑76 (MQ=255)
                                                                       tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:825339/1‑76 (MQ=255)
                                                                       ccccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    <  1:6285822/74‑1 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  1:5817405/1‑76 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  1:3725496/1‑76 (MQ=255)
                                                                         ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  1:3662379/1‑76 (MQ=255)
                                                                         |                                                                             
AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720391‑720539

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