Assembly conversion of large BED files
2
0
Entering edit mode
4.0 years ago
bioinfo456 ▴ 150

Hi all, I have a set of BED files containing certain peak locations obtained using NCBI Build GRCh37/UCSC Build hg19 genome assembly. I want the same in GRCh38 assembly. Ensembl's assembly convertor has an upload restriction of 50MB. Any alternate solutions please? Thanks.

Assembly ensembl • 3.2k views
ADD COMMENT
1
Entering edit mode
4.0 years ago

UCSC liftOver: https://genome.ucsc.edu/cgi-bin/hgLiftOver

Run on the command-line to get around web file size restrictions.

ADD COMMENT
1
Entering edit mode
4.0 years ago
bernatgel ★ 3.3k

You can use liftOver in your computer with no restrictions.

We have used both the UCSC liftOver executable and the R Bioconductor liftOver function in the rtracklayer package and both have worked with no problems. In any case you will need a chain file for your specific conversion going from your your original assembly to your target assembly.

ADD COMMENT
0
Entering edit mode

Thanks. Just to cross confirm, the chain file needed for the BED file conversion in my case is to be downloaded under Human genome -> (GRCh37/hg19) -> LiftOver files -> hg19ToHg38.over.chain.gz. Right?

ADD REPLY
1
Entering edit mode

Yes, that's the one you need

ADD REPLY

Login before adding your answer.

Traffic: 2568 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6