Bamcoverage unique mapping reads

0

Entering edit mode

3.8 years ago

liang795 ▴ 10

Hi Guys

Would bamcoverage only calculate read-pairs if any pair maps to unique locations? As I my bam files were generated with multiple mappings, but for visualization I only want unique mapping locations. Does bamcoverage have this option?

Thanks!

RNA-Seq bamcoverage multiplemapping • 1.4k views

ADD COMMENT • link updated 2.9 years ago by tim.ivanov.92 ▴ 40 • written 3.8 years ago by liang795 ▴ 10

0

Entering edit mode

You can try the bamCoverage options --minMappingQuality, --samFlagInclude or --samFlagExclude to take into account only unique properly paired reads. How to set these parameters will probably depend on the program and arguments you used for mapping.

ADD REPLY • link 3.8 years ago by h.mon 35k

0

Entering edit mode

How can i specify the SAM file flag NH=1 into option --samFlagInclude? what is the SAM flag code for this?

ADD REPLY • link 2.9 years ago by tim.ivanov.92 ▴ 40

Login before adding your answer.