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 desired read coverage. The resulting output of the sequencing is 2 or 3 fastq files for one individual sample.
If one has to mark duplicates (for example using Picard's MarkDuplicates) should the sub-samples be merged at the fastq level or at the bam file level (post alignment) after flagging duplicates before the merge?
If one takes the 2-3 sub-fastq files align them and then mark the duplicates in the bam files if the sequencing is shallow one could end up with not calling reads that are actually PCR duplicates of the library preparation. While in the other scenario, a fastq merge might result in reads to be marked as "duplicates" just because they appear 2-3 times, but maybe that read is a single reads that I sequenced 2-3 times, individually in each run.
My point is that I would expect this 2 different approached to lead to different outcomes. Probably this is a negligible effect in a overall analysis, but maybe it might actually have an impact in certain biological conditions (I'm thinking of heterochromatin silencing of repeat elements in the genome). I'm going to try it out as soon as I have a change, but I was wondering if this was discussed already somewhere else (I couldn't find much relevant) or if someone has some ideas on how to tackle this.
Let's assume that the reads are 50bp single end.