Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,313,864 Δ16,405 bp yjdPrdcB 18 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4313864 4330268 16405 38 [0] [0] 41 yjdP–rdcB yjdP,phnP,phnO,phnN,phnM,phnL,phnK,phnJ,phnI,phnH,phnG,phnF,phnD,phnC,yjdN,yjdM,rdcA,rdcB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 43138630 (0.000)38 (0.520) 37/490 0.5 100% intergenic (+64/‑5) rpiB/yjdP allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B/protein YjdP
?NC_000913 4330269 = 0 (0.000)intergenic (+31/‑233) rdcB/proP regulator of diguanylate cyclase RdcB/osmolyte:H(+) symporter ProP

GAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313614‑4314111
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
gAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGa                                                                                                                                                                                                                                                         <  1:517770‑M1/251‑2 (MQ=255)
   ccttcttacgctccgctttcgcggaagcataACGACATCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACGGGCAAAAATGATTTTGGATGCGTGGCTGGGCGCACAGTTCGAATGCGGTCGTCATCAACAACTCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTTAGATTCATCCACTACTTGTATGGATAAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaat                                                                                                                                                                                                                                                      <  2:321906‑M1/236‑5 (MQ=255)
            gcgcAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtg                                                                                                                                                                                                                                             <  1:499712‑M1/251‑14 (MQ=255)
                    ttCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttg                                                                                                                                                                                                                                       <  1:280872‑M1/249‑20 (MQ=255)
                      cgcgGCAGCATAACGCCACCAACGGGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaagccatgtgtgaagttgatc                                                                                                                                                                                                                                    <  2:433732‑M1/250‑23 (MQ=255)
                      cgcgGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatc                                                                                                                                                                                                                                    <  2:505201‑M1/250‑23 (MQ=255)
                                     acacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactg                                                                                                                                                                                                                    <  1:449507‑M1/251‑39 (MQ=255)
                                          aaCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagg                                                                                                                                                                                                               >  1:574540‑M1/1‑208 (MQ=255)
                                                      tttGGTTCACGAGTGGTTTGCCTCGAACTNGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggt                                                                                                                                                                                                     <  2:272225‑M1/249‑54 (MQ=255)
                                                          gTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgggtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaa                                                                                                                                                                                               <  1:463013‑M1/251‑60 (MQ=255)
                                                                           cTCGAACTCGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGGAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatcaatgtgagaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcacaattcaggcgt                                                                                                                                                                              >  1:405130‑M1/1‑175 (MQ=255)
                                                                             cGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaagtttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgt                                                                                                                                                                              <  1:320421‑M1/249‑77 (MQ=255)
                                                                                  tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCTTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagt                                                                                                                                                                                                                                         <  2:232878‑M1/185‑18 (MQ=255)
                                                                                  tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagt                                                                                                                                                                                                                                         >  1:232878‑M1/1‑168 (MQ=255)
                                                                                            gATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgt                                                                                                                                                              >  1:1310‑M1/1‑158 (MQ=255)
                                                                                                      gCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattcctta                                                                                                                                                    >  1:478461‑M1/1‑148 (MQ=255)
                                                                                                        gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaacc                                                                                                                                                 <  2:1310‑M1/251‑106 (MQ=255)
                                                                                                                          cGAAGTCGGTCGTCAGCAACATCGCGTGGAGGCGATTACTGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTTGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactgttagggtaaaaaggtcattaacggcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacc                                                                                                                               <  1:414530‑M1/251‑124 (MQ=255)
                                                                                                                                gggtCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaattaaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctac                                                                                                                         <  2:131540‑M1/250‑130 (MQ=255)
                                                                                                                                         tcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttac                                                                                                                >  1:475784‑M1/1‑113 (MQ=255)
                                                                                                                                             caacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagag                                                                                                            <  1:404955‑M1/251‑143 (MQ=255)
                                                                                                                                                      gAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcc                                                                                                   >  1:12882‑M1/1‑100 (MQ=255)
                                                                                                                                                                  gCAATAGAGCAGCGGAGAAATTGAGATTCATCCTCTACTTTCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaatcaaacccgctacgcttgttacagagattgcatcctgcaattcccgc                                                                                       <  2:169528‑M1/251‑164 (MQ=255)
                                                                                                                                                                     aTAGAGCAGCGGAGAAATTGAGATTCCTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgggtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcc                                                                                    <  2:26423‑M1/251‑167 (MQ=255)
                                                                                                                                                                                    aaTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtc                                                                     <  1:399259‑M1/251‑182 (MQ=255)
                                                                                                                                                                                       tGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgc                                                                   <  1:334350‑M1/250‑184 (MQ=255)
                                                                                                                                                                                        gagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcg                                                                  <  2:9805‑M1/250‑185 (MQ=255)
                                                                                                                                                                                                         tGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtt                                                 <  2:201716‑M1/250‑202 (MQ=255)
                                                                                                                                                                                                            aTGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcg                                             >  1:544328‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                              ggATGAGTAAGGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcgg                                            <  2:469520‑M1/250‑207 (MQ=255)
                                                                                                                                                                                                                     tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggc                                    >  1:403567‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                       aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggcc                                   >  1:42858‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                        tGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcaaaaatttaaaaactggtagggtaaaaaggtcattaactgcccaattaaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttatggcgtttgcggaaatgggcca                                  >  1:295612‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                            tAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacacgggtagggtaaaaaggtcattaactgcccaatacaggcgacaactggtttgattgtacattccttaaccggcgggtgtaagcaaacccgctacgcttgttacagagactgcatcctgcaattcccgctccccttttgcggccggcgcgctgatttttctggagtttgcggaaatgtgccaactcc                             >  2:41390‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                                                  ggATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctngcgtttgcggaactgggccaactctgcgagg                       >  2:291214‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                                                          tttCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctat               >  1:202062‑M1/1‑16 (MQ=255)
                                                                                                                                                                                                                                             cgttgtgaggtTGaaatcaatgtgggaagttgatcacaaatgtaaacactgttagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgatacgcttgttacagagatttcatcctgcaatgcccgctccccttttgcggccgtagcgcggatttttctggcgtttgcggaaatgggccatctcttcgaggaaagctatgc             <  2:16398‑M1/245‑238 (MQ=255)
                                                                                                                                                                                                                                                        tGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggcaagctatgctgaaaaggaaa  >  2:331610‑M1/1‑2 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
GAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313614‑4314111

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