Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 837,880 (G)5→6 coding (226/2226 nt) ybiO ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606837,8751.G98.0% 395.5 / 0.5 98coding (231/2226 nt)ybiOpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (55/41);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 122 total aligned reads displayed.

GACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑GGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACAT  >  REL606/837800‑837945
                                                                            |                                                                       
gACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTcagca                                                                          <  1:214837/76‑1 (MQ=255)
gACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACAGGTTcagca                                                                          >  2:5876025/1‑76 (MQ=255)
gACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGGACCGGTTcagca                                                                          >  1:353130/1‑76 (MQ=255)
gACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTACTTTTGGTACCGGTTcagca                                                                          >  2:5056290/1‑76 (MQ=255)
 aCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTcagca‑‑g                                                                       <  1:4852121/76‑1 (MQ=255)
 aCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTcagca‑‑g                                                                       <  1:1439943/76‑1 (MQ=255)
 aCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTcagca‑‑g                                                                       <  2:2432875/76‑1 (MQ=255)
   ttCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑ggg                                                                     <  2:1845848/76‑1 (MQ=255)
   ttCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑ggg                                                                     >  1:3206423/1‑76 (MQ=255)
    tCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCC‑‑gggg                                                                    >  1:2791974/1‑76 (MQ=255)
      ggTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑gggggg                                                                  >  2:1571671/1‑75 (MQ=255)
      ggTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑gggggg                                                                  <  2:2477009/76‑2 (MQ=255)
      ggTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑gggggg                                                                  <  2:4708010/76‑2 (MQ=255)
      ggTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACAGGTTCAGCA‑‑gggggg                                                                  >  2:5979052/1‑75 (MQ=255)
           ccTGTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAg                                                              <  2:1176953/76‑1 (MQ=255)
            cTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTATTGTACCGGTTCAGCAGGGGGGCGTAGc                                                             >  2:6162224/1‑76 (MQ=255)
             ttttGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCg                                                            <  1:4377217/76‑1 (MQ=255)
              tttGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCgg                                                           >  1:5735934/1‑76 (MQ=255)
               ttGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGc                                                          >  1:1975374/1‑76 (MQ=255)
                 gcagcaCGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCaa                                                        >  2:824152/1‑76 (MQ=255)
                  cagcaCGATTTGCTCTTCGACCAGCGTCGGCGGCAATATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAAc                                                       >  2:1804794/1‑76 (MQ=255)
                    gcaCGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAAccg                                                     >  2:4288764/1‑74 (MQ=255)
                     caCGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGt                                                    <  2:1107661/76‑1 (MQ=255)
                      aCGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCCACGGTg                                                   >  2:5813507/1‑76 (MQ=255)
                      aCGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACgggg                                                   >  1:2856103/1‑74 (MQ=255)
                       ctgTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTgc                                                  >  1:6060846/3‑76 (MQ=255)
                       cGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTgc                                                  <  2:4171265/76‑1 (MQ=255)
                        ggTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGGCGTAGCGGCAACGGGgc                                                  >  2:4437228/1‑76 (MQ=255)
                        ggTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTgcg                                                 >  2:1682215/1‑76 (MQ=255)
                        ggTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGGgcg                                                 >  2:1998438/1‑76 (MQ=255)
                        ggTTTGCGCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTgca                                                 <  2:6110164/76‑2 (MQ=255)
                         gTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTgcgc                                                >  1:3603087/1‑76 (MQ=255)
                         gTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCCGCAAAGGTgcgc                                                >  2:5792071/1‑76 (MQ=255)
                           ttGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGGGCGCaa                                              >  1:4618524/1‑76 (MQ=255)
                           ttGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAAAGGTGCGCaa                                              >  2:721576/1‑76 (MQ=255)
                            tGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAAc                                             >  2:6046170/1‑76 (MQ=255)
                            tGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAAACGTGCGCCAc                                             >  1:1147064/1‑76 (MQ=255)
                             gCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACt                                            <  1:121911/76‑1 (MQ=255)
                             gCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACt                                            >  1:940355/1‑76 (MQ=255)
                             gCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑GGGGGCGTAGCGGCAACGGTGCGCAACTg                                           >  1:4919590/1‑76 (MQ=255)
                                ctTCGGCCAGCGTCGGCGGTACTATTTTTGGTACCGGTTAAGCAGGGGGGCGTAGCCGCAACAGGGCGCAACTGGt                                         >  2:341247/1‑76 (MQ=255)
                                ctTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGGGCGCAACTggc                                         >  2:2977978/1‑75 (MQ=255)
                                 ttCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTc                                        >  1:6147315/1‑76 (MQ=255)
                                 ttCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTc                                        <  1:5884520/76‑1 (MQ=255)
                                  tCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCg                                       >  2:1487727/1‑76 (MQ=255)
                                   cGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGa                                      <  1:900554/76‑1 (MQ=255)
                                    gACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGAt                                     <  1:5961818/76‑1 (MQ=255)
                                     aCCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGCTc                                    >  1:3555847/1‑76 (MQ=255)
                                      ccAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCa                                   <  1:3822869/76‑1 (MQ=255)
                                      ccAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCa                                   <  1:136506/76‑1 (MQ=255)
                                      ccAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCa                                   <  2:3159228/76‑1 (MQ=255)
                                      ccAGCGTCGGCGGCACTATGTTTGGTACCGATTCAGCAGGGGGGCTTAGCGGCAACGGTGCGCAACTGGTCGATCa                                   <  2:307813/76‑1 (MQ=255)
                                       cAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCaa                                  <  2:4948947/76‑1 (MQ=255)
                                       cAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGGGCGCAACTGGTCGATCaa                                  >  2:2004017/1‑76 (MQ=255)
                                        aGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGGGCGCCACTGGTCGATCAAc                                 >  1:778020/1‑76 (MQ=255)
                                        aGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGGGCGCAACTGGGCGATCAAc                                 >  2:5991353/1‑76 (MQ=255)
                                         gCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAAct                                <  1:824152/76‑1 (MQ=255)
                                           gTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAActct                              <  1:1737525/76‑1 (MQ=255)
                                           gTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAActct                              <  1:480213/76‑1 (MQ=255)
                                           gTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGGCGATCAActct                              >  1:3944430/1‑76 (MQ=255)
                                            tCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCtt                             <  2:2592257/76‑1 (MQ=255)
                                             cggcggGACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGTGGCAACGGTGCGAAACTGGTCGATCAACTCttt                            <  2:4552995/76‑1 (MQ=255)
                                             cggcggCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCttt                            >  2:1340345/1‑76 (MQ=255)
                                             cggcggCACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCttt                            <  1:5442666/76‑1 (MQ=255)
                                               gcggcACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTAc                          <  1:5429046/76‑1 (MQ=255)
                                                cggcACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTAcg                         >  1:3319770/1‑76 (MQ=255)
                                                  gcccTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTAcgcg                       <  1:2452937/73‑1 (MQ=255)
                                                  gcACTATTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTAcgcg                       <  1:1117861/76‑1 (MQ=255)
                                                       aTTTTTGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGgtg                  >  2:3833486/1‑76 (MQ=255)
                                                         ttttGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTc                >  1:3995814/1‑76 (MQ=255)
                                                          tttGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGTTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCa               >  1:4628474/1‑76 (MQ=255)
                                                          tttGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCa               >  1:2614032/1‑76 (MQ=255)
                                                           ttGGTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCAt              >  2:5593638/1‑76 (MQ=255)
                                                             ggTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTa            >  2:2783431/1‑76 (MQ=255)
                                                             ggTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTa            >  1:5237959/1‑76 (MQ=255)
                                                             ggTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTa            <  1:1163379/76‑1 (MQ=255)
                                                              gTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTAt           <  1:5417227/76‑1 (MQ=255)
                                                              gTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTAt           <  1:5249673/76‑1 (MQ=255)
                                                              gTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTAt           >  1:3097448/1‑76 (MQ=255)
                                                              gTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTAt           >  2:1140155/1‑76 (MQ=255)
                                                              gTACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGGGTCATTAt           >  2:5816001/1‑76 (MQ=255)
                                                               tACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCCTTATc          >  1:5538514/1‑76 (MQ=255)
                                                               tACCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATc          >  1:504187/1‑76 (MQ=255)
                                                               tACCGGTTCAGCAGGGGGGCGAAGCGGCAACGGAGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATc          >  2:4718902/1‑76 (MQ=255)
                                                                aCCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATcc         <  1:2344652/76‑1 (MQ=255)
                                                                aCCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATcc         <  1:1998438/76‑1 (MQ=255)
                                                                aCCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATcc         <  2:3860806/76‑1 (MQ=255)
                                                                aCCGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATcc         >  2:513198/1‑76 (MQ=255)
                                                                 ccGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCa        <  1:4134981/76‑1 (MQ=255)
                                                                  cGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAg       <  1:3978389/76‑1 (MQ=255)
                                                                  cGGTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAg       >  2:2276208/1‑76 (MQ=255)
                                                                   ggTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGc      <  1:2394313/76‑1 (MQ=255)
                                                                    gTTCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCAAGca     >  2:4576630/1‑76 (MQ=255)
                                                                     ttCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGcac    <  2:5953237/76‑1 (MQ=255)
                                                                     ttCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGcac    >  2:5894949/1‑76 (MQ=255)
                                                                     ttCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGcac    <  1:299675/76‑1 (MQ=255)
                                                                     ttCAGCAGGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGcac    >  2:1327046/1‑76 (MQ=255)
                                                                       cagcagGGGGGCGTAGCGGCAACGGTGCGCCACTGGTCGAACAACCTTTTACGCGCGGTGGCATTATACAGAACAt  >  2:343867/1‑76 (MQ=255)
                                                                       cagcagGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACAt  <  2:205614/76‑1 (MQ=255)
                                                                       cagcagGGGGGCGTAGCAGCACCGGTGCGCAACTGGGCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACAt  <  2:5649523/76‑1 (MQ=255)
                                                                            |                                                                       
GACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCA‑‑GGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACAT  >  REL606/837800‑837945

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