I want to add paired end support to my chip-seq software*, but to do that I need some paired end data. Do you know of any public PE chip-seq data? If you have some you could share a tiny subset of, that would also be greatly appreciated.
If you know a way to create fake PE bam files I'll accept that too.
•
written
by
endrebak
3 answers
Did you read the paper ? or saw where the dataset was uploaded ? The dataset are uploaded on GEO with GSE27679 id. Read the description for each sample. It says : Average pair-end library size: ~300bp . Go thru the paper and it is written quite explicitly about the paired end samples.
There's some good quality PE mouse ChIP seq in this H3K9me3 study
Hi,
We used paired end for a bunch of histone modifications here:
http://biorxiv.org/content/early/2016/05/19/054387
I have not uploaded it to SRA yet.
You can email me if you want to use it mbusby at broadinstitute dot org and motivate me to do it if it is what you need. I think most of them are 2x36 but they are all paired end.
Paired end is better!
Michele
Log in to answer this question.
-
Differential binding ChIP-seq/RNA-seq analysis
written by dcheng1 •
-
Dataset request: RNA-seq experiment that uses both single and paired end data with multiple experim…
written by endrebak
-
Strange peaks in paired-end ChIP-Seq independent on Abs but not single-end or input
written by greulich.fr •
-
Current status of EPIC peak caller
written by ariel.balter
-
Create bigwigs in python?
written by endrebak
-
ChIP-Seq: When reducing paired end reads to one coordinate, should I use the midpoint?
written by endrebak
-
Different mappability/effective genome fraction with paired end reads? (ChIP-Seq)
written by endrebak
-
Tool: Script to compute the effective genome size: epic-effective
written by endrebak
-
Paired-end data with different length
written by marimaria963 •
-
Sorting Fastq Files After Trimming (Orphans And Pe)
written by neal.platt
Here you go: http://www.ebi.ac.uk/ena/data/view/SRX1485421
http://www.ebi.ac.uk/ena/data/view/SRP062397
Human: http://www.ebi.ac.uk/ena/data/view/SRR1573750
I have tried a few datasets from sra, but they were not PE, despite the labels. I'll look into those, thanks.
The last one looks interesting, but I could not find the input data (using mobile, that might be it).
No guarantee they are something you could use. They are labelled ChIP-seq and are PE :-)
Yeah, those archive data are often terribly badly labeled. Takes ages to download and then you see that they cannot be used/aren't what you expected after all. Will perhaps try making my own, unless anyone has any datasets with both chip and input (only one of those has to be pe).
Score! The last file was input. Thanks! Ps. other answers still appreciated :)
Damn, with hisat2 you can just give the SRA-nb on the command line and it downloads and aligns it for you. Sweet:
hisat2 -x hg38 --sra-acc SRR1573750 -S chip_SRR1573750.sam