Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,927,357 (T)8→9 coding (74/1266 nt) rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,927,3491.T96.5% 565.1 / 11.7 142coding (82/1266 nt)rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (47/89);  total (49/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 162 total aligned reads displayed.

CCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATG  >  REL606/3927274‑3927418
                                                                             |                                                                     
ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAccc                                                                        >  2:4726524/1‑76 (MQ=255)
ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAccc                                                                        >  2:1029955/1‑76 (MQ=255)
   gCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttt                                                                    >  1:826125/1‑76 (MQ=255)
   gCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttt                                                                    >  1:1690810/1‑76 (MQ=255)
     cgggCTACGCCCCGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttt                                                                  <  1:56693/74‑1 (MQ=255)
     ccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttt                                                                  >  2:14466/1‑76 (MQ=255)
       ggCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttttt                                                                >  1:2722596/1‑76 (MQ=255)
        gCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttt                                                               >  2:6881423/1‑76 (MQ=255)
        gCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGGGCCTGAATGGGCGTACAGTTATGAAACCCtttttttt                                                               >  1:2268145/1‑76 (MQ=255)
              tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAgg                                                          >  2:1572634/1‑76 (MQ=255)
              tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAgg                                                          >  2:1749747/1‑76 (MQ=255)
               cACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAACCTTTTTTTTTCAAggg                                                         >  1:1247065/1‑76 (MQ=255)
                 ccccccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCt                                                       <  2:4025700/73‑1 (MQ=255)
                  cccccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCtt                                                      <  2:2354830/74‑1 (MQ=255)
                  cccccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCta                                                      <  2:3971746/74‑2 (MQ=255)
                   aCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTc                                                     >  2:2143215/1‑76 (MQ=255)
                    cccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCt                                                    <  1:1064308/76‑1 (MQ=255)
                     ccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTa                                                   <  2:1900192/76‑1 (MQ=255)
                       gCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACa                                                 <  2:5085468/76‑1 (MQ=255)
                       gCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACa                                                 <  2:510326/76‑1 (MQ=255)
                        ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACaa                                                <  2:2262610/76‑1 (MQ=255)
                         cAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAc                                               >  2:2638942/1‑76 (MQ=255)
                         cAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAc                                               <  1:517334/76‑1 (MQ=255)
                         cAGCGTCAGCGGAAGGCCCAGTGCCTGAATGGGCGTACAGTTATGAATCCCTTTTTTTTTCAAGGGCTTCTACAAc                                               <  2:6154789/76‑1 (MQ=255)
                            cGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAACtt                                            >  2:4145097/1‑76 (MQ=255)
                             gTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTc                                           >  2:4904631/1‑76 (MQ=255)
                             gTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTc                                           >  1:4101571/1‑76 (MQ=255)
                              tCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCg                                          <  1:1876138/76‑1 (MQ=255)
                              tCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCg                                          >  1:6423533/1‑76 (MQ=255)
                               cAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCgg                                         <  2:1012590/76‑1 (MQ=255)
                               cAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCgg                                         <  2:3825800/76‑1 (MQ=255)
                                aGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGa                                        <  1:5604075/76‑1 (MQ=255)
                                 gCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACACCCTTCGGAt                                       <  2:5659088/76‑1 (MQ=255)
                                   gggaGGTCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                     <  2:1394253/73‑1 (MQ=255)
                                   ggAGGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                     <  1:2005197/76‑1 (MQ=255)
                                   ggAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                     <  2:4494453/76‑1 (MQ=255)
                                    gAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                    <  2:1191547/76‑1 (MQ=255)
                                    gAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                    >  1:5489984/1‑76 (MQ=255)
                                    gAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                    <  2:3299789/76‑1 (MQ=255)
                                      aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                  <  1:1885667/76‑1 (MQ=255)
                                      aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                  <  1:4575610/76‑1 (MQ=255)
                                      aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                  <  2:2691419/76‑1 (MQ=255)
                                       gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAggg                                 <  2:888097/76‑1 (MQ=255)
                                       gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAggg                                 <  1:1275123/76‑1 (MQ=255)
                                           cAGTGCCTGAATGGTCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                             <  1:5073691/76‑1 (MQ=255)
                                           cAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGGGaa                             >  2:5553053/1‑76 (MQ=255)
                                           cAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                             >  2:2888618/1‑76 (MQ=255)
                                           cAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                             <  2:2484588/76‑1 (MQ=255)
                                            aGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAg                            >  1:1734752/1‑76 (MQ=255)
                                            aGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAg                            <  1:4732517/76‑1 (MQ=255)
                                             gTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                           <  1:2693407/76‑1 (MQ=255)
                                              tGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc                          >  2:6274746/1‑76 (MQ=255)
                                              tGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc                          >  1:1516060/1‑76 (MQ=255)
                                              tGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc                          >  2:5466860/1‑76 (MQ=255)
                                               gCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCg                         >  1:2713994/1‑76 (MQ=255)
                                               gCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCg                         <  1:4168953/76‑1 (MQ=255)
                                                ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                        <  2:3676847/76‑1 (MQ=255)
                                                ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                        >  2:6258846/1‑76 (MQ=255)
                                                 cTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGa                       <  2:1286026/76‑1 (MQ=255)
                                                  tGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                      <  2:6653896/76‑1 (MQ=255)
                                                  tGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                      <  2:3438725/76‑1 (MQ=255)
                                                  tGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                      >  2:2078688/1‑76 (MQ=255)
                                                    aaTGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                    <  2:275933/76‑1 (MQ=255)
                                                    aaTGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                    <  1:2223642/76‑1 (MQ=255)
                                                    aaTGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                    >  1:3559064/1‑76 (MQ=255)
                                                     aTGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAAc                   >  1:5010223/1‑76 (MQ=255)
                                                      tGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACt                  <  2:4673766/76‑1 (MQ=255)
                                                      tGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACt                  >  2:160501/1‑76 (MQ=255)
                                                        ggCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttc                <  2:4978776/76‑1 (MQ=255)
                                                         gCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct               <  2:5991761/76‑1 (MQ=255)
                                                         gCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct               <  2:4358048/76‑1 (MQ=255)
                                                         gCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttct               <  1:6704039/76‑1 (MQ=255)
                                                          cGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg              <  1:1132328/76‑1 (MQ=255)
                                                          cGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg              <  2:2359255/76‑1 (MQ=255)
                                                          cGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg              >  2:1747090/1‑76 (MQ=255)
                                                             aCAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc           >  2:3395256/1‑76 (MQ=255)
                                                              cAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct          >  1:3067462/1‑76 (MQ=255)
                                                              cAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct          <  2:26814/76‑1 (MQ=255)
                                                              cAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct          <  2:1943232/76‑1 (MQ=255)
                                                               aGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg         >  1:3679409/1‑76 (MQ=255)
                                                               aGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg         <  1:3525456/76‑1 (MQ=255)
                                                               aGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg         <  1:6175280/76‑1 (MQ=255)
                                                                 ttATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt       >  1:6735496/1‑76 (MQ=255)
                                                                 ttATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt       <  2:3245579/76‑1 (MQ=255)
                                                                 ttATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt       >  2:6177573/1‑76 (MQ=255)
                                                                 ttATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt       >  1:6021642/1‑76 (MQ=255)
                                                                 ttATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa      >  1:379397/1‑76 (MQ=255)
                                                                  tATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa      <  1:776003/76‑1 (MQ=255)
                                                                  tATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa      >  1:6019594/1‑76 (MQ=255)
                                                                  tATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa      <  2:4804169/76‑1 (MQ=255)
                                                                  tATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa      >  2:4439175/1‑76 (MQ=255)
                                                                   aTGAACCCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa     <  1:2101468/76‑1 (MQ=255)
                                                                   aTGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa     <  1:4195170/76‑1 (MQ=255)
                                                                   aTGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa     <  1:6515065/76‑1 (MQ=255)
                                                                   aTGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa     <  2:733425/76‑1 (MQ=255)
                                                                    tGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa    >  1:2556708/1‑76 (MQ=255)
                                                                     gAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt   >  1:6227012/1‑76 (MQ=255)
                                                                     gAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt   <  2:3759323/76‑1 (MQ=255)
                                                                       aaaCCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtg  <  1:5929398/76‑1 (MQ=255)
                                                                       aaaCCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtg  >  2:1955837/1‑76 (MQ=255)
                                                                             |                                                                     
CCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATG  >  REL606/3927274‑3927418

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