Thank you for your help. I used Bowtie2.
• 0 views
•
link
written
by
lisadavic66 •
Remove some reads from a BAM file
Hi,
I have a paired-end BAM file and I found that some pairs/mates have identical positions on the chromosome, i.e., the numbers in the POS and PNEXT fields are exactly the same. I want to remove these reads.
I also found that some reads present more than 2 times in my BAM file (same QNAME) (I think paired-end reads should have two records/rows in a BAM file). I want to keep the first two records/rows and remove the others.
I used samtools rmdup to remove duplicates.
Any help will be greatly appreciated. Thank you.
• 4,382 views
•
link
updated
by
Ram
•
written
by
lisadavic66 •
•
written
by
lisadavic66 •
1 answer
For your first task, picard's Mark Duplicates might do better job than rmdup. For your second task, maybe you have some supplementary alignments? Those would show up as multiple rows. Try removing them with samtools:
samtools view -bh -F 2048 input.bam >output.bam
Also please add which mapper you used, that will be helpful.
• 0 views
•
link
written
by
Biomonika (Noolean)
Is there any simple way, e.g., a Linux command solution, for these problem? Thanks.
• 0 views
•
link
written
by
lisadavic66 •
Log in to answer this question.
More posts like this
-
Remove all entries with duplicate names from fastq file?
written by wormball •
-
Count the fragments from BAM file
written by Lolla •
-
Removal of PCR dupliactes from trimmed reads
written by deepti1rao •
-
Picard vs samtools rmdup
written by xd_d
-
How do tools parse the QNAME?
written by John
-
Identifying mate pair reads from SAM file without using QNAME
written by StarCute
-
Calculate the number of paired-end reads their middle size is greater than # KB
written by lisadavic66 •
-
Find mates between all alignments in a SAM file
written by Pfs
-
Remove Duplicates Workflow
written by Anjan
-
Removing Pcr Duplicates From Bam Without Collapsing Sequences By Position?
written by user
