SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   ATAC-seq cutting frequency analysis (http://seqanswers.com/forums/showthread.php?t=62269)

cloverliang 08-26-2015 08:00 AM

ATAC-seq cutting frequency analysis
 
1 Attachment(s)
Hi everyone!

I'm analyzing ATAC-seq data. But I got problems in observing nice peaks as shown in Fig. 1c. in Buenrostro 2013 nature method paper.

I have the following steps:
1. Map to hg19 with bowtie2 (-X 2000).
2. Remove reads with MAPQ < 10 and sort bam file using samtools.
3. Remove duplicates using picard tools. (60%~70% duplication rate)
4. Merge sorted bam of four replicates together.
5. Get cutting position for each fragment and visualize cutting frequency (per basepair) in IGV.

I got the attached figure in the same position as shown in Fig. 1c. The first track is GM12878_ATACseq_50k. Data ranges from 0 to 6. There are some peaks but the noise-signal ratio seems really high. The second track is from ENCODE U Washington GM12865. Data range from 0 to 16. I wonder if I did something wrong or further smoothing/normalizing steps are needed.

Thank you!!


All times are GMT -8. The time now is 04:53 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.