How to remove Mate mapped to a different chromosome
1
0
Entering edit mode
4.4 years ago

Hi all!

When I generated the alignment statistics using Samtools Flagstat, I came across "mate mapped to a different chr". I was advised to remove these reads by my senior because then such reads will might introduce ambiguity at the time of peak calling. I checked out the various posts that underline the usefulness of this step on Biostars. However, none of them mention how to remove them.

Can someone help?

ChIP-Seq alignment samtools next-gen • 1.4k views
ADD COMMENT
1
Entering edit mode
4.4 years ago
ATpoint 80k

Extracting Paired-End Reads Sitting In Different Chromosomes

Peak calling by macs won't /should not be affected because macs only counts properly-paired reads which these aren't.

ADD COMMENT
0
Entering edit mode

This post does not address their removal from a bam file; instead, it addresses their extraction to a new file.

ADD REPLY

Login before adding your answer.

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