Removing Duplicate Reads Post Alignment
Hi,
I'm analyzing RNA-Seq data and I want to check for snp in it. For that, I want to remove duplicated reads from the bam file. Duplicate reads are reads that are aligned at the same position on the genome. Could anyone help me with that ?
Thanks,
N.
duplicates
reads
• 4,811 views
•
link
written
by
Nicolas Rosewick
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Search A Dna Sequence Motif In Bam File?
written by newDNASeqer<p>I need to find a short DNA sequence motif from a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file, because the wet biologist found a potential mutation from Sanger sequencing. …
-
Removing Duplicates Reads (>= 5)
written by vj<p>Hello,</p> <p>I am having trouble finding softwares that can remove duplicate reads from the aligned files (single-end <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a>/bam or bed files). I intend to …
-
The Reference Genome Of The Reading In Bam File
written by shl198<p>Dear all, I have a question about analyzing <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files with R. In <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file, it includes all the reads. Are the reads …
-
Picard Mark Duplicate Reads In Galaxy
written by Tonyzeng •<p>HI, I want to use picard mark duplicate reads to just mark my read with potential duplication. from the manual of picard, we need to …
-
Samtools Filter Reads Cigar Field
written by dfernan<p>Hi,</p> <p>I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file that I'd like to filter. I want to filter out all the reads that were aligned to intronic …
-
Rna-Seq : How To Check For Very High Coverage Region
written by Nicolas Rosewick<p>Hi,</p> <p>How can I check for very high coverage region in a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. I've several samples with good alignment rate (~90%) but when …
-
How Do I Pull Out One Read (Given Its Id) From A Bam?
written by zam.iqbal.genome<p>Hi there</p> <p>Does anyone know how I can pull just one read from a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file, when I know the read-id (and sequence) of …
-
How To Fix A Bam-Readcount Sm Error: "Couldn'T Grab Single-End Mapping Quality For Read"
written by amnonim •<p>Hi I'm using <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a>-readcount to get coverage per nucleotide as part of our exome-seq variant annotation pipeline, following <a href='http://bio-bwa.sourceforge.net/'>bwa</a> and gatk when i …
-
Chip-Seq Peak Calling Using Macs For Replicates In Galaxy
written by anon111 •<p>I'm new to Galaxy. I have Chip-seq data for two cell types. Each cell type has two replicates. I already ran an alignment workflow in …
-
Running A Somatic Sniper With Marked Duplicate And Sorted Bam File.
written by kulkarni.aishwarya87 •<p>Hi for running the <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file through somatic sniper, do we need a marked duplicate <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file? Thanks</p>
please read Removing duplicates from RNA-seq data
furthermore, please search "Removing duplicates from RNA-seq data", and please avoid duplicate posting
This has been answered several times on the site. Briefly, two options for software are samtools and picard. As far as I know, If you have paired-end data, picard is preferred. You should have a think about whether it's appropriate to remove duplicates in your situation. See link provided by Nandini.
Closing this question as a duplicate.