Why is there a big difference between the outputs of `geneBody_coverage` of `RSeQC` when using different bed files?
0
3
Entering edit mode
18 months ago
Dan ▴ 180

Hi

I run geneBody_coverage.py of RSeQC (https://rseqc.sourceforge.net/#contact) using different bed files and the same batch of bam files:

geneBody_coverage.py -r long_short_genes/changeformat/hg19.RefSeq.union.bed \
-i split_bam/short_genes/bams/ \
-o geneBody_coverage/short_genes/RefSeq

The output is very different.

What is difference between union.bed files and other bed files which I downloaded from the following links:

https://sourceforge.net/projects/rseqc/files/BED/Human_Homo_sapiens_merge_transcripts/

https://sourceforge.net/projects/rseqc/files/BED/Human_Homo_sapiens/

Which bed file should I use?

Thanks

RSeQC • 912 views
ADD COMMENT
1
Entering edit mode

If you scroll down the pages you linked above there are README contents displayed that explain how the "union" files are made.

ADD REPLY
0
Entering edit mode

Thanks. I still don't understand why the output is so different: the output of hg19_Ensembl_gene.bed enter image description here the output of hg19.Ensembl.union.bed enter image description here

ADD REPLY
0
Entering edit mode

I also run geneBody_coverage.py using different .bed files, hg19_Ensembl_gene.bed and hg19_RefSeq.bed and the same batch of bam files, but the outputs are very different:

geneBody_coverage.py -r hg19_Ensembl_gene.bed \
-i split_bam/short_genes/bams/ \
-o geneBody_coverage/short_genes/Ensembl_gene

enter image description here

geneBody_coverage.py -r hg19_RefSeq.bed \
-i split_bam/short_genes/bams/ \
-o geneBody_coverage/short_genes/RefSeq

enter image description here

I am very confused about why they are so different.

Which bed file should I use?

Thanks

ADD REPLY

Login before adding your answer.

Traffic: 2424 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