Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4660GC55.0% ‑3.9 / 18.4 20A23P (GCC→CCC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base C (0/11);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 102 total aligned reads displayed.

AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGT  >  REL606/720391‑720541
                                                                             |                                                                            
aGGAGTTTTATGATTACCCATCTGTTGTGGTTTCGCCAGGTTTTACGTCTGCACAATAATCTTGCAGTGGctgccg                                                                              <  2:5661244/76‑1 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                              >  1:3419563/1‑76 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                              >  1:2285388/1‑76 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                              >  1:5067937/1‑76 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                              >  1:5728356/1‑76 (MQ=255)
aGGAGTTTTATGATGGCCCATCTGGTCTGGTTTCGCTAGGATTTACGTCGGCACGATAATCTCGCACTGGctgccg                                                                              <  2:5020856/76‑1 (MQ=255)
 ggATTTTTATGATTCCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgc                                                                             <  2:2185953/76‑1 (MQ=255)
  gAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcc                                                                            >  2:18767/1‑76 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgcggccc                                                                           >  2:6731898/1‑75 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                           <  1:5274649/76‑1 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                           <  1:1559702/76‑1 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                           <  2:195848/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                          <  1:3362271/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGGCTGCACGATAATCTCGCACTGGctgccgcctg                                                                          <  2:490993/76‑1 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgc                                                                        >  1:5743107/1‑76 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCCCGCCCGGCctgccgcctgc                                                                        <  2:6440301/76‑1 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgcc                                                                       >  1:4905349/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgcc                                                                       >  1:4271095/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgcc                                                                       >  1:6890716/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctgcc                                                                       <  2:3387839/76‑1 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcccctcccc                                                                      >  1:822972/1‑75 (MQ=255)
        tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcctccccc                                                                     >  2:2506122/1‑74 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcctgccccc                                                                    >  1:551216/1‑73 (MQ=255)
               aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCGGCCCAATTcct                                                              >  1:6805228/1‑74 (MQ=255)
               aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCACTTCGt                                                              >  1:1326927/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGCGCCTCCGCCCTTCGTCt                                                            >  2:6233444/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCGGCCGCGAATCGtat                                                            >  1:6840840/1‑74 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCATTTGTCt                                                            >  2:1242500/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCCACTCGTCt                                                            >  2:3203180/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCt                                                            >  2:5540163/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCGTCt                                                            >  1:4338118/1‑76 (MQ=255)
                  cATCTGGTCTGGTTTCGCCAGGATTTCCGTCGGCACGATAATCCCGCACTGCCTCCCCCCTGCCGCAATTCGTCTg                                                           <  2:5441934/76‑1 (MQ=255)
                  cATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCAATCGTCTg                                                           >  1:1089607/1‑76 (MQ=255)
                   aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGc                                                          >  1:3289937/1‑76 (MQ=255)
                    tctggtctggTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGcc                                                         >  1:4549788/1‑75 (MQ=255)
                     ctggtcggttTTCGCCAGGTTTTCCGTCTCCACGATAACCCCCCCCGGGCTGCCCCCTGCCGCAATTCGTCTGCAc                                                        <  1:6762227/76‑1 (MQ=255)
                      tggtttgttTGCGCCAGGATTTACGTCTGCACAATAATCTCGCACGGGCTCCCCCCTGCCGCAATTCGTCTGCAcg                                                       <  2:3084328/76‑1 (MQ=255)
                      tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGGCGCCAATCGTCTGCCcg                                                       >  2:3393254/1‑76 (MQ=255)
                      tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGcccc                                                       >  1:2389126/1‑73 (MQ=255)
                      ggttcttgtTTCGCCAGTATTTACTTCAGCACGATAATCCCGCCGTGGCTGCCGCCTGCCGCAATTCGTCTGCAcg                                                       <  1:6850326/73‑1 (MQ=255)
                       ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCAATTCGTCTGCAcac                                                      >  1:5569052/1‑74 (MQ=255)
                       ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCAATCGTCTGCAcgc                                                      >  1:3736649/1‑76 (MQ=255)
                        gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCAATTGTCCGCAcgcg                                                     >  1:330995/1‑76 (MQ=255)
                        ggctggtTCCGCCAGGTTTTGCGTCCGCACGATAACCCCGCCCTGCCCCCCCCCTGCCGCAATTCGTCTGCAcgcg                                                     <  1:4497098/74‑1 (MQ=255)
                          ggtgtTTCCCCAGTTTTTCCTTCTGCAGAATACTCCCGCACGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTg                                                   <  1:3589679/73‑1 (MQ=255)
                          ctgttTCCGCCGGGTTTTACGTCTGCACAATAATCTCCCACTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTg                                                   <  1:302751/76‑1 (MQ=255)
                          ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCCCGCGTg                                                   >  1:2287106/1‑76 (MQ=255)
                           tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCAATTCTCTGCACCCGTGc                                                  >  1:6138492/1‑76 (MQ=255)
                           tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCCCGCGTGc                                                  >  1:2835341/1‑76 (MQ=255)
                            ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCACTCGTCTGCACGCGTGCt                                                 >  1:5208994/1‑76 (MQ=255)
                               ttCGCCAAGGTTTTCCGTCTGCAGGATAATCTCGCACGGGCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTgg                                               <  1:6214465/76‑1 (MQ=255)
                               ttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTCCACGCGCGGTGGc                                              >  1:2598951/1‑76 (MQ=255)
                                tCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCGCCGCCAATCTTCTGCCCGCGGGCCGGCg                                             >  2:2169692/1‑76 (MQ=255)
                                 cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCGCCCGCCATTCGTCTGCCCCCGGGCTGGCGt                                            >  2:5323320/1‑76 (MQ=255)
                                    cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCAATTCGTCTGCACGCCTGGCGGCGTtat                                         >  2:1746544/1‑76 (MQ=255)
                                      aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGCCCGCGTGGTGGGGGtata                                        >  1:2900253/1‑76 (MQ=255)
                                      aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCCCGCGGGCTGGCGTtatt                                        >  1:3593775/1‑75 (MQ=255)
                                      aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAAATCGTCTGCACGCGTGCTGGCGTtata                                        >  1:6923089/1‑76 (MQ=255)
                                        tttttGCGCCTGCAGGATAATCCCGCCCGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatata                                      <  1:13578/74‑1 (MQ=255)
                                        gttttCCGTCTGCAAGATAACCCCGCCCTGGCGGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatata                                      <  1:6935646/74‑1 (MQ=255)
                                         aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCGTCTGCCCGCCGGCCGGCGTTatatat                                     >  1:6461401/1‑76 (MQ=255)
                                         aTTTACGTCTGCACGATAATCTCGCACTGGCCCCCCCCTGCCGCAATTCGTCTGCACACGTGCTGGCGTTatatat                                     <  2:2013650/76‑1 (MQ=255)
                                          tttCCGTCGCCACGATAACCGCGCAGGGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                    <  2:4209034/76‑1 (MQ=255)
                                          tttCCGGCTGCACGATATCCCCGCACGGCCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                    <  1:6559018/76‑1 (MQ=255)
                                          tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTTCCCGCGTGCTGGCGTTATATATc                                    >  1:1233691/1‑76 (MQ=255)
                                          tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATTGTCTCCCCCCCTGCGGGCGGTttatatt                                    >  2:1475303/1‑75 (MQ=255)
                                          tttACGTCTCCACGATACTCCCCCCCGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                    <  1:6724007/76‑1 (MQ=255)
                                           ttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATTGTCTGCACGCGTGCTGGGGGTatatattg                                   >  1:207186/1‑74 (MQ=255)
                                            tACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCGTCTGCCCGCGGGCTGGCGTTttatattcc                                  >  1:2507670/1‑73 (MQ=255)
                                             ccttCTGCTCGAAACTCTCGCGCGGCCCCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt                                 <  2:5516758/73‑1 (MQ=255)
                                               gTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGCACGCGGGCTGGGGTTATATATCGCTac                               >  1:2906760/1‑76 (MQ=255)
                                               gTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACGCGGGGTGGCGTTATATATCGCTac                               >  1:3380107/1‑76 (MQ=255)
                                                tCTGCAGGATAATCTCGCCCGGGCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTaca                              <  1:2123269/76‑1 (MQ=255)
                                                tCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCAATCGTCTGCACCCCCGCCGGCGGTATATATCTGTaca                              >  1:1887485/1‑76 (MQ=255)
                                                  tGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACGCGGGCTGGCGGTATATATCGCTAcacc                            >  1:3030808/1‑76 (MQ=255)
                                                   gCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTTGTCTGCACCCGTGCCGGGGTTTTATTTTGCTAcaccc                           >  2:4018512/1‑75 (MQ=255)
                                                   ccacaATAATCTCGCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacca                           <  1:2134600/75‑1 (MQ=255)
                                                    cACGATAATCTCGCACTGGCTGGCGCCTGCCGCAATTTGTCTGGACGCGCGCTGGGGCTTTTTATCGCTAcacccc                          >  1:6225561/1‑74 (MQ=255)
                                                    cACGATAATCCCGCACTGGCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcaccac                          <  1:1645536/76‑1 (MQ=255)
                                                    acagATACTCTCCCTCCGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcaccac                          <  1:157437/73‑1 (MQ=255)
                                                      cGATACTCTCGCACTGGCGCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc                        <  2:2881932/76‑1 (MQ=255)
                                                         tAATCTCGCACTGGCTGGCGGCTGCCCCCATTTGTCTTCACCCGTGGTGGCGTTTTATATCGCTACACCAcccccg                     >  2:5963388/1‑74 (MQ=255)
                                                          aaTCTCGCACTGGCTGGCGCCCGCCGCCATTTGTCTTCACCCGGGGTGGCGGTTTATTTTGCTTCACCCCCCCacg                    >  2:5676855/1‑74 (MQ=255)
                                                             ctcGCACTGGCTGGCGCCCGCCCCCAATCGTCTGCACGCGGGCTGGGGGTTTATATCGCTTCAACACACCCGCggg                 >  2:4203399/1‑76 (MQ=255)
                                                             cccGCACTGCCGCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTggg                 <  1:6305712/74‑1 (MQ=255)
                                                                   cTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGGGGGCGACGc           >  2:4929909/1‑76 (MQ=255)
                                                                     ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTCTATATCGCTACACCACGCCCGTGGGCGGAGCAt         >  2:3716533/1‑76 (MQ=255)
                                                                     ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt         >  2:4542223/1‑76 (MQ=255)
                                                                       ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa       >  1:1474768/1‑76 (MQ=255)
                                                                       ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa       >  2:16773/1‑76 (MQ=255)
                                                                        tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaaa      >  2:3007964/1‑75 (MQ=255)
                                                                        tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  2:4507328/1‑76 (MQ=255)
                                                                        tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  1:1758652/1‑76 (MQ=255)
                                                                        tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  2:5776236/1‑76 (MQ=255)
                                                                        tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  1:3891775/1‑76 (MQ=255)
                                                                        tgccccctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      <  2:3617871/76‑1 (MQ=255)
                                                                          gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa     >  2:4650880/1‑76 (MQ=255)
                                                                           ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt    >  2:3891556/1‑76 (MQ=255)
                                                                            cgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATg   >  2:159756/1‑76 (MQ=255)
                                                                             gcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGt  >  2:4471472/1‑76 (MQ=255)
                                                                             |                                                                            
AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGT  >  REL606/720391‑720541

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