Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 79,133 (T)5→4 coding (971/1656 nt) yabN ← DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL60679,1290T.99.3% 720.6 / ‑2.1 147coding (975/1656 nt)yabNDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (67/79);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 164 total aligned reads displayed.

CTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCA  >  REL606/79054‑79200
                                                                           |                                                                       
cTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAAActct                                                                         >  2:2867707/1‑76 (MQ=255)
 tGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTCGCAACCTTCCTCCAGGCGGCGTTCAATCTCtt                                                                        >  2:3977567/1‑76 (MQ=255)
 tGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGGTTTCAATCTCtt                                                                        >  1:3729300/1‑76 (MQ=255)
  gATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCttt                                                                       <  1:3979214/76‑1 (MQ=255)
   aTTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCtttt                                                                      <  2:6039712/76‑1 (MQ=255)
   aTTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCtttt                                                                      >  2:5931139/1‑76 (MQ=255)
   aTTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCtttt                                                                      >  2:3851531/1‑76 (MQ=255)
    ttCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCCATCTCttttc                                                                     >  1:3587542/1‑75 (MQ=255)
     tCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCttttcc                                                                    >  1:289372/1‑74 (MQ=255)
         cccGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCAc                                                               <  2:2249926/76‑1 (MQ=255)
            gcgATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCt                                                            <  1:1358733/76‑1 (MQ=255)
            gcgATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCt                                                            >  1:2654795/1‑76 (MQ=255)
             cgATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCtg                                                           <  2:1732423/76‑1 (MQ=255)
              gATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCtgt                                                          >  1:2374821/1‑76 (MQ=255)
              gATGGGTGCGGCTGTCGAATAGTAAATAGGAGCAACCTTCCTCCAGGCGGCTTTCAATATCTTTTTCTCACCCcgt                                                          >  2:1759199/1‑76 (MQ=255)
              gATGGGTGCGGCTGGCGAATAGTAAAAAGTAGCAACCTTCCTCCAGGCGGCGTTCAATCTCTTTTCCTCCCCCtct                                                          >  2:1617433/1‑74 (MQ=255)
               aTGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCtgtg                                                         <  1:798909/76‑1 (MQ=255)
               aTGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCtgtg                                                         <  2:4930025/76‑1 (MQ=255)
               aTGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCtgtg                                                         >  1:3220776/1‑76 (MQ=255)
               aTGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCtgtg                                                         <  1:1895128/76‑1 (MQ=255)
                tGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTgg                                                        >  1:3146135/1‑76 (MQ=255)
                tGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTgg                                                        <  1:3461722/76‑1 (MQ=255)
                tGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTgg                                                        <  2:3837831/76‑1 (MQ=255)
                tGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTgg                                                        <  1:5108597/76‑1 (MQ=255)
                  ggTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGAc                                                      <  1:3024753/76‑1 (MQ=255)
                  ggTGCGGCTGGCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGGTTTCAATCTCTTTTTCTCACCCTGTGGAc                                                      >  2:3170165/1‑76 (MQ=255)
                   gTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGAcc                                                     >  1:7556/1‑76 (MQ=255)
                     gCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGGTTTCAATCTCTTTTCCTCACCCTGTGGACCtt                                                   >  2:4222280/1‑76 (MQ=255)
                     gCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCCCACCCTCTGGACCtt                                                   >  2:6293692/1‑76 (MQ=255)
                       ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCtttt                                                 <  2:2687786/76‑1 (MQ=255)
                       ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCtttt                                                 <  2:3989781/76‑1 (MQ=255)
                       ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCtttt                                                 <  1:4248213/76‑1 (MQ=255)
                       ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCtttt                                                 <  2:5345454/76‑1 (MQ=255)
                       ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCtttt                                                 <  2:4129407/76‑1 (MQ=255)
                         cTGTGGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCTCCCTGTGGACCTTTTAg                                               <  2:5212014/76‑1 (MQ=255)
                         cTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAg                                               <  1:2128451/76‑1 (MQ=255)
                         cTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAg                                               >  2:3427313/1‑76 (MQ=255)
                         cTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAg                                               >  1:5667686/1‑76 (MQ=255)
                          tGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGc                                              >  1:4081647/1‑76 (MQ=255)
                          tGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGc                                              >  1:3626761/1‑76 (MQ=255)
                           gTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCa                                             >  2:6005871/1‑76 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            >  1:2247681/1‑76 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            <  1:1002103/76‑1 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            <  1:2484618/76‑1 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            <  2:2642479/76‑1 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            >  2:5494036/1‑76 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            >  1:2263342/1‑76 (MQ=255)
                            tCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCAt                                            <  1:4868810/76‑1 (MQ=255)
                             cGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATc                                           <  2:3692412/76‑1 (MQ=255)
                              gAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCa                                          <  2:4368803/76‑1 (MQ=255)
                               aatagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAg                                         <  1:3846770/76‑1 (MQ=255)
                               aatagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAg                                         <  2:5284859/76‑1 (MQ=255)
                                atagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGa                                        >  2:4764021/1‑76 (MQ=255)
                                atagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGa                                        >  1:1023755/1‑76 (MQ=255)
                                atagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGa                                        <  1:5773911/76‑1 (MQ=255)
                                 tagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGAc                                       >  1:6172968/1‑76 (MQ=255)
                                 tagtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGAc                                       >  2:1981182/1‑76 (MQ=255)
                                  agtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGAcc                                      >  2:1854227/1‑76 (MQ=255)
                                  agtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGAcc                                      <  2:190891/76‑1 (MQ=255)
                                   gtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCt                                     <  1:2458558/76‑1 (MQ=255)
                                   gtaaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCt                                     <  1:976813/76‑1 (MQ=255)
                                    taaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTc                                    <  2:1584441/76‑1 (MQ=255)
                                    taaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTc                                    <  1:4708666/76‑1 (MQ=255)
                                    taaatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTc                                    >  1:4720068/1‑76 (MQ=255)
                                      aatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCa                                  <  1:2689308/76‑1 (MQ=255)
                                      aatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCa                                  <  1:5244566/76‑1 (MQ=255)
                                      aatagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCa                                  >  1:309896/1‑76 (MQ=255)
                                       atagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAg                                 >  1:4808562/1‑76 (MQ=255)
                                       atagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAg                                 >  2:2415077/1‑76 (MQ=255)
                                       atagtaGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAg                                 <  2:529799/76‑1 (MQ=255)
                                          gtagCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTg                              >  1:6025795/1‑76 (MQ=255)
                                           tagCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTgg                             >  2:4266258/1‑76 (MQ=255)
                                           tagCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTgg                             <  2:1046014/76‑1 (MQ=255)
                                           tagCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTgg                             >  2:1321006/1‑76 (MQ=255)
                                            agCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGCCCTCCAGCTGGc                            >  2:1638769/1‑76 (MQ=255)
                                            agCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGc                            <  1:4541039/76‑1 (MQ=255)
                                            agCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGc                            <  2:4426301/76‑1 (MQ=255)
                                             gCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCt                           <  2:3614848/76‑1 (MQ=255)
                                                aCCTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCtcgt                        <  1:1625628/76‑1 (MQ=255)
                                                 ccTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCtcgtc                       <  1:500401/76‑1 (MQ=255)
                                                 ccTTCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCtcgtc                       <  2:5887217/76‑1 (MQ=255)
                                                   ttCCTCCAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCgg                     >  2:4246270/1‑76 (MQ=255)
                                                     cctccAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGCCGGCa                   >  1:6243364/1‑76 (MQ=255)
                                                        ccAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGGaaat                >  2:4848862/1‑74 (MQ=255)
                                                        ccAGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAAtt                <  1:2795900/76‑1 (MQ=255)
                                                          aGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTc              >  2:5085903/1‑76 (MQ=255)
                                                          aGGCGGCTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTc              >  1:1104127/1‑76 (MQ=255)
                                                           ggcggcTTTCAATCGCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTcc             <  2:1487149/76‑1 (MQ=255)
                                                              ggcTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGc          <  2:1236251/76‑1 (MQ=255)
                                                                cTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGGCCTCCAGCTGGCTCGTCGGCAATTTCCGGCAg        <  1:5501117/76‑1 (MQ=255)
                                                                cTTTCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAg        <  1:659147/76‑1 (MQ=255)
                                                                 tttCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGa       >  1:4090563/1‑76 (MQ=255)
                                                                  ttCANTCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCCGAAAGTTTCGGCAGaa      >  1:3173083/1‑76 (MQ=255)
                                                                  ttCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGaa      >  2:4709763/1‑76 (MQ=255)
                                                                  ttCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGaa      <  1:5033342/76‑1 (MQ=255)
                                                                  ttCAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGaa      <  1:1473116/76‑1 (MQ=255)
                                                                    cAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAAcc    <  1:4322709/76‑1 (MQ=255)
                                                                    cAATCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCCTCGGCCATTTCCCGCAGaaac    >  1:2734839/1‑74 (MQ=255)
                                                                      aTCTCTTTTCCTCTCCCTGTGGCCCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGACCCCa  <  1:1182042/76‑1 (MQ=255)
                                                                      aTCTCTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCa  <  1:4625629/76‑1 (MQ=255)
                                                                           |                                                                       
CTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCA  >  REL606/79054‑79200

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