Tool: Sorting bed files with bash sort
I've been implementing this command for the past few days now and I wished I did it earlier.
It functions similar to sortBed in bedtools but uses the bash sort which can be used in pipes
Find your .bashrc file in your home directory
$ cd $HOME
$ vi .bashrc # vi ~/.bashrc should work fine from any directory
Add this to .bashrc
alias sortbed="sort -k1,1 -k2,2g "
Save and source your .bashrc to get it to work
$ source ~/.bashrc
Example of use
$ intersectBed -a in.bed -b /segDup_unmappable.bed -wao | sortbed |uniq | cut -f 1,2,3,8 >in_segDup_unmappable_overlap.txt
This intersects a bed file of chr, start, end to a list of segmental duplications and unmappable regions in hg19. It also pipes to bash commands to only remove the in positions and the number of base pairs overlapping it.
sortbed is used to sort the output and uniq is applied to return only unique lines. You can treat sortbed like sort. Just a nice shortcut I thought others might like.
• 6,176 views
•
link
updated
by
Ram
•
written
by
QVINTVS_FABIVS_MAXIMVS
•
written
by
QVINTVS_FABIVS_MAXIMVS
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
rMATS-turbo build_rmats --conda issue: Run 'conda init' before 'conda activate'
written by Angelos •
-
HISAT2 error: Use of uninitialized value $l
written by sanchezgil.juanjose •
-
Tutorial: Creating workflows with snakemake and conda
written by finswimmer
-
PBS script not able to locate the tool
written by MAPK
-
Tutorial: A Wet-Lab Guide to 'what I wish I had known' for Computational Biology
written by Sinji
-
Tool: Sorting bed files with bash sort
written by QVINTVS_FABIVS_MAXIMVS
Like I said, bedtools sort gives the UNIX command to sort bed files; I don't see any advantage in your approach.
and set
LC_ALL=Cto make things faster.BEDOPS sort-bed works faster at sorting BED files than GNU sort, and you can pipe data in and out via standard UNIX streams.
Unlike other tools, it also handles arbitrary numbers of columns and can be assigned a chunk of memory, to sort very large BED files that will not otherwise fit into system memory.
Add semantic version sort to the first key for the chromosomes to be sorted correctly.