Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,180,2490TG40.8% 227.4 / 55.7 147G24G (GGT→GGGyjaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/54);  new base G (60/0);  total (93/54)
Fisher's exact test for biased strand distribution p-value = 3.24e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 268 total aligned reads displayed.

GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324
                                                                           |                                                                            
ggCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGt                                                                              <  2:3478602/76‑1 (MQ=255)
 gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTg                                                                             <  2:93920/76‑1 (MQ=255)
   aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                                                                           <  1:7762693/76‑1 (MQ=255)
     tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCt                                                                         >  2:140652/1‑76 (MQ=255)
      gAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTg                                                                        >  2:6632955/1‑76 (MQ=255)
         cAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGgtt                                                                     >  1:2675188/1‑76 (MQ=255)
          aGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtc                                                                    <  1:6983410/76‑1 (MQ=255)
            ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTtttc                                                                  >  2:4834428/1‑76 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCa                                                                 <  1:8145209/76‑1 (MQ=255)
              ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGGTGGTtttccc                                                                >  1:3314791/1‑74 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTTACt                                                               >  2:3209667/1‑76 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACt                                                               >  2:1716636/1‑76 (MQ=255)
                aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTg                                                              <  1:8081539/76‑1 (MQ=255)
                aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTctctc                                                              >  1:148600/1‑75 (MQ=255)
                  cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCGTTTTTTACTAcc                                                            >  2:3594685/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                           <  2:5183991/76‑1 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                           <  2:4081618/76‑1 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTTccc                                                           >  2:2731669/1‑75 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTg                                                          <  1:1043010/76‑1 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGt                                                         <  2:7906228/76‑1 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTTACTGcccgg                                                         >  2:6188779/1‑73 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGCCCGt                                                         >  1:3932082/1‑76 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTTTCACTGCCTGTa                                                        >  1:6467894/1‑76 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                        <  2:2153954/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                        <  1:4435648/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                        <  2:7455454/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTa                                                        >  2:1330804/1‑76 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                       <  1:5843920/76‑1 (MQ=255)
                          tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCca                                                    <  2:5060413/76‑1 (MQ=255)
                           ggtTTCCCCCCTTTTGTTTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTTTAGCcac                                                   <  2:3913831/76‑1 (MQ=255)
                           ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCCGTAGcccc                                                   >  2:860756/1‑74 (MQ=255)
                            gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTTCTTACCGCCTGTAGCcaca                                                  >  2:8467018/1‑76 (MQ=255)
                            gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTCTCACTGCCTGTCGCcaca                                                  >  1:3473667/1‑76 (MQ=255)
                             tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTTCCTCTAGCCACaa                                                 >  2:2697878/1‑76 (MQ=255)
                              ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGCCTGTAGGCACAAc                                                >  1:5027825/1‑76 (MQ=255)
                               tccccccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAAct                                               <  1:659493/73‑1 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTTCTTCCTGCCTTTAGCCACAAct                                               >  1:6325102/1‑76 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTTTCACTGCCTGTAGCCACaaac                                               >  1:670598/1‑74 (MQ=255)
                                cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTTCATGCCTGTAGCCCCAActc                                             >  1:7192437/1‑76 (MQ=255)
                                 aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCAGTTTTTCACTGCCTGCAGCCCCAActct                                            >  2:8239837/1‑76 (MQ=255)
                                    ctttttGTCTTTCTGGCGTCGATCATTGATGCTGGCGGGTGCGCTTTTTCTCACTGCCTTTAGCCACAACTCTTCa                                         <  2:8917711/74‑1 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCa                                         <  1:6191286/76‑1 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCa                                         >  1:2011513/1‑76 (MQ=255)
                                     cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGCCACAACTCTTccc                                        >  2:5431848/1‑74 (MQ=255)
                                      ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACt                                       <  1:1332859/76‑1 (MQ=255)
                                      ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTTCCTGTAGCCCCCACCCCTCACt                                       >  1:2275759/1‑76 (MQ=255)
                                       tttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACtt                                      >  1:414961/1‑76 (MQ=255)
                                         tGTTTTTCGGGCGTCGATCATTGATGCTGGCGGGCGCGCTGTTCCTCCCTGCTTGTAGCCACAACTCTTCACTtcc                                    <  2:5260221/76‑1 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGCCACCACCCTTCACTtcc                                    >  2:4660130/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTtcc                                    >  1:850568/1‑76 (MQ=255)
                                          ttCTTTCGGGCGTCGATCATTGATGCTGGCAGGTGCGCTTTTTCCCCCTGCCTTTAGCCACAACTCTTCACTtcct                                   <  1:984997/75‑1 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTTCCTGTAGCCCCAACACTTCACTtcccc                                  >  1:4242584/1‑74 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTCACCACAACTCTTCACTtcctc                                  >  1:2093676/1‑76 (MQ=255)
                                           tCTTTCTGGCGTCGACCAGTGATGCTGGCAGGTGCGTTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTtcctc                                  <  1:3157155/76‑1 (MQ=255)
                                            ttttCGGGGGTCGTTAATGGATGTGGGCGGGTGCGCTTTTTCTCCCTGCTTTTAGCCACAACTCTTCACTtcctcc                                 <  2:3349232/75‑1 (MQ=255)
                                            cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGGCCGTAGCCACAACTCTTCACTttcttc                                 >  1:5176316/1‑74 (MQ=255)
                                            cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGACACAAATCTTCACCtcctcc                                 >  2:2507802/1‑76 (MQ=255)
                                            cTTTCTGGCGTAGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAAATCTTCACTtcctcc                                 >  2:320485/1‑76 (MQ=255)
                                             tttCGGGCGTCGATCATTGATGCGGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTccc                                <  1:3019903/76‑1 (MQ=255)
                                              ttgtGGCGTCTTTCTTTGATGCTGGCAGGTGCGTTTTTTCTCACTCCCTTTAGCCACAACTCTTCACTTCCTCCCt                               <  2:4270489/73‑1 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCt                               >  1:5988687/1‑76 (MQ=255)
                                               tCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACCGCCTGTAGCCACAACTCTTTACTTCCTCCCtt                              >  2:5875941/1‑76 (MQ=255)
                                               tCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCCGTAGCCACAACTCTTCACTTCCTCCCtt                              >  2:6149089/1‑76 (MQ=255)
                                                cTGGCGTGGATTTTTGATGCTGGCGGGTGCTCTTTTTCTCCCTTCCTGTAGCCACAACTCTTCACTTCCTCCCttt                             <  2:7288994/76‑1 (MQ=255)
                                                cTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTTGCCACAACTCTTCACTCTCACCCttt                             >  2:4221268/1‑76 (MQ=255)
                                                cTGGCGTCGATCATTGATGCTGGCAGGGGGGCGGTTTTTCACTGCCTGTAGCCACCACTCTTCACTTCCTCCCctt                             >  1:7610387/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGTGGGCTGTTTCTTACTGCCTGTTGGCAAAACTCTTCACCTCCTCCCtttt                            >  1:7403111/1‑75 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCcttt                            >  2:7329269/1‑75 (MQ=255)
                                                  ggCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTTACTGCCTGTAGCCCCAAATCCTCCACTTCTCCCttttc                           >  1:6776311/1‑74 (MQ=255)
                                                  ggCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGGCACAAATCTTCACTTCCTCCCTTTAc                           >  2:771558/1‑76 (MQ=255)
                                                   gCGTGGATCTTTAAGCCTGGCGGGTCCCTTTTTTCCCCCTCCTTGTAGCCACAACTCTTCACTTCCTCCCTTTAcc                          <  1:7370165/76‑1 (MQ=255)
                                                   gCGTCGATCTTTGCTGCTGGGAGGTGCGTTGTTTCCCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTAcc                          <  1:1853130/76‑1 (MQ=255)
                                                    cGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGCCCGTAGCCACAAATCTTCACTTCCTCCCTTTTccc                         >  2:7961844/1‑75 (MQ=255)
                                                     gTCGATTTTTGATGCTGGCGGGGGCGCGTTTTCCCCCTGCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGc                        <  2:7640782/76‑1 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTGCCCGTAGGCACAACTCTTCACTTTCTCCCTTTACCGc                        >  1:1806155/1‑76 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCa                      >  2:3091773/1‑76 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGccc                      >  1:4677411/1‑75 (MQ=255)
                                                        gATCATTGATGCTGGCAGGGGCGCTGTTTTTCACCGCCTGTAGCCACCAACCTTCACTTTCTCCCCTTTCCCCCac                     >  2:594703/1‑75 (MQ=255)
                                                            aTTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGa                 >  1:4823569/1‑76 (MQ=255)
                                                              tGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTTCTCCCTTTACCGCCCGTGGAtt               >  1:7973826/1‑76 (MQ=255)
                                                              tGATGCTGGCAGGGGCGCGGTTTTTCACTTCCCGTGGCGCCCACACTTCCACTCCTCCCCTTACCCCCCGCGGAtt               >  2:4706375/1‑76 (MQ=255)
                                                                   cTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGccg          >  2:4018688/1‑74 (MQ=255)
                                                                   cTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  1:6029128/1‑76 (MQ=255)
                                                                   cTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  1:7460086/1‑76 (MQ=255)
                                                                   cTGGCAGGGGCGCTGTTTCTCACTGCCTATAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  2:1544583/1‑76 (MQ=255)
                                                                    tGGCAGGTGCGCTGTTTCTCACTGCCTGTAGGCACAACTCTTCACTTCCTCCCTTTACCGCCAGCGGATTTGCTGa         >  2:2016161/1‑76 (MQ=255)
                                                                    tGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGa         >  2:8421196/1‑76 (MQ=255)
                                                                    tGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGa         >  2:2925360/1‑76 (MQ=255)
                                                                     ggCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGaa        >  2:2298007/1‑76 (MQ=255)
                                                                     ggCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGaa        >  2:7174762/1‑76 (MQ=255)
                                                                     ggCAGGGGCGCCTTTTCGCACTGCTTGTAGCTCCGACTCTTCACTTCCTCCTCTTACGCCCAGTGGACTTGCTgca        >  2:4857252/1‑74 (MQ=255)
                                                                      gCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAg       >  2:4555301/1‑76 (MQ=255)
                                                                      gCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAg       >  1:7052992/1‑76 (MQ=255)
                                                                        aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:3838301/1‑76 (MQ=255)
                                                                         ggTGCGTTTTTCCTCCCGGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc    <  1:8810470/76‑1 (MQ=255)
                                                                          gTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCa   >  1:1426217/1‑76 (MQ=255)
                                                                           tGCTCTGTTTCCCCCCCCTTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  <  1:1518438/76‑1 (MQ=255)
                                                                           tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:8870246/1‑76 (MQ=255)
                                                                           tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:8374273/1‑76 (MQ=255)
                                                                           tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:5439938/1‑76 (MQ=255)
                                                                           |                                                                            
GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324

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