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

Only 100 of 127 total aligned reads displayed.

GAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720393‑720539
                                                                        |                                                                              
gAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcc                                                                           >  1:5456477/1‑76 (MQ=255)
gAGTTTTATGATTACCCATCTGGTCTGGTTGCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcc                                                                           <  1:62057/76‑1 (MQ=255)
 ggTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                          <  2:6740361/75‑1 (MQ=255)
 aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                          <  1:4087897/76‑1 (MQ=255)
 aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct                                                                          <  1:3762803/76‑1 (MQ=255)
  gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg                                                                         >  2:4814164/1‑74 (MQ=255)
  gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgcggccgc                                                                         >  2:3881296/1‑74 (MQ=255)
  gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                         <  2:3906749/76‑1 (MQ=255)
     ttATAATAACCCATCTGGTCTGGTTCCGCCAGGTTTTACTTCTCCACAATAATCGCGCCCTGGctcccccctgccg                                                                     <  2:332883/76‑1 (MQ=255)
       aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccccc                                                                   >  1:1628506/1‑73 (MQ=255)
        tGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccgcca                                                                 >  1:1474421/1‑74 (MQ=255)
          aTTCCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCGTAATAATCCCGCACGGGCTCCCGCCTGCCGCAAtt                                                               <  2:5798714/76‑1 (MQ=255)
           ttACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCGCCCCAAttt                                                              >  1:5104896/1‑75 (MQ=255)
           ttACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAAATc                                                              >  1:6802849/1‑76 (MQ=255)
           tcccccATCTGGTCTGGTTCCGCCAGAATTTCCGTCTCCACGATAACCTCGCCCGGCCGCCCCCCTGCCGCAATTc                                                              <  2:6952095/73‑1 (MQ=255)
            tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCACTCg                                                             >  2:3969474/1‑76 (MQ=255)
             aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGt                                                            >  1:5465381/1‑76 (MQ=255)
             aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCCGCCGCAACTCGt                                                            >  1:1155687/1‑76 (MQ=255)
              cccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCACTCGTc                                                           >  2:6437505/1‑76 (MQ=255)
               ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGGCCGGAATTctcc                                                           >  1:6368457/1‑73 (MQ=255)
               ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCCCCGCAAATCGTCt                                                          >  2:2696312/1‑76 (MQ=255)
               ccATCTGTTTGGTTTCCGCGGGGTTTTCCGTCTGCACGATACTCTCGCCCTGGCCCCCCCCTGCCGCAATTCGTCt                                                          <  2:6917553/76‑1 (MQ=255)
                cATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCCCCGCCATTCGTccg                                                         >  1:58999/1‑74 (MQ=255)
                 ttctggttgggTTTCGCCGGGTTTTCGGTCTCCAGAATATCCTCGCCCGGCCTCCCCCCTGCCGCAATTCGTCTGc                                                        <  2:185913/75‑1 (MQ=255)
                 aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCCCAATTCGTCTGc                                                        >  1:1967393/1‑76 (MQ=255)
                  tctggtctggTTCCGCCGGTATTGACGTCTGCACAATAGTCTCGCGCGGGCCCCCGCCTGCCGCAATTCGTCTGCa                                                       <  2:6973233/76‑1 (MQ=255)
                   ctggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCCAATCGTCTGCAc                                                      >  1:6522306/1‑76 (MQ=255)
                   ctggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCAc                                                      >  1:3683079/1‑76 (MQ=255)
                    tgttctggtTTCCCCAGGTTTTACGTCTGCACGATAATCTCGCACGGGCGGCCCCCTGCCGCAATTCGTCTGCAcg                                                     <  1:1648199/73‑1 (MQ=255)
                     ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCAcgc                                                    >  1:7353727/1‑76 (MQ=255)
                      gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCGCCCCACATCGTCTTCAcccc                                                   >  2:6481680/1‑73 (MQ=255)
                       tctgttTTCGCCAGGTTTTCCTTCGGCACGATAATCTCCCCCTGGCGCCCCCCTGCCGCAATTCGTCTGCACGCGt                                                  <  1:2234999/76‑1 (MQ=255)
                       tctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCGCCGCCATTCGCCTTCAcacgg                                                  >  2:7778770/1‑75 (MQ=255)
                       tctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAACTCGTCTGCCcgcct                                                  >  2:1488831/1‑74 (MQ=255)
                        ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGGCGCAATTAGACTCCAcacgcg                                                 >  1:4828091/1‑74 (MQ=255)
                        ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCAcgcgcg                                                 >  1:4426366/1‑74 (MQ=255)
                        ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAACTCGTCTCCCCGCGTg                                                 >  1:1269512/1‑76 (MQ=255)
                         tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCGGCCGCCATTCGTCTGCCCGCGTGc                                                >  2:5888788/1‑76 (MQ=255)
                          ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCGCCATTTCCTCGCACGCGTGcg                                               >  2:1017563/1‑75 (MQ=255)
                          ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCGCCCGCGATTCGTCTCCCCGCGGGcc                                               >  1:5603694/1‑75 (MQ=255)
                           gtTTCGCCAGGATTTACGTCGGCACGATATCCTCGCCCTGGCGCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTg                                              <  1:2522464/76‑1 (MQ=255)
                            tttCGCCGTGATTTCCGCCGGCACGATATCCCCCCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTgg                                             <  1:4821311/76‑1 (MQ=255)
                            tttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGGCGCCCTTTGTCTGCCAGCGGGCTgg                                             >  1:1034435/1‑76 (MQ=255)
                             ttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCCAATTGTCTGCACCCGGGGTggg                                            >  1:6140825/1‑75 (MQ=255)
                                cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTTGTCTGCACGCGggcgggcgg                                          >  1:6543213/1‑75 (MQ=255)
                                 gcgaGGTTTTGCGTCGCCAGGATAACCCCGCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGtt                                         <  2:6317960/73‑1 (MQ=255)
                                 gCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGGCCCCATTTGTCTGGCAGCGGGCTGGCggt                                         >  1:1414290/1‑74 (MQ=255)
                                 ccaaGGTTTTCCGTCTGCAGGATACTCTCGCCCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGtt                                         <  2:3157627/73‑1 (MQ=255)
                                  ccAGGTTTTCCGTCTCCACAATAATCCCGCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTta                                        <  1:1570561/76‑1 (MQ=255)
                                   cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGCCCGCGTGCTGGCGGtat                                       >  1:1066023/1‑76 (MQ=255)
                                   cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCCCCGCCATTCGTCTTCACGCGGGCTGGCGGtat                                       >  2:4317243/1‑76 (MQ=255)
                                   cAGGATTTACGTCTCCACGATAATCGCGCCCTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtat                                       <  1:7175735/76‑1 (MQ=255)
                                    aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCCGCCCCCATTTCTCTGCACCCGCGCGGGGGGtata                                      >  1:646221/1‑76 (MQ=255)
                                      gATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCCCCGCAAATTCTCTGTCCGCGTGGTGGCGTgatata                                    >  1:5351427/1‑76 (MQ=255)
                                       aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCAAATCATCTTCCCCCGCGCTGGCGttatttat                                   >  1:51721/1‑76 (MQ=255)
                                       aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCCATTGTCTGCACCCGCGCCGGCGTTatatat                                   >  2:3505309/1‑76 (MQ=255)
                                       aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCGTCTTCACCCGTGCTGGCGGtgtattt                                   >  2:6121015/1‑74 (MQ=255)
                                        tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCCCGCGTGCCGGGGTTATATATc                                  >  1:7040742/1‑76 (MQ=255)
                                        tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACGCGGGCTGGGGTTATATATc                                  >  1:7397900/1‑76 (MQ=255)
                                          tACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCCCCAATTGTCTGCCCGCGCGCCGGGGGTTTATATCGc                                >  2:2217843/1‑76 (MQ=255)
                                           aCGTCTGCAGAATACTCCCGCCGGGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt                               <  1:2165400/76‑1 (MQ=255)
                                            cGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCCATTCGCCTTCCCGCGCGGTGGCGGTATATATCGCTa                              >  1:6476993/1‑76 (MQ=255)
                                             gTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAAATTGTCTGCACGCGTGCTGGCGGTTTATATAGCtca                             >  1:1618321/1‑74 (MQ=255)
                                              tCTGCACTATACTCTCGCACGGGCTCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTaca                            <  2:3943489/76‑1 (MQ=255)
                                               cTGCAGAATAACCTCGCACTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTacac                           <  1:1933531/76‑1 (MQ=255)
                                                tGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACGCGTGCTGGCGTTTTATATCCCTAcacc                          >  1:6940487/1‑76 (MQ=255)
                                                 gCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGCACGCGGGCTGGGGTTTTATATCGCTAcaccc                         >  1:4675020/1‑75 (MQ=255)
                                                 gCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCACCCGCGCTGGCGGTATATATCGCTCcaccc                         >  2:5902886/1‑75 (MQ=255)
                                                 gCACGATAATCTCGCACTGGCTGGCGCCCCCCGCCAATCGTCTGCCACCGCGCTGGCGGTTTATATTACTAcacca                         >  2:7308785/1‑76 (MQ=255)
                                                   aggATAATCGCGCCCGGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg                       <  2:2291501/74‑1 (MQ=255)
                                                   aCGATATCCTCGCCCTGCCCGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg                       <  1:5384292/76‑1 (MQ=255)
                                                   aCGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTTCTCTGCACGCCTGGCGGCGGTTTTTATTGCTAcaccacc                       >  2:6529214/1‑75 (MQ=255)
                                                   aCGATAATCCCGCACGGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg                       <  2:2030947/76‑1 (MQ=255)
                                                    caaTAATCCCCCACTGCCGGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc                      <  2:6641347/74‑1 (MQ=255)
                                                    agATAATCTCCCCCGGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc                      <  1:573052/75‑1 (MQ=255)
                                                       tAATCTCGCACTGGCTCCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg                   <  2:5170626/76‑1 (MQ=255)
                                                       ataTCTCGCTCTGCCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg                   <  2:1348046/74‑1 (MQ=255)
                                                       aaaTCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg                   <  2:2346061/75‑1 (MQ=255)
                                                        aaTCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACGCGGGCTGGCGGTATATATCGCTACACCAccccccg                  >  2:2581693/1‑73 (MQ=255)
                                                        aaTCTCGCACTGGCTGGCGCCCGCCCCAAATCGTCTGGCCGCGTGCTGGGGTTATATATTGCTACACCACCCCacg                  >  2:7560845/1‑74 (MQ=255)
                                                           ctcGCAGTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTggg               <  2:6649312/76‑1 (MQ=255)
                                                           ctcGCACTGGCTGGCGCCTGCCCCCAATTCTTTTCCCCCGTGGTGGCGCTTTATATAGCTACACCACGCCACTggg               >  2:5465914/1‑76 (MQ=255)
                                                            tcGCCCGGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGc              <  2:4328116/76‑1 (MQ=255)
                                                             cGCCCTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCg             <  1:7825954/76‑1 (MQ=255)
                                                             cGCACTGGCTGGGGCCCGCCGCAATTCGTCTGCAAGCGGGCTGGCGGTATATATCGCTACACCCCCCCCCTGGGCg             >  2:3658054/1‑76 (MQ=255)
                                                             cGCACTGGCTGGCGCCCGCCCCCATTCGTCTGCACGCGTGCTGGGGGTATATATTGCTACACCACCCCCGTGGGCg             >  2:4454638/1‑76 (MQ=255)
                                                                  tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCa        >  1:2620211/1‑76 (MQ=255)
                                                                   ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt       >  2:2572206/1‑76 (MQ=255)
                                                                   ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt       >  2:3589421/1‑76 (MQ=255)
                                                                     ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa     >  1:4092071/1‑76 (MQ=255)
                                                                     cgcccgcctgccgcAATTCGTCGGCACGCGTGCTGGCGTTATATATCGCTACGCCACGCCAGTGGGCGACGCATaa     <  2:2067915/73‑1 (MQ=255)
                                                                      tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:6057206/1‑76 (MQ=255)
                                                                      tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:4626971/1‑76 (MQ=255)
                                                                      tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:1934252/1‑76 (MQ=255)
                                                                      tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:4885697/1‑76 (MQ=255)
                                                                      tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:2655632/1‑76 (MQ=255)
                                                                       gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  1:1070014/1‑76 (MQ=255)
                                                                       gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  1:7129302/1‑76 (MQ=255)
                                                                        ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  1:1466490/1‑76 (MQ=255)
                                                                        ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  2:5631628/1‑76 (MQ=255)
                                                                        |                                                                              
GAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720393‑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: