This is a test version of Biostars. For the public version, visit https://www.biostars.org.
multiBamSummary

Hi

I want to count reads in defined intervals using multiBamSummary from deepTools. I want to count reads only on the reverse strand. Will --samFlagInclude 16 work? And when I want to count reads on the forward strand will --samFlagExclude 16 work ?

deeptools multibamsummary reverse-strand

1 answer

Yes, --samFlagExclude 16 will exclude reads in the reverse strand, you need --samFlagInclude 16 (primary hits, reverse strand) or some combination for the flags

Log in to answer this question.