Mercurial > repos > fubar > jbrowse2

changeset 24:fb6cc7bc24df draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/usegalaxy-eu/temporary-tools/tree/master/jbrowse2 commit cde4b8a574ded34a0ff8df3ecafc1a057787dcfb-dirty
author fubar
date Sat, 03 Feb 2024 11:49:56 +0000
parents 39b717d934a8
children 172acf9c7a53
files cool_in_jb2.jpg dm4_in_jb2.png jbrowse2.py readme.rst
diffstat 4 files changed, 9 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
Binary file cool_in_jb2.jpg has changed
Binary file dm4_in_jb2.png has changed
--- a/jbrowse2.py	Sat Feb 03 10:17:27 2024 +0000
+++ b/jbrowse2.py	Sat Feb 03 11:49:56 2024 +0000
@@ -461,7 +461,7 @@
                 self.genome_firstcontig = None
                 fl = open(fapath, "r").readline().strip().split(">", 1)
                 if len(fl) > 1:
-                    self.genome_firstcontig = fl[1].strip()
+                    self.genome_firstcontig = fl[1].split().strip()
         if self.config_json.get("assemblies", None):
             self.config_json["assemblies"] += assemblies
         else:
@@ -1148,7 +1148,6 @@
                     outputTrackConfig,
                 )
             elif dataset_ext in ("cool", "mcool", "scool"):
-                hictempd = tempfile.mkdtemp()
                 hic_path = os.path.join(
                     self.outdir, "%s_%d_%s.hic" % (track_human_label, i, dataset_ext)
                 )
@@ -1159,8 +1158,6 @@
                         "-f",
                         "--output-fmt",
                         "hic",
-                        "--tmpdir",
-                        hictempd,
                         dataset_path,
                         hic_path,
                     ]
@@ -1169,7 +1166,6 @@
                     hic_path,
                     outputTrackConfig,
                 )
-                shutil.rmtree(hictempd)
             elif dataset_ext in ("bed",):
                 self.add_bed(
                     dataset_path,
@@ -1283,9 +1279,9 @@
                 refName = loc_match.group(1)
                 drdict["refName"] = refName
                 if loc_match.group(2) > "":
-                    drdict["start"] = int(loc_match.group(2))
+                    drdict["start"] = int(loc_match.group(2).replace(',',''))
                 if loc_match.group(3) > "":
-                    drdict["end"] = int(loc_match.group(3))
+                    drdict["end"] = int(loc_match.group(3).replace(',',''))
             else:
                 logging.info(
                     "@@@ regexp could not match contig:start..end in the supplied location %s - please fix"
--- a/readme.rst	Sat Feb 03 10:17:27 2024 +0000
+++ b/readme.rst	Sat Feb 03 11:49:56 2024 +0000
@@ -1,6 +1,12 @@
 JBrowse2 in Galaxy
 ==================
 
+Added Feb 3: cool/mcool -> hic 
+
+.. image:: dm4_in_jb2.png
+
+Other tracks
+
 .. image:: jbrowse8.png
 
 JBrowse2 is a fast, embeddable genome browser built completely with