Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,338,860 A→C R230S (AGA→AGC yhdG → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,338,8600AC99.2% 1093.8 / ‑5.0 264R230S (AGA→AGCyhdGtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (103/159);  total (104/160)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

CCCGCTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGACCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCAGAG  >  REL606/3338785‑3338932
                                                                           |                                                                        
  cGCTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGccc                                                                        <  1:2169186/76‑1 (MQ=255)
   gCTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGcccc                                                                       >  1:3566607/1‑76 (MQ=255)
    cTTAACGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCCCAGCTCAGGGAAGCCCCt                                                                      <  1:529894/76‑1 (MQ=255)
    cTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCCGGGAAGCCCCt                                                                      >  1:1360314/1‑76 (MQ=255)
    cTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCt                                                                      <  1:5772299/76‑1 (MQ=255)
    cTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCt                                                                      <  1:1261973/76‑1 (MQ=255)
    cTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCt                                                                      <  2:4714422/76‑1 (MQ=255)
     ttAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATTATAAGGAACGCCGCTCAGGGAACCCCCtc                                                                     >  2:5399173/1‑75 (MQ=255)
      tAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTgg                                                                    <  2:1137035/76‑1 (MQ=255)
      tAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTgg                                                                    <  1:3843776/76‑1 (MQ=255)
      tAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTgg                                                                    >  1:4834080/1‑76 (MQ=255)
       aaaGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGa                                                                   <  1:2016616/76‑1 (MQ=255)
       aaaGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGa                                                                   >  1:6353222/1‑76 (MQ=255)
        aaGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGAt                                                                  <  2:1308984/76‑1 (MQ=255)
        aaGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGAt                                                                  <  2:6851786/76‑1 (MQ=255)
         aGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCCGGGAAGCCCCTGGATc                                                                 >  2:5925281/1‑76 (MQ=255)
            cAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCttt                                                              >  2:1735779/1‑76 (MQ=255)
            cAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCttt                                                              >  2:7025376/1‑76 (MQ=255)
             agagCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTc                                                             <  2:6544407/76‑1 (MQ=255)
             agagCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTc                                                             <  1:625824/76‑1 (MQ=255)
             agagCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTc                                                             <  1:1855243/76‑1 (MQ=255)
                gCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCggg                                                          >  1:5565238/1‑76 (MQ=255)
                gCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCggg                                                          <  1:1596930/76‑1 (MQ=255)
                 cTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGa                                                         >  1:3584750/1‑76 (MQ=255)
                 cTGTGCTCGACTATACAGGGGCGGAAGCCCCGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCCTTCgggg                                                         >  1:4598942/1‑75 (MQ=255)
                   gtgCTCGACTATACAGGGGCGGATGCCCTGATGATAGGGCGCGCCACTCCGGGAAGGCCCTGGATCTTTCCGGgaa                                                       >  2:6496714/1‑76 (MQ=255)
                   gtgCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGaaa                                                       >  1:2330875/1‑76 (MQ=255)
                   gtgCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGaaa                                                       >  1:6423064/1‑76 (MQ=255)
                    tgCTCGACTATACAGGTGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAAt                                                      >  1:4040421/1‑76 (MQ=255)
                    tgCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGGAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAAt                                                      >  1:6089069/1‑76 (MQ=255)
                    tgCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAAt                                                      <  1:4002934/76‑1 (MQ=255)
                       tCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCa                                                   <  1:5917904/76‑1 (MQ=255)
                       tCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGCGAAATCCa                                                   <  1:4173624/76‑1 (MQ=255)
                        cGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAg                                                  <  2:2133597/76‑1 (MQ=255)
                        cGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAg                                                  <  1:4330743/76‑1 (MQ=255)
                         gACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATACAGc                                                 >  1:5341992/1‑76 (MQ=255)
                          aCTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCa                                                <  2:419459/76‑1 (MQ=255)
                          aCTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCa                                                <  2:6631229/76‑1 (MQ=255)
                           cTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCAt                                               <  2:817761/76‑1 (MQ=255)
                            tataCAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCAtt                                              <  1:2265296/76‑1 (MQ=255)
                             ataCAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTa                                             <  1:56126/76‑1 (MQ=255)
                             ataCAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTa                                             >  2:6649378/1‑76 (MQ=255)
                             ataCAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTa                                             <  2:5860743/76‑1 (MQ=255)
                              taCAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTAt                                            <  1:2945532/76‑1 (MQ=255)
                                cAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCt                                          >  1:5429651/1‑76 (MQ=255)
                                cAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCt                                          <  1:6565010/76‑1 (MQ=255)
                                cAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCt                                          <  2:4064059/76‑1 (MQ=255)
                                 aGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTg                                         <  1:6127212/76‑1 (MQ=255)
                                 aGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTg                                         <  2:1497629/76‑1 (MQ=255)
                                   gggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGa                                       <  2:4868479/76‑1 (MQ=255)
                                   gggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGa                                       <  1:728076/76‑1 (MQ=255)
                                   gggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGa                                       <  2:2366220/76‑1 (MQ=255)
                                    ggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGac                                      >  1:824731/1‑76 (MQ=255)
                                    ggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGac                                      >  1:5422678/1‑76 (MQ=255)
                                    ggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGac                                      <  1:2291176/76‑1 (MQ=255)
                                    ggCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGac                                      <  1:3126858/76‑1 (MQ=255)
                                      cGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGacac                                    >  1:1799086/1‑76 (MQ=255)
                                       ggATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACt                                   <  2:3747037/76‑1 (MQ=255)
                                         aTGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTgg                                 >  2:5574817/1‑76 (MQ=255)
                                          tGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTggg                                >  1:750291/1‑76 (MQ=255)
                                           gCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTgggg                               >  1:1515691/1‑76 (MQ=255)
                                           gCCCTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTgggg                               >  2:4890854/1‑76 (MQ=255)
                                            cccTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGa                              <  2:5791506/76‑1 (MQ=255)
                                            cccTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGa                              <  2:2533691/76‑1 (MQ=255)
                                             ccTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAg                             >  2:4512570/1‑76 (MQ=255)
                                             ccTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAg                             <  1:4521126/76‑1 (MQ=255)
                                              cTGATGATAGGCCGCGCAGCTCAGGGAAGCGCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGt                            >  2:3888146/1‑76 (MQ=255)
                                              cTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGt                            >  1:5218197/1‑76 (MQ=255)
                                              cTGATGATAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGt                            >  2:4565205/1‑76 (MQ=255)
                                               tgatgaTAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGtt                           >  2:2592694/1‑76 (MQ=255)
                                                 atgatAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGATACTGGGGAGTtgc                         >  1:3398441/1‑76 (MQ=255)
                                                 atgatAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTtgc                         <  2:4958761/76‑1 (MQ=255)
                                                  tgatAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTtgct                        >  1:3701913/1‑76 (MQ=255)
                                                    atAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTtgctgc                      >  1:4089139/1‑76 (MQ=255)
                                                    atAGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTtgctgc                      >  2:4010811/1‑76 (MQ=255)
                                                      aGGCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGccc                    <  2:4296060/76‑1 (MQ=255)
                                                        gCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGccccc                  <  2:5733752/76‑1 (MQ=255)
                                                        gCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGccccc                  <  1:1892477/76‑1 (MQ=255)
                                                        gCCGCGCAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGccccc                  <  1:268975/76‑1 (MQ=255)
                                                           gcgcAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCt               <  2:2697198/76‑1 (MQ=255)
                                                            cgcAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTg              <  1:6797181/76‑1 (MQ=255)
                                                            cgcAGCCCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTg              <  1:3344179/76‑1 (MQ=255)
                                                              cAGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGcc            >  1:885684/1‑76 (MQ=255)
                                                               aGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCt           >  2:3940958/1‑76 (MQ=255)
                                                               aGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCt           >  2:337198/1‑76 (MQ=255)
                                                               aGCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCt           >  2:5986893/1‑76 (MQ=255)
                                                                gCTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCtt          <  2:3551048/76‑1 (MQ=255)
                                                                 cTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCttt         <  2:5247036/76‑1 (MQ=255)
                                                                 cTCAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCttt         <  1:1011216/76‑1 (MQ=255)
                                                                   cAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTgg       >  2:4822275/1‑76 (MQ=255)
                                                                   cAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTgg       >  1:2535818/1‑76 (MQ=255)
                                                                   cAGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTgg       >  2:6822776/1‑76 (MQ=255)
                                                                    aGGGAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGc      >  2:3079842/1‑76 (MQ=255)
                                                                     gggAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCa     <  1:2820597/76‑1 (MQ=255)
                                                                      ggAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCag    <  1:4656222/76‑1 (MQ=255)
                                                                      ggAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCag    <  2:1997552/76‑1 (MQ=255)
                                                                      ggAAGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTNGGCag    <  2:2253035/76‑1 (MQ=255)
                                                                        aaGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCagag  <  1:1043286/76‑1 (MQ=255)
                                                                        aaGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCagag  <  2:6321934/76‑1 (MQ=255)
                                                                        aaGCCCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCagag  <  2:4622584/76‑1 (MQ=255)
                                                                           |                                                                        
CCCGCTTAAAGCCAGAGCTGTGCTCGACTATACAGGGGCGGATGCCCTGATGATAGGCCGCGCAGCTCAGGGAAGACCCTGGATCTTTCGGGAAATCCAGCATTATCTGGACACTGGGGAGTTGCTGCCCCCGCTGCCTTTGGCAGAG  >  REL606/3338785‑3338932

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