How to sort chromosome name in bed file

Hello Friends!

I have one stupid question. I have a bed file

chr1 955541 955763 AGRN.chr1.955552.955753
chr1 957569 957852 AGRN.chr1.957580.95784

in my bed file my chromosome name are not sorted I want output in

chr1 chr2 ..... chr22 chrX chrY chrM

Can anyone help me out

Thanks!

bed

• 4,590 views • link

updated 3.4 years ago by Ram

• written 10.1 years ago by raya.girish •

use the sortBed

• 1 views • link

written 10.1 years ago by Prasad

can you tell me how should i sort using sortbed i mean commands

• 1 views • link

written 10.1 years ago by raya.girish •

i used bedtool/sort -chrThenSizeA -i my.bed > new.bed

it not working

• 1 views • link

written 10.1 years ago by raya.girish •

Please see the answer below and provide feedback.

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. Upvote|Bookmark|Accept

• 0 views • link

written 8.6 years ago by Ram

1 answer

sortBed -i inputfile >outfile

By default it sorts by chromosome and then by start position in ascending order. It has many other options, kindly go through the manual

• 1 views • link

updated 10.1 years ago by GenoMax

• written 10.1 years ago by Prasad

@Prasad - It would be better if you change it answer from comment and format the code, so that OP can accept the answer and the question will not remain as unanswered, and also helpful for future googlers.

• 1 views • link

written 10.1 years ago by venu

@Venu - Sure i will keep that in mind. @genomax - thanks for the format

• 1 views • link

written 10.1 years ago by Prasad

More posts like this

How to get chrM chr1 ...chr22 chrX chrY order in my bam file
written 6.8 years ago by Ada •

Hi I went BWA alignment and sorted my bam file when I was doing my WES analysis to find SNP and Indel, however, my bam …
How can I get this sorted by ascending order? (when the number and character is together)
written 7.1 years ago by woojoy14 •

samtools view NA12878_dedup.bam | cut -f 3 | sort -u > chr.list8 > chr1 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr2 …
Reordering Bam files based on chromosome number, when the chromosome format is chr1..etc
written 8.0 years ago by jparker4 •

When trying to call variants I get the error message: ##### ERROR MESSAGE: Lexicographically sorted human genome sequence detected in reads. Please see https://software.broadinstitute.org/gatk/documentation/article?id=1328for more …
Removing a chromosome from a fasta file
written 8.3 years ago by t86dan •

Hello, I have a fasta file for human Genome GRCh37 Reference Assembly, but for some reason chromosome 20 is repeated and I would like to …
coverageBed with -g option ?
written 10.9 years ago by bioguy24

I am running a coverageBed and getting the below error. ``` cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -d -sorted -a /home/cmccabe/Desktop/NGS/bed/bedtools/xgen_targets_sorted.bed -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ERROR: Database file /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam …
coverage bed error but produces output file
written 10.9 years ago by bioguy24

The below command produces a file (674 MB). ``` coverageBed \ -d \ -sorted \ -a /home/cmccabe/Desktop/NGS/bed/sorted_unix_5column_xgen_targets.bed \ -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam \ > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ``` Output …
add strand info to a bed file
written 11.1 years ago by bioguy24

Is there a script/tool that can add strand (+/-) information to the fourth column of a bed file? I have a bed file (first few …
How To Sort Bed Format File
written 13.5 years ago by lyz10302012

Is there any tools or commands to sort bed format file? The correct order is chr1, chr2, chr3, chr4,..., chr22, chrX, chrY.
How To Reorder Bam File Contigs For The Gatk Pipeline
written 13.7 years ago by win

i have a sorted & read group added bam file and when i try to run the following GATK command on it <pre><code>java -jar GenomeAnalysisTK.jar …
Cuffdiff Not Seeming To Work With Ucsc Gtf
written 15.6 years ago by Sam •

Hello, I've been working with a public dataset from the SRA. I've gone through conventional RNAseq workflow of alignment and sorting the SAM files using …