Questions tagged [markduplicate]
A tool from the Picard toolkit to mark or remove optical or PCR duplicates in sequencing data.
6
questions
2
votes
1
answer
228
views
MarkDuplicatesSpark failing with cryptic error message. MarkDuplicates succeeds
I have been trying to follow the GATK Best Practice Workflow for 'Data pre-processing for variant discovery' (https://gatk.broadinstitute.org/hc/en-us/articles/360035535912).
This has all been run on ...
- 21
0
votes
1
answer
94
views
Impact of merging ChIP-seq runs of the same sample on PCR duplicates identification?
I'd like to a follow up question to this question related to merging fastq files for ChIP-seq.
Let's assume we have one sample library that is re-sequenced two or three times in order to achieve a ...
- 1
0
votes
1
answer
119
views
Diploid Consensus Sequences
My goal is to generate a diploid consensus genome. By "diploid consensus", I mean that I want to merge all supporting reads for a genomic region into a single haplome such that the resulting ...
- 179
4
votes
1
answer
257
views
How does picard's MarkDuplicate handle unmapped reads?
Our BAM files are created according to a "lossless" alignment procedure [1] from the Broad Institute GATK documenation and involves re-adding the ...
- 319
4
votes
2
answers
372
views
Targeted NGS, up to 99% of reads have been marked as duplicates
Currently I'm performing whole analysis (pipeline from *.fastq to *.vcf) of 41 samples (targeted NGS). I rely on GATK best practices, however with some modifications. I decided to use the following ...
- 206
5
votes
1
answer
176
views
Does picard markduplicate toggle PCR duplicate samflag
I have a RNA-seq bam file and there are few reads that are puzzling me.
According to the bam header, this bam file is sorted by coordinate, created using tophat and markduplicate step is not done. ...
- 151