Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4640CG51.8% 5.1 / 41.8 27A22G (GCC→GGC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (2/11);  new base G (14/0);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 5.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 101 total aligned reads displayed.

CAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720390‑720539
                                                                            |                                                                            
cAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgcc                                                                             <  1:6416049/76‑1 (MQ=255)
 aGGAGTTTTATGATTACCCATCTGTTGTGGTTTCGCCAGGTTTTACGTCTGCACAATAATCTTGCAGTGGctgccg                                                                            <  2:5661244/76‑1 (MQ=255)
 aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                            >  1:5728356/1‑76 (MQ=255)
 aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                            >  1:5067937/1‑76 (MQ=255)
 aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                            >  1:3419563/1‑76 (MQ=255)
 aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                            >  1:2285388/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)
                     tctggtctggTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCCAATCGTCTGcc                                                       >  1:6687528/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)
                                                                            |                                                                            
CAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720390‑720539

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: