Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Given BAM/SAM file, how to see if it's single-end or paired-end sequencing? xxatbio Bioinformatics 2 08-11-2013 02:51 AM
How to count number of mapped paired-end and single-end rna-seq reads repinementer Bioinformatics 8 01-06-2013 05:06 AM
50 bp Single end Vs 100 bp Single end Vs 50 bp Paired end dhanapala RNA Sequencing 4 06-08-2012 05:09 PM
RNA-seq: Replicates, single-end, paired-end story pasta Bioinformatics 2 07-04-2011 11:51 PM
Single or Paired-End ChIP ScottC Illumina/Solexa 0 07-06-2009 02:38 PM

Thread Tools
Old 11-22-2013, 05:24 AM   #1
Location: UK

Join Date: Jun 2011
Posts: 61
Default Combination of paired-end and single-end samples in Chip-seq TF study

I have 2 batches of chip-seq samples:

(A) One biological SE replicate

This batch, actually, consists of 4 SE Chip-seq samples - one treatment, one control, both have IP controls. All are SE sequenced at a lab A

(B) Two biological PE replicates

This batch consists of 8 PE Chip-Seq samples - 2 treatment, 2 controls, all have IP controls. All are PE sequenced at a lab B.

The idea is to establish differential binding between the treatment and the controls. I plan to use macs and diffbind but I am open to suggestions.

My first specific question is how to treat the PE samples. I see 3 options:

(1) Treat all forward and reverse reads as SE reads

(2) Take just the forward or the reverse reads

(3) Treat each PE sample as, essentially, 2 technical replicates (one consisting of forward reads, one of reverse)

I find (3) intuitively attractive but, again, I am open to suggestions. If I go for it than the model would need to introduce 2 additional factors - one accounting for the technical replicates, another for the batch effect. So far I have dealt with batch effects only - I presume that adding additional factor shall be straight-forward but, please, let me know if there are things that I need to pay attention to.

Last edited by feralBiologist; 11-23-2013 at 10:51 AM.
feralBiologist is offline   Reply With Quote

chip-seq, diffbind, macs, paired end, single end

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 12:10 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO