This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get the count.matrix from Chip-Seq BAM files

Hi guys, is there a way to get the count matrix from BWA-BAM files of a Chip experiment? Normally this is the standard pipeline in RNA-Seq but this is the first time my boss asked me to do this for Chip-Seq data. normally I run MACS and then I go on. I tried using featureCounts but as far as I know I have to specify the GTF.featureType (eg. exon). How if I don't want to specify the type of feature? Thank you in advance

chip-seq
• 2,089 views • link
updated by ATpoint ATpoint • written by elb elb

1 answer

featureCounts accepts the SAF format:

awk 'OFS="\t" {print $1"_"$2"_"$3, $1, $2, $3, "."}' in.narrowPeak > out.saf

The first column is a name column, I simply use a concatenation of chr_start_end. This can go into featureCounts with -F SAF.

• 0 views • link
written by ATpoint ATpoint

Okk I'm getting the point. Is it equal with the *broad_peaks.bed file (as from MACS) instead of narrowPeak to generate the SAF file? I suspect yes but I'm not sure. The coordinates (i.e. start-end) should be the same.

• 0 views • link
written by elb elb

SAF is Name-Chr-Start-End-Strand

• 0 views • link
written by ATpoint ATpoint

Yes, sure! thank you very much!

• 0 views • link
written by elb elb

Log in to answer this question.

More posts like this