Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,180,2490TG40.2% 180.4 / 35.5 107G24G (GGT→GGGyjaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/45);  new base G (42/1);  total (61/46)
Fisher's exact test for biased strand distribution p-value = 1.25e-13
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 193 total aligned reads displayed.

GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324
                                                                            |                                                                           
   aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                                                                          <  1:2321774/76‑1 (MQ=255)
     tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                        <  2:1166385/76‑1 (MQ=255)
       aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGt                                                                      >  2:347676/1‑76 (MQ=255)
          aGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtc                                                                   <  2:3746943/76‑1 (MQ=255)
           gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtct                                                                  <  2:761793/76‑1 (MQ=255)
           gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtct                                                                  <  1:3334381/76‑1 (MQ=255)
           gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtct                                                                  <  1:2668720/76‑1 (MQ=255)
           gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtct                                                                  <  2:5966623/76‑1 (MQ=255)
           gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTtct                                                                  >  1:2033031/1‑76 (MQ=255)
            ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtctc                                                                 <  2:623029/76‑1 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTtttct                                                                >  1:3062011/1‑75 (MQ=255)
              ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  1:6545081/76‑1 (MQ=255)
                  cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGcc                                                           >  1:1717680/1‑76 (MQ=255)
                  cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTccc                                                           >  1:1407582/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                          >  2:6971992/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                          <  1:6918378/76‑1 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCt                                                          >  2:1637924/1‑76 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTg                                                         <  2:3366889/76‑1 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTTcccg                                                         >  2:4177284/1‑74 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTg                                                         >  2:6857432/1‑76 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGt                                                        <  1:471272/76‑1 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGt                                                        <  1:1760566/76‑1 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACCGCCTGt                                                        >  2:1581064/1‑76 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                       <  1:261386/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                       <  2:6108167/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTa                                                       >  1:3435236/1‑76 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTTCCTGTa                                                       >  2:2166879/1‑76 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  2:2339773/76‑1 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  1:96184/76‑1 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  1:3745742/76‑1 (MQ=255)
                        ggtggtTTCACGCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGc                                                     <  2:203439/76‑1 (MQ=255)
                         gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGcc                                                    <  1:4170618/76‑1 (MQ=255)
                         gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCGGTTTTTACTGCCCGTTGcc                                                    >  1:5129521/1‑76 (MQ=255)
                          tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGGGGTGGTTTTCTCTTCCTTTCGccc                                                   >  1:5502194/1‑75 (MQ=255)
                           ggtTTCCCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCcac                                                  <  1:6796486/76‑1 (MQ=255)
                           ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCcac                                                  <  1:6170855/76‑1 (MQ=255)
                           ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGGcac                                                  >  1:4817563/1‑76 (MQ=255)
                            gtTTCCCCCTTTTTGTTTTTCTGGCGTCGATCATTGATTCTGGCAGGTGCTCTGTTTCTCCCTGCCTTTAGCcaca                                                 <  2:1493611/76‑1 (MQ=255)
                              ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGTGCTGTTTCTCACTTCCTGTAGCCACAAc                                               >  2:2750712/1‑76 (MQ=255)
                              ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAAc                                               >  2:4794508/1‑76 (MQ=255)
                               tccccctTTTTGTCTTTCTGGCGCCGATCTTTGATGCTGGCGGGTGCCTTGTTTCTCCCTGCTTTTAGCCACAAct                                              <  2:3035055/73‑1 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCTCTGCCTGTAAGcacacat                                              >  2:5649020/1‑73 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGCCTGTAGCCACaaac                                              >  1:5925626/1‑74 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACaaat                                              >  2:1973362/1‑74 (MQ=255)
                                cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTCCGCGTAGCCCCaaacc                                             >  1:2189435/1‑73 (MQ=255)
                                cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCcaccaccc                                             >  2:5113273/1‑74 (MQ=255)
                                 aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACCActct                                            >  1:1648417/1‑76 (MQ=255)
                                 aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAAcact                                            >  2:4058311/1‑76 (MQ=255)
                                  ccctTTTTGTCTTTCGGGCGTCGATCATTGATGCTGGCAGGTGCGCTTTTCCTCACTGCCTGTAGCCACAACTCtt                                           <  2:5182204/76‑1 (MQ=255)
                                  ccccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCCGTAGGCCCAACCCtt                                           >  1:5542023/1‑76 (MQ=255)
                                  ccccTTTTGTCTTTCGGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTTTAGCCACAACTCtt                                           <  2:4505269/76‑1 (MQ=255)
                                   cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTc                                          <  1:904771/76‑1 (MQ=255)
                                   cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTTGCCACAACTCTTc                                          >  1:6610465/1‑76 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTTCTCACTGCCTGTAGCCACAACTCTTcc                                         >  2:5051239/1‑75 (MQ=255)
                                     cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTTGCCACAACTCTTCAc                                        >  1:2683744/1‑76 (MQ=255)
                                     cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGGCACAACTCCTCAc                                        >  2:862220/1‑76 (MQ=255)
                                     cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAc                                        >  2:6904990/1‑76 (MQ=255)
                                      ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGAGCGCTCTTTCTCCGTGCCCGTAGCCACAACTCCTCACt                                       >  2:6172798/1‑76 (MQ=255)
                                      ttttGCCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACt                                       <  1:2550720/76‑1 (MQ=255)
                                       tttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTTCCTTTAGCCACAActcttctctt                                      >  2:6387067/1‑76 (MQ=255)
                                        ttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCAATGCCTGTCGACACAACTCTGCAGTtc                                     >  1:5323114/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTtcc                                    >  2:3025844/1‑76 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGTTTTTTACCGTCTGTAGCAACAACTCATCACTccctc                                  >  1:1046448/1‑76 (MQ=255)
                                            cTTTCTGGCGTCGATCATTGATGCTGGCAGGCGCGCCGTTTTTCCCTCCCTCTAACCACAAACCTTCACAtccacc                                 >  1:379359/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTCCTccc                                >  1:5216809/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTccc                                <  1:1450492/76‑1 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAAATCTTCACTTCCTCCCt                               >  1:5352339/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACCGCCCCTAGGCACAACACTTCACTTCCTCCCt                               >  1:3976784/1‑76 (MQ=255)
                                              ttCGGGCGTCGTTCATTGATGCGGGGGGGTGCGCGTTTTCTCGCGGCTTTTAGCCACAACTCTTCACTTCCTCCCt                               <  2:3259467/76‑1 (MQ=255)
                                               tCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCtt                              >  2:6320260/1‑76 (MQ=255)
                                                cTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCttt                             >  1:2978034/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGTGCGCTTTTTCTCACTGCGTGTAGCCACAACTCTTCACTTCCTCCCTTTa                            <  1:1224708/76‑1 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGGCACAACTCTTCACTTCCTCCCTTTa                            >  1:2836598/1‑76 (MQ=255)
                                                 gggCGTCGATCATTGATGCGGGCGGGGGCGTTTTTTCCCACTCCCTTTAGCCACAACTCTTCACTTCCTCCCTTTa                            <  2:2478799/75‑1 (MQ=255)
                                                  tgCGTCGGTCATTGATGCTGGGCGGGGGCGTTTTTCCCCCCTCCCTTTAGCCACAACTCTTCACTTCCTCCCCTTa                            <  1:3527043/75‑1 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCCCAAATCTTCACTTCCTCCCTTTAcacc                        >  1:2495374/1‑73 (MQ=255)
                                                      tCGCTCATTGATGCTGGGGGGTGCTCTTTTTCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGcc                       <  1:6431282/76‑1 (MQ=255)
                                                      tCGATCATTGATCCTGGCAGGTGCGCTGTTTCCCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGcc                       <  1:5684100/76‑1 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGTGCGCTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCa                      <  1:2266803/76‑1 (MQ=255)
                                                        gATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAg                     >  1:4383761/1‑76 (MQ=255)
                                                         aTCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTTCCGCCCGt                    >  2:5476339/1‑76 (MQ=255)
                                                          tCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGGCAAGAATCTTCACTTCCTCCCTTTACCGCCACTg                   >  2:3902035/1‑76 (MQ=255)
                                                          tCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTg                   >  2:6513326/1‑76 (MQ=255)
                                                             ttGATGTGGGCGGGTGCCTTGTTTCTCCCTGCTTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAt                <  1:1082391/76‑1 (MQ=255)
                                                               gATGCTGGCGGGTGCGCTTTTTCTCCCTCCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAttt              <  1:2473790/76‑1 (MQ=255)
                                                                aTGCTGGCAGGGGGGCTGTTTCTCACCGCCTGTAGCCACAACTCTTCCCTTCCTCCCTTTACCCCCCGTGGAtttt             >  2:1854186/1‑75 (MQ=255)
                                                                aTGCTGGCAGGGGGGCCGTTTCTCACCTCCTGTAGGCACAACTCTTCCCTTCCTCCCTTTACCGCCAGTGGATTTg             >  2:559733/1‑76 (MQ=255)
                                                                      ggCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGaa        >  1:347114/1‑76 (MQ=255)
                                                                       gCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAg       >  1:5647058/1‑76 (MQ=255)
                                                                        cgggTGCGTTTTTCCCCCCTGCTGGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa      <  1:5364288/74‑1 (MQ=255)
                                                                        cAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa      >  1:6597141/1‑76 (MQ=255)
                                                                         aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:2124728/1‑76 (MQ=255)
                                                                         aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:10982/1‑76 (MQ=255)
                                                                         aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:4522605/1‑76 (MQ=255)
                                                                         aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:367573/1‑76 (MQ=255)
                                                                         aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  2:1694332/1‑76 (MQ=255)
                                                                          ggTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc    >  1:4440412/1‑76 (MQ=255)
                                                                            tGCGTTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  <  1:5068590/76‑1 (MQ=255)
                                                                            tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:4090893/1‑76 (MQ=255)
                                                                            tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  2:1913044/1‑76 (MQ=255)
                                                                            |                                                                           
GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324

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