view convertMAF.sh @ 100:666c07b60cd2 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/jbrowse2 commit a1537aea75fc902d0e38c0b7c698830a939648b1-dirty
author fubar
date Fri, 21 Jun 2024 23:49:10 +0000
parents e7a6f7a7148d
children
line wrap: on
line source

#!/usr/bin/env bash
# https://github.com/cmdcolin/jbrowse-plugin-mafviewer/blob/master/bin/convert.sh
# maf2bed modified to work right as a python script by ross lazarus in desperation
#  MAF file must contain the species name and chromosome name
#  e.g. hg38.chr1 in the sequence identifiers with hg38 passed in as $2
python $3/maf2bed.py $2 < $1 | sort -k1,1 -k2,2n > $4.sorted.bed
bgzip -c $4.sorted.bed > $4.sorted.bed.gz
tabix -p bed $4.sorted.bed.gz