This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Adjusting for uneven Nextseq loading

Ran a ChIP-seq experiment with cells treated with vehicle and a drug, so sequenced 4 samples (vehicle input, vehicleIP, drugInput, drugIP). During multiplexing the vehicle input was unevenly added to the pooled libraries, and following sequencing therefore had a much bigger fastq file. I am wondering about usability of these results and how to potentially adjust during analysis for this uneven loading. Thank you.

p.s. The 3 smaller fastq files still gave a good amount of aligned reads so thats not an issue. Just thinking in terms of peak calling, differential peak analysis etc.

nextseq fastq sequencing chip-seq

1 answer

Just downsampling the big fastq might be a good plan. You can use samtools for this. I'd visualize all the results too before analysis.

I find Chip-seq to be notoriously dodgy.

So you would just subsample at different % until it hits around the size of the other BAM files or is there another downstream metric?

Yep, sounds about right!

Log in to answer this question.