Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4660GC55.0% ‑4.0 / 22.7 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 = 8.00e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGT  >  REL606/720391‑720541
                                                                           |                                                                           
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg                                                                             >  1:6004999/1‑76 (MQ=255)
  gAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcc                                                                           >  1:3320771/1‑76 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg                                                                         <  2:303657/76‑1 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccgccc                                                                      >  2:5438815/1‑75 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggctcccc                                                                      >  2:4606201/1‑75 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccg                                                                      >  1:5191560/1‑76 (MQ=255)
       tgaTGATTCCCCACCTGTTCTGGTTTCGCGAGGATTTCCGTCTGCCCGATAATCTCGCACTGGctgccccctgccg                                                                      <  2:753464/74‑1 (MQ=255)
        tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgccgc                                                                     >  1:5981766/1‑76 (MQ=255)
        tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccgcccc                                                                     >  1:2417215/1‑74 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggctgccgcc                                                                    >  1:2597410/1‑75 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggccgcccgc                                                                    >  2:6270811/1‑73 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcctgccgcc                                                                    >  2:4502308/1‑75 (MQ=255)
            aTTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAtt                                                                 >  1:903980/1‑76 (MQ=255)
               aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCGCCGCCAATCGt                                                              >  2:3670554/1‑76 (MQ=255)
                 ccATCGGGTCGGGTTGCGCCAGGTTTTACGTCGGCACGATAATCTCGCACTGGCTCCCCCCTGCCGCAATTCGTCt                                                            <  1:4721238/76‑1 (MQ=255)
                  cATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCGCCGCAAATCCTccg                                                           >  2:189830/1‑74 (MQ=255)
                   aTCTGGTTTGGTCTCGCCGGGTTTTCTTCCTGCACGATAATCTGCCATTGTCTCCCCCCTGCCGCAATTCGTCTGc                                                          <  1:2959261/76‑1 (MQ=255)
                    tgtggtgtgttTCCACCGGGATTTGCTGCTGCACGATACCCCCGCCCGGGCTCCCCCCTGCCGCAATTCGTCTGCa                                                         <  2:5074188/74‑1 (MQ=255)
                    tctgttcttgTTCCGCCGGATTTTCCGTCTCGACGATATCCCCGCCCGGGCTCCCCCCTGCCGCAATTCGTCTGCa                                                         <  2:2457630/76‑1 (MQ=255)
                     ttggtgtgttTCCGCCAGGTTTTGCGTCGGCCGAATAACCCCGCCCTGGCGGCCCCCTGCCGCAATTCGTCTGCAc                                                        <  1:1640565/75‑1 (MQ=255)
                      tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCCCAATTCGTCTGCAcg                                                       >  1:3044656/1‑76 (MQ=255)
                       gttctgttTTCGCCGGGTTTTCCGTCTGCACGATACTCTCGCACTGGCTCCCCCCTGCCGCAATTCGTCTGCAcgc                                                      <  1:2471360/74‑1 (MQ=255)
                         tctggtTCCGCCAGTTTTTCCGTCTGCACGATAATCTCCCCCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGt                                                    <  2:961962/76‑1 (MQ=255)
                          ctgttTTCGCCAGGTTTTCCGTCTCCACGATAACCCCGCACTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTg                                                   <  1:5129772/76‑1 (MQ=255)
                          ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCCATTCGTCTGCACGCGTg                                                   >  1:2845556/1‑76 (MQ=255)
                          ctggtCTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCAcccgcg                                                   >  1:3959086/1‑74 (MQ=255)
                            ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACCCGTGCt                                                 >  1:6333848/1‑76 (MQ=255)
                            ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCACTCGTCTGCCCGCGCGCt                                                 >  1:3399671/1‑76 (MQ=255)
                            ggtTTCGCCAGGATTTACGTCTGCACGATAAGCTCGCACTGGCTGGGGGCCGCCCCCAATCGTCTCCACgcgcgcg                                                 >  2:4902792/1‑75 (MQ=255)
                             gtTCCGCCGGGTTTTCCGCCTGCCGAATATTCCCCCACGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTg                                                <  2:666423/76‑1 (MQ=255)
                               ttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATCGCCTCCACGCGTGCTGGc                                              >  1:5850544/1‑76 (MQ=255)
                                tCGCCAGTTTTTACGTCTGCACGATACTCTCGCCCTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCg                                             <  1:2145485/76‑1 (MQ=255)
                                  gCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAACTTGTCTTCACCCCGGCTGGCGtt                                           >  2:2186969/1‑76 (MQ=255)
                                  gCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCCCCCCCATTCCTCCGCCCGCCGGCTGGCggt                                           >  2:5574007/1‑74 (MQ=255)
                                     aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCCCCAATCGTCTTCCCGCGCGCCGGCGGtatt                                        >  2:4400043/1‑75 (MQ=255)
                                     aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCCCGCGTGCTGGCGTtata                                        >  1:4392892/1‑76 (MQ=255)
                                      ggATTTACGTCTGCACGATAATCTCGCACGGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatat                                       <  2:6285907/76‑1 (MQ=255)
                                        ttttACGTCTGCAGAATAATCCCGCGCTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTatatat                                     <  1:4967096/75‑1 (MQ=255)
                                        aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCAATTGTCTGCACGCGGGCTGGGGGTttatat                                     >  2:5601887/1‑76 (MQ=255)
                                         tttCCGCCTGCAGAATATTCCCGCACTGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                    <  1:1045422/76‑1 (MQ=255)
                                          ttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAATTCGTCTGCCCCCGTGCTGGGGTTatatattg                                   >  1:330181/1‑74 (MQ=255)
                                           tACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCGCCGCAATTCGTCTGCACCCGGGCCGGCGGTTTATATCGc                                  >  2:5151256/1‑76 (MQ=255)
                                            aCGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCCCGCGTGCTGGGGTTTTATATCGcc                                 >  1:3128322/1‑75 (MQ=255)
                                            aCGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCAAATTGTCTGCACGCGGGCTGGCGGTATATATCGCt                                 >  1:879291/1‑76 (MQ=255)
                                             cGTCTGCACGATAATCTCGCACTGGCTGGCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATTGCTa                                >  1:303963/1‑76 (MQ=255)
                                             cGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCCATTCGTCTGCACGCGGGCTGGCGTTATATATCGCTa                                >  1:4728041/1‑76 (MQ=255)
                                             cGTCTGCACGATAATCTCGCACTGGCTGCCGCCCGCCGCCAATCGTCTGCACCCCGGCTGGCGGTTTATATCTCTa                                >  1:678741/1‑76 (MQ=255)
                                             cGTCGGCACAATAACCTCGCCCTGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTa                                <  2:2844782/76‑1 (MQ=255)
                                               tctcCACGATAATCTCGCCCGGGCGCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTaca                              <  1:4055942/76‑1 (MQ=255)
                                               tCTGCACGATAATCTCGCACTGGCTGGCGCCCGCCGCCATTCGTCTGCACGCGGGCTGGCGTTTTATATCGCttcc                              >  1:5745822/1‑73 (MQ=255)
                                                cTGCACGATAATCTCGCACTGGCTGGCGCCTGCCCCAATTCCTCTTCACGCGCGCCGGCGGTTTATTTTGCTacac                             >  2:1499519/1‑76 (MQ=255)
                                                 tGCCCGATACCCCCGCCCTGGCTGCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                            <  2:4055693/76‑1 (MQ=255)
                                                 tGCACGATACTCGCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc                            <  2:2472516/76‑1 (MQ=255)
                                                  gCACGATAATCTCGCACTGGCTGGCGCCCCCCGCCACTTGTCTTCACGCGTGCTGGCGGTTTATATTTCTAcccca                           >  2:4671187/1‑76 (MQ=255)
                                                  gCACGATAATCTCGCACTGGCTGCCCCCCCCCCCCAATCGTCTTCACGCGCGCTGGCGGTATATATCCCcacaacc                           >  1:2416977/1‑73 (MQ=255)
                                                   caagATAACCTCGGGCGGGCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcaccac                          <  1:678900/73‑1 (MQ=255)
                                                    aCGATATTCTCGCCCGGGCCCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg                         <  1:5562703/76‑1 (MQ=255)
                                                     cGATAATCTCGCACTGGCTGGCGCCCGCCGCACTTCGGCTGCACGCGCGGTTGCGGTTTATATTCCTACCACACGc                        >  1:4246210/1‑76 (MQ=255)
                                                        tAATCTCGCACTGGCTGGCGCCCGCCGCCATTCCTCTGCACGCGGGCCGGCGGTATATATTGCTTCACCCCGcccc                     >  2:158118/1‑74 (MQ=255)
                                                             ccGCCCTGCCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGc                <  1:4935821/75‑1 (MQ=255)
                                                              cGCACTGGCTGGCGCCTGCCCCCCATCGTCTGCACCCGGGGTGGCGGTATATATCGCTACACCACACCCGCggggg               >  2:3186783/1‑74 (MQ=255)
                                                                 cggggCTCCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACg            <  2:2838510/73‑1 (MQ=255)
                                                                  cTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc           >  1:2328420/1‑76 (MQ=255)
                                                                   tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCATGCCAGTGGGCGACGCa          >  2:4234682/1‑76 (MQ=255)
                                                                   tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGGGACGCa          >  2:5543685/1‑76 (MQ=255)
                                                                   tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCa          >  2:4352319/1‑76 (MQ=255)
                                                                     gCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATa        >  1:611515/1‑76 (MQ=255)
                                                                      ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa       >  1:4505161/1‑76 (MQ=255)
                                                                      ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa       >  2:946454/1‑76 (MQ=255)
                                                                       tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  2:3669775/1‑76 (MQ=255)
                                                                       tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  1:5632657/1‑76 (MQ=255)
                                                                       tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc      >  2:3762404/1‑76 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa     >  2:4041438/1‑76 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa     >  2:5214040/1‑76 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa     >  2:2072114/1‑76 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa     >  1:5372097/1‑76 (MQ=255)
                                                                        gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa     >  2:905892/1‑76 (MQ=255)
                                                                         ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt    >  2:1867758/1‑76 (MQ=255)
                                                                           gcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGt  >  2:5588280/1‑76 (MQ=255)
                                                                           gcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGt  >  2:602259/1‑76 (MQ=255)
                                                                           |                                                                           
AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATGT  >  REL606/720391‑720541

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