Hi, Biostars.
I have a SAM file with some PCR duplicates mapped by bowtie2.
The original fastq reads of these PCR duplicates are not same sequence length,
so I cannot use "samtools rmdup" command.
I would like to extract unique position sam lines with all fields of orginal sam. Anyone has a solution? thanks.
----Example---
name1 0 chr1 124344 . . . . . AGTAGGTGGGG FFFFFFFFFFF AS:i:0 XN:i:0 XM:i:0
name2 0 chr1 124344 . . . . . AGTAGGTGGGGGATT FFFFFFFFFFFFFFF AS:i:0 XN:i:0 XM:i:0
↓
name1 0 chr1 124344 . . . . . AGTAGGTGGGG FFFFFFFFFFF AS:i:0 XN:i:0 XM:i:0
1 answer
The original fastq reads of these PCR duplicates are not same sequence length, so I cannot use "samtools rmdup" command.
Are you sure you cannot use samtools rmdup (or even better picard MarkDuplicates)? With single end reads these tools look only for the start position of a read to call it duplicate. (In fact I wrote my own program to look at both start and end position, MarkDupByStartEnd).
Log in to answer this question.
-
Best mapped reads from Bowtie2 alignment: AS>XS
written by Tash •
-
Duplicate reads remain in bam after picard MarkDuplicates --REMOVE_DUPLICATES true
written by kalavattam •
-
Command line .sam file processing
written by lydiaiglesias •
-
bowtie2 '-a' does not return all alignments
written by shim •
-
Extract unique and multimapped reads from bowtie 1 sam/bam output
written by Konstantinos Yeles
-
read strand identification from sam file
written by alireza346 •
-
bowtie2 messes up headers with sequences in output
written by chefarov
-
Extract uniquely mapped reads from bam file
written by User000
-
How to find Significant Trans- Interactions in HiC data
written by komal.rathi
-
Finding Reads Joining Exons From Tophat Output
written by hbw
