Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4640CG55.2% ‑0.1 / 38.4 29A22G (GCC→GGC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (3/10);  new base G (16/0);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 1.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-04
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.

AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720391‑720539
                                                                           |                                                                           
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccg                                                                           <  1:4895371/76‑1 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccg                                                                           <  2:5873046/76‑1 (MQ=255)
aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccg                                                                           <  1:2405041/76‑1 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggcgccc                                                                        >  2:4423356/1‑75 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcct                                                                        <  2:5334557/76‑1 (MQ=255)
   aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGCATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcct                                                                        <  2:5231103/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcctg                                                                       <  2:2432451/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcctg                                                                       <  1:729977/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcctg                                                                       <  2:2377495/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcctg                                                                       <  2:3580059/76‑1 (MQ=255)
    gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑gccgcctg                                                                       <  2:4279758/76‑1 (MQ=255)
     ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCtggccgcctc                                                                       >  2:2444270/1‑75 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggggcctgcc                                                                     >  1:6657149/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgcc                                                                     >  1:1193378/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgcc                                                                     >  1:3738663/1‑76 (MQ=255)
      tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCtggccgccggc                                                                      >  1:4734734/1‑76 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggggcctgccc                                                                    >  1:1356495/1‑75 (MQ=255)
       ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgccg                                                                    >  1:8142506/1‑76 (MQ=255)
        tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggggcccgcccc                                                                   >  2:2995653/1‑74 (MQ=255)
        tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggcgcctgcaac                                                                   >  2:5995392/1‑73 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggggccggccccc                                                                  >  1:7094912/1‑73 (MQ=255)
         aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGct‑ggggcccgccgcc                                                                  >  1:4876582/1‑75 (MQ=255)
           gATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCAt                                                                >  1:2401111/1‑76 (MQ=255)
             ttACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCTGCCGCAAtac                                                              >  1:1929406/1‑74 (MQ=255)
             ttACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCCCACTTc                                                              >  2:7995760/1‑76 (MQ=255)
               aCCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGTTAATCTCGTACTGGCT‑GGGGGCTGCCCTCATTCGt                                                            >  2:1942137/1‑76 (MQ=255)
               aCCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGt                                                            >  1:8055241/1‑76 (MQ=255)
                 ccATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGGCGCACTTCGTcc                                                          >  1:4500865/1‑75 (MQ=255)
                   aTCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCGCTTCGTCTGc                                                        >  2:4033556/1‑76 (MQ=255)
                   aTCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCCCCCCAATTCGTCtcc                                                        >  2:3250394/1‑74 (MQ=255)
                    tctggtctggTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGGCCCGCCCCCAATCGTCTGcc                                                       >  1:7630587/1‑75 (MQ=255)
                    tctggtctggTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCCCCGCCACTCCTCTGcc                                                       >  2:3035338/1‑75 (MQ=255)
                     gtggtctggtTTCGCGAGGATTTA‑CGTCTGCACGATAATCTCGCAGTGGCT‑CCCCCCTGCCGCAATTCGTCTGCAc                                                      <  2:1987259/75‑1 (MQ=255)
                     cgggtctggtTCCCCCAGGATTTA‑CTTCTGCAGAATACTCCCGCCCGGGCT‑CCCCCCTGCCGCAATTCGTCTGCAc                                                      <  1:650880/74‑1 (MQ=255)
                      tggtctggtTTCGCGAGTTTTTC‑CGTCTGCACAATAATCTCCCCCTGGCT‑GCCGCCTGCCGCAATTCGTCTGCAcg                                                     <  1:7139138/76‑1 (MQ=255)
                      tggtctggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCTGCCGCCACTTGTCTGCAcc                                                     >  1:7309249/1‑75 (MQ=255)
                       ggtctggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGGGCCGCCCCCCTTCGTCTGTAcgc                                                    >  2:1002719/1‑76 (MQ=255)
                       ggtctggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCCCCCGCCGCAATTCCTCCGCAccc                                                    >  1:131263/1‑74 (MQ=255)
                         tttggtTTCCCCAGGTTTTA‑CGCCGCCACAATAATCCCGCCCTGGCC‑GCCGCCTGCCGCAATTCGTCTGCACGCGt                                                  <  1:2587443/74‑1 (MQ=255)
                         tctggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGGCCGCCCCCACTCCTCTTCAcgcgg                                                  >  2:18059/1‑75 (MQ=255)
                          ctggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCTCCCGCCATTCGTCTGCAcccggg                                                 >  1:4753110/1‑74 (MQ=255)
                          ctggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGTCTGCAcgcgcg                                                 >  1:5189976/1‑74 (MQ=255)
                           tggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGGCCCGCCGCCAATCGTGCGCCCCCGTGc                                                >  2:1379894/1‑76 (MQ=255)
                            ggtTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGGCCCGCCCCCCTTCCTCTGCCCCCGCGCt                                               >  2:109720/1‑76 (MQ=255)
                             ttttCGCCGGTTTTTC‑CGTCTGCACGATAATCTCCCCCGGGCT‑GCCCCCTGCCGCAATTCGTCTGCACGCGTGCTg                                              <  2:1083307/75‑1 (MQ=255)
                              tttCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCAATCGTCTGCACGCGGGCTgg                                             >  1:3763947/1‑76 (MQ=255)
                              tttCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCAATTCGTCTGCACGCGGGGTgg                                             >  1:5674589/1‑76 (MQ=255)
                               ttCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGGCCCGCCCCCACTCCTCTTCCCGCGtggtggg                                            >  2:2746005/1‑75 (MQ=255)
                                tCGCCAGTATTTA‑CGTCTGCCCGATAATCTCGCACTGGCT‑GCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCg                                           <  1:1999788/76‑1 (MQ=255)
                                  gCCAGGATTTA‑CGTCTGCAGGAAATCCTCGCCCTGGCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGtt                                         <  2:2673245/76‑1 (MQ=255)
                                  gCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCAATTTGTCTGCACGCGTGCGGGCggt                                         >  1:4104350/1‑74 (MQ=255)
                                  gCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GCCGGCCGCCGCATTTTTTCTGCCCGCGGGCGGGCGta                                         >  2:5760847/1‑75 (MQ=255)
                                  cgcAGGTTTTACCGTCGGCAGAATAACCGCGCGCTGGCT‑GCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGt                                          <  2:740443/74‑1 (MQ=255)
                                   ccAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGGCTGCCGGGATTACTCTGTACACGCGTTGGCgccg                                        >  2:8161023/1‑73 (MQ=255)
                                     aGGTTTTA‑CGTCTGCACAATATCCCCGCCCGGCCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtata                                      <  2:4813583/76‑1 (MQ=255)
                                     aGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCGCCATTCGGCTGCCCGCGTGCTGGCGTttta                                      >  1:73995/1‑76 (MQ=255)
                                      ggttttA‑CGTCGGCACAATAACCTCGCCCTGCCT‑GCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatat                                     <  2:6771199/73‑1 (MQ=255)
                                       gttttG‑CGTCTGCAGGATAATCGCGCCCGGGCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTtatata                                    <  2:5688089/74‑1 (MQ=255)
                                        ttttA‑CGTCGCCAGAATACTCCCGCACTGGCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTatatat                                   <  1:1630563/75‑1 (MQ=255)
                                        aTTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGCCCCGCCGCAATTCGTCTGCACGCGTGCCGGCGTtataatt                                   >  1:3452567/1‑73 (MQ=255)
                                        aTTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGTCTGCCCGCGCGCTGGCGCTatatat                                   >  1:6217807/1‑76 (MQ=255)
                                         tttA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGGGGCCGCCGCCAATCCGCTTCCCACGGGCTGGGGTTttatatt                                  >  2:4067564/1‑75 (MQ=255)
                                         tttA‑CGTCGGCACGATATTCCCGCCCGGGCT‑GCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc                                  <  2:4324710/76‑1 (MQ=255)
                                          ttA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCCCAATTTGTCTGCCCGCGGGGTGGGGGTTTATATcc                                 >  1:3290086/1‑75 (MQ=255)
                                          ttA‑CGTCGCCACGATATTCCCCCCCTGGCT‑CCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCg                                 <  1:8044529/76‑1 (MQ=255)
                                           tA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCTGGCGCCATTCCTCTTCACGCGTGCTGGCGTTATTTATTGc                                >  1:8084352/1‑76 (MQ=255)
                                           tA‑CGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGTCTGCCCGCGGGCTGGCGGTATATATCGc                                >  1:187817/1‑76 (MQ=255)
                                              cGTCTGCACGATAATCTCGCCCGGCCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTa                              <  1:3267195/76‑1 (MQ=255)
                                              cGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGTCTGCACGCGTGGTGGGGTTATATATCGCTa                              >  1:5090124/1‑76 (MQ=255)
                                               gTCTGCACGATAATCTCGCACTGGCT‑GGCGCCTGCCGCAATTCGTCTGCACGCGGGCTGGCGGTTTTTATCGCtcc                             >  1:4025114/1‑74 (MQ=255)
                                               gTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCAATTCGTCTGCACGCGTGCTGGGGGTATATATTGCTac                             >  1:3111574/1‑76 (MQ=255)
                                                 cTGCACGATAATCTCGCACTGGCT‑GGGGCCCGCCCCCCATCGTCTGCCCGCGGGCTGGCGGTTTATTTTGCTaccc                           >  2:5068507/1‑74 (MQ=255)
                                                 cTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTTGTTTGCACCCGTGCTGGCGGTATATATTTCTacac                           >  1:2056935/1‑76 (MQ=255)
                                                 cTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCCCCACTCTTCTGCACGCGCGCTGGCGTTTTATATTTCTacac                           >  2:6645037/1‑76 (MQ=255)
                                                 cTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCCCCAATCTTCTGCCACCGCGCTGGGGGTATATTACGCCacac                           >  2:1527327/1‑76 (MQ=255)
                                                 cTGCACGATAATCTCGCACTGGCT‑GCCGCCCGCCGCAATTCGTCTGCCCGCGTGCTGGCGGTATATATCGCTacac                           >  1:3432420/1‑76 (MQ=255)
                                                    cACGATAATCTCGCACTGGCT‑GGCGCCCGCCGCAATTCGTCTGCACGCGGGCTGGCGTTATATATTGCTAcaccac                        >  1:603884/1‑76 (MQ=255)
                                                     aCGATATCCCCGCGCTGGCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg                       <  1:5769430/76‑1 (MQ=255)
                                                      cGATAATCTCGCACTGGCT‑GGGGCCCGCCGCAAATCGTCTGCACGCGGGCGGGCGTTTTATATCGCTACACCACGc                      >  2:7728640/1‑76 (MQ=255)
                                                      cGATAATCTCGCACTGGCT‑GCCGCCTGCCGCACTTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACCCCACGc                      >  2:4662915/1‑76 (MQ=255)
                                                        aTAATCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGTCTGCCCCCGTGGTGGCGGTTTATATCGCCACACCAccccc                    >  2:3974696/1‑75 (MQ=255)
                                                        aTAACCTCCCTCGGGCT‑CCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCa                    <  2:960716/76‑1 (MQ=255)
                                                         tAATCTCGCACTGGCT‑GGCGCCCGCCCCCATTCCTCTGCACGCGGGCTGGCGGTATATATCGCTACACCAcccccg                   >  2:4577808/1‑74 (MQ=255)
                                                          aaTCTCGCACTGGCT‑GGCGCCCGCCGCCATTCGTCTGCACGCGCGCTGGCGGTTTATATCGCTACACACCGcccgc                  >  2:7251726/1‑73 (MQ=255)
                                                          aaTCTCGCACTGGCT‑GGCGCCCGCCGCCAATCCTCCGCACGCGGGCCGGCGGTTTATATATCTACACCACGcccca                  >  2:2622147/1‑73 (MQ=255)
                                                           ccctcGCGCGGCCT‑GCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTg                 <  2:4640098/74‑1 (MQ=255)
                                                             ctcGCGCTGGCT‑GCCCCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTggg               <  2:7247817/76‑1 (MQ=255)
                                                                   cTGGCT‑GCCGCCTGCCGCAATTCGTCTGCACGCGGGCTGGCGTTATATATCGCTACACCACGCCAGTGGGGGACGc         >  2:1226070/1‑76 (MQ=255)
                                                                    tGGCT‑GCCGCCCGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGCCGCa        >  2:7956085/1‑76 (MQ=255)
                                                                       ct‑gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa     <  2:7866114/76‑1 (MQ=255)
                                                                        t‑gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:5729648/1‑76 (MQ=255)
                                                                        t‑gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  2:4173633/1‑76 (MQ=255)
                                                                        t‑gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc    >  1:2032072/1‑76 (MQ=255)
                                                                        ccgccccctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa     <  2:2526483/74‑1 (MQ=255)
                                                                          gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  2:2302074/1‑76 (MQ=255)
                                                                          gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  2:1132998/1‑76 (MQ=255)
                                                                          gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa   >  1:439680/1‑76 (MQ=255)
                                                                           ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  2:3294950/1‑76 (MQ=255)
                                                                           ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  2:2262519/1‑76 (MQ=255)
                                                                           ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt  >  1:4162156/1‑76 (MQ=255)
                                                                           |                                                                           
AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTA‑CGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT  >  REL606/720391‑720539

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