This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tab delineated file containing the number of reads mapping to each gene on chromosome 22. How I can do this with bedtools and shell commands?

After an RNA-seq experiment, I got the alignment results in 6 bed format files located in a folder. I'm trying to make a tab delineated file containing the number of reads mapping to each gene on chromosome 22 in each experiment. I will appreciate very much your help. How I can do this with bedtools and shell commands?

Thanks

rna-seq alignment sequencing rna-seq
• 1,739 views • link
updated by Benn Benn • written by Kevin_Smith • Kevin_Smith

1 answer

What you need is the bam files, and a bed file with the coordinates of interest. Then in the unix shell:

bedtools coverage -abam file.bam -b file.bed > results.cov.txt
• 0 views • link
written by Benn Benn

Log in to answer this question.

More posts like this