Merge replicate fastq files of the same sample together?
1
0
Entering edit mode
15 days ago
foxiw • 0

I have 13 samples which were initially sequenced and expected to yield around 30 million reads per sample, however we only managed to get half of that. Therefore we re-pooled the samples and ran them again, getting more reads this time. Therefore I have 2 files per sample, run 1 (with low reads) and run 2 (with high reads). I tried merging these files together (using cat function) to get even greater depth but my QC analysis shows a lot of duplications. Is there a better way of merging these files whilst avoiding high levels of duplications? Thanks!

fastq RNA-seq • 157 views
ADD COMMENT
2
Entering edit mode
15 days ago
GenoMax 103k

RNAseq data is expected to have some duplication. So if you are going based on FastQC result of something showing up as "overrepresented" then that is likely not a big cause of worry. You should go ahead and analyze the data.

Is there a better way of merging these files whilst avoiding high levels of duplications?

If you do have a problem with duplicates (coming from one too many PCR cycles) then nothing is going to fix that at this stage.

ADD COMMENT
0
Entering edit mode

Thanks.

After using multiqc on my html files, my Sequence Counts shows me that there is a lot of duplications per sample. However, it is making me think that multiqc is considering the sequences from run 2 as 'duplicates' of run 1, which means I haven't merged the files properly. Is is normal then?

ADD REPLY
0
Entering edit mode

cat'ing the files simply tacks contents of second file at end of first. So as far as FastQC/MultiQC is concerned that is just one set of data. That in itself is going to have no effect on duplicates per se. Your sample is going to contain that duplication (if present). At this point you can't fix that part (if sample was overamplified for example). Hopefully all samples in this dataset underwent an identical treatment so there will be no experimental bias.

ADD REPLY
0
Entering edit mode

Thank you!

ADD REPLY

Login before adding your answer.

Traffic: 2581 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6