Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 3,843,787 Δ54 bp coding (251‑304/1404 nt) dnaA ← chromosomal replication initiation protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3843787 3843840 54 134 [6] [1] 131 dnaA chromosomal replication initiation protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 38437866 (0.040)104 (0.790) 58/132 0.8 97.1% coding (305/1404 nt) dnaA chromosomal replication initiation protein
?REL606 3843841 = 1 (0.010)coding (250/1404 nt) dnaA chromosomal replication initiation protein

Only 100 of 131 total aligned reads displayed.

CACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGAT  >  REL606/3843766‑3843916
                                                                           |                                                                           
tcggttctgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgt                                                                             >  1:399469‑M2/76‑76 (MQ=255)
  ggttctgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCa                                                                           <  1:45225‑M2/3‑1 (MQ=255)
  ggttctgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCa                                                                           >  2:3266834‑M2/74‑76 (MQ=255)
  ggttctgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCa                                                                           >  1:1881616‑M2/74‑76 (MQ=255)
  ggttctgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcggttgccTCa                                                                           >  2:3148762‑M2/74‑76 (MQ=255)
     tctgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCGCCg                                                                        >  1:873730‑M2/71‑76 (MQ=255)
       tgccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgGCACcggg                                                                      >  2:2195819‑M2/69‑75 (MQ=255)
         ccggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCCGttt                                                                    >  2:3868152‑M2/67‑76 (MQ=255)
          cggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcggttgcgTCACCGGGTTg                                                                   >  2:3099508‑M2/66‑76 (MQ=255)
          cggggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcggttgcgGCACCGGGTTg                                                                   >  1:680586‑M2/66‑76 (MQ=255)
             ggccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgGCACCGGGTTGGTg                                                                >  2:2285561‑M2/63‑76 (MQ=255)
               ccgggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGcc                                                              <  2:1812929‑M2/16‑1 (MQ=255)
                 gggacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgttttcgTCACCCGTTTGGTGCCGa                                                            >  2:485994‑M2/59‑76 (MQ=255)
                   gacgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACt                                                          <  2:4674863‑M2/20‑1 (MQ=255)
                    acgttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACtt                                                         >  1:587575‑M2/56‑76 (MQ=255)
                       ttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCaa                                                      <  1:3574321‑M2/24‑1 (MQ=255)
                       ttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCaa                                                      <  1:1229682‑M2/24‑1 (MQ=255)
                       ttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGGGCCGACttcca                                                      >  2:3839219‑M2/53‑74 (MQ=255)
                       ttatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGGGCCGACTTCaa                                                      >  1:1295502‑M2/53‑76 (MQ=255)
                        tatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTTGGGCCGACTTCaaa                                                     >  2:3378772‑M2/52‑76 (MQ=255)
                        tatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCaaa                                                     <  2:4970194‑M2/25‑1 (MQ=255)
                        tatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCaaa                                                     <  2:2270455‑M2/25‑1 (MQ=255)
                        tatcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCCGTTTGGTGCCGACTTCaaa                                                     >  2:3157493‑M2/52‑76 (MQ=255)
                         atcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCaaaa                                                    >  1:4325581‑M2/51‑76 (MQ=255)
                         atcccaacctgagcgcgtagaaggcgcagcacgttgcggctgcggttgcgTCACCGGTTTGGGGGCGACTTCaaaa                                                    >  2:3570353‑M2/51‑76 (MQ=255)
                           cccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTtcaaaaac                                                  >  1:2303152‑M2/49‑74 (MQ=255)
                            ccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGc                                                 <  2:421080‑M2/29‑1 (MQ=255)
                            ccaacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGc                                                 <  2:1438538‑M2/29‑1 (MQ=255)
                             caacctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCa                                                <  1:3833904‑M2/30‑1 (MQ=255)
                               acctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGc                                              <  1:2233248‑M2/32‑1 (MQ=255)
                               acctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGc                                              <  1:4153552‑M2/32‑1 (MQ=255)
                                 ctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCtg                                            <  2:2275454‑M2/34‑1 (MQ=255)
                                 ctgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAAAGCCGCtg                                            >  1:990050‑M2/43‑76 (MQ=255)
                                 ctgagcgagtagaaggcgcagcgcgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCtg                                            <  2:4280404‑M2/34‑1 (MQ=255)
                                  tgagcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCtgt                                           <  2:4434285‑M2/35‑1 (MQ=255)
                                  tgagcgcgtagaaggcgcagcacgttgcggctgcctttgcgTCACCGGGTTGGTGGCGGCTTCAAAAAGCAgcttt                                           >  1:2895301‑M2/42‑74 (MQ=255)
                                     gcgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGGTTGGTGCCGACTTCAAAACGCAGCTGTggg                                        >  2:1528057‑M2/39‑76 (MQ=255)
                                      cgcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCCGTTTGGTGCCGGCTTCAAAACGCAGCTgggggg                                       >  2:1324436‑M2/38‑76 (MQ=255)
                                       gcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCCGCTGTGGGGc                                      >  1:611490‑M2/37‑76 (MQ=255)
                                       gcgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGGGCCGACTTCAAAACGCAGCTGTGGGGc                                      >  1:4935781‑M2/37‑76 (MQ=255)
                                        cgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGGCTTCAAAACGCCGCTgtggggga                                     >  2:4937349‑M2/36‑74 (MQ=255)
                                        cgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTgtggggga                                     >  2:4475887‑M2/36‑74 (MQ=255)
                                        cgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCa                                     <  2:120504‑M2/41‑1 (MQ=255)
                                        cgtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCa                                     <  2:3652960‑M2/41‑1 (MQ=255)
                                         gtagaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAt                                    >  1:2387975‑M2/35‑76 (MQ=255)
                                           agaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTTGGGCCGACTTCAAAACGCAGCTGTGGGGCAacc                                  >  2:3143023‑M2/33‑76 (MQ=255)
                                           agaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtcc                                  >  2:2466533‑M2/33‑76 (MQ=255)
                                           agaaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtcc                                  >  1:3931131‑M2/33‑76 (MQ=255)
                                           agaaggcgcagcacgttgcggctgcgtttgcgTCACCGGGTTGGTGCAGCCTTCAAAAAGCCACTGTCGGGCAtcc                                  >  2:4426222‑M2/33‑76 (MQ=255)
                                             aaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgc                                >  2:2360089‑M2/31‑76 (MQ=255)
                                             aaggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGGGCCGACTTCAAAACGCAGCTGTGGGGCAtccgc                                >  2:3080267‑M2/31‑76 (MQ=255)
                                              aggcgcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgct                               <  1:3917856‑M2/47‑1 (MQ=255)
                                                  gcagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgctccgc                           <  2:611013‑M2/51‑1 (MQ=255)
                                                   cagcgcgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCa                          <  2:3186541‑M2/52‑1 (MQ=255)
                                                   cagcacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCa                          >  2:1970344‑M2/25‑76 (MQ=255)
                                                      cccgtggcgtctgcgtttgcgTCACCGGTTTTGTGCCGACTTCAAAACGCAGCTGTGGGTCATCCGCTCCACAGaa                       <  1:3339471‑M2/55‑1 (MQ=255)
                                                      cacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGaa                       >  2:2570851‑M2/22‑76 (MQ=255)
                                                      cacgttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGaa                       <  2:1231616‑M2/55‑1 (MQ=255)
                                                         gttgcggctgcgtttgcgTCACCGGTTTGGTGCCGCCTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACt                    <  1:2982340‑M2/58‑1 (MQ=255)
                                                         gttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACt                    <  2:1823717‑M2/58‑1 (MQ=255)
                                                         gttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCCGTGGGGGATCCGCTCCCCagaaaac                    >  2:2786208‑M2/19‑74 (MQ=255)
                                                          ttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGACACTg                   >  2:1234702‑M2/18‑76 (MQ=255)
                                                          ttgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTg                   >  2:4159024‑M2/18‑76 (MQ=255)
                                                          ttgcggctgcgtttgcgTCACCGGTTTGGGGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTg                   <  2:2570422‑M2/59‑1 (MQ=255)
                                                           tgcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCCGCTGTGGGGCATCCGCTCCGCAGAAACTgg                  >  2:4968005‑M2/17‑76 (MQ=255)
                                                            gcggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGt                 <  2:784260‑M2/61‑1 (MQ=255)
                                                            gcggctgcgtttgcgTCACCCGTTTGGTGCCGACTTCAAAACGCAGCCGTGGGGCATCCGCTCCGCAGAAACTGGt                 >  2:3377109‑M2/16‑76 (MQ=255)
                                                             gggctgcgtttgcgTCACGGGTTTGGTGCCTCCTTCAAAACGCAGCTGTGGGGCATCCTCTCCGCAGAAACTGGtt                <  2:3292163‑M2/62‑1 (MQ=255)
                                                             cggctgcgtttgcgTCACTGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGtt                <  2:1558114‑M2/62‑1 (MQ=255)
                                                             cggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGtt                >  2:434401‑M2/15‑76 (MQ=255)
                                                             cggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGtt                <  2:626625‑M2/62‑1 (MQ=255)
                                                             cggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGtt                <  1:10445‑M2/62‑1 (MQ=255)
                                                             cggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGtt                <  1:4878536‑M2/62‑1 (MQ=255)
                                                             cggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGtt                <  2:4142648‑M2/62‑1 (MQ=255)
                                                              gtctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTa               <  2:161860‑M2/63‑1 (MQ=255)
                                                              ggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTa               <  1:2035264‑M2/63‑1 (MQ=255)
                                                              ggctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTa               <  1:875905‑M2/63‑1 (MQ=255)
                                                              ggcggCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTa               <  1:4743528/73‑1 (MQ=255)
                                                               gctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAg              <  1:4555241‑M2/64‑1 (MQ=255)
                                                                ctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGc             <  2:4897598‑M2/65‑1 (MQ=255)
                                                                ctgcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGc             >  1:2955919‑M2/12‑76 (MQ=255)
                                                                  gcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAg           >  2:4463340‑M2/10‑76 (MQ=255)
                                                                  gcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAg           >  2:76724‑M2/10‑76 (MQ=255)
                                                                  gcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAg           >  2:2145005‑M2/10‑76 (MQ=255)
                                                                  gcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAg           >  2:3515254‑M2/10‑76 (MQ=255)
                                                                  gcgtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCAGCAGAAACTGGTTAGCAg           >  2:4868550‑M2/10‑76 (MQ=255)
                                                                    gtttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTc         >  2:3566674‑M2/8‑76 (MQ=255)
                                                                    gtttgcgGCACCGGTTTGGTGCCGCCTGCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTc         <  1:3744432‑M2/69‑1 (MQ=255)
                                                                    gttggcgTCACCGGTTTGGTCCCAACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTc         <  1:1080502‑M2/69‑1 (MQ=255)
                                                                     tttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTcc        >  1:3499091‑M2/7‑76 (MQ=255)
                                                                     tttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTcc        <  1:2012769‑M2/70‑1 (MQ=255)
                                                                      ttgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCa       <  1:1489399‑M2/71‑1 (MQ=255)
                                                                       tgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAt      >  2:4813021‑M2/5‑76 (MQ=255)
                                                                       tgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAt      >  2:4989400‑M2/5‑76 (MQ=255)
                                                                       tgcgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAt      >  2:67153‑M2/5‑76 (MQ=255)
                                                                         cgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTg    <  1:45012‑M2/74‑1 (MQ=255)
                                                                         cgTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTg    <  1:2183730‑M2/74‑1 (MQ=255)
                                                                          gTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGa   <  1:3891237‑M2/75‑1 (MQ=255)
                                                                           tCACCGTTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGat  <  2:4043058/76‑1 (MQ=255)
                                                                           tCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGat  <  1:1579998/76‑1 (MQ=255)
                                                                           |                                                                           
CACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGAT  >  REL606/3843766‑3843916

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