Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
ATAC-seq primers Fasteno Sample Prep / Library Generation 14 02-06-2020 06:06 AM
ATAC-seq protocol wen yuan Sample Prep / Library Generation 69 02-22-2018 08:01 AM
ATAC Seq data analysis Sciurus Bioinformatics 22 05-04-2016 12:15 PM
How to use DANPOS2 with ATAC Seq data? Sciurus Bioinformatics 0 06-22-2015 05:34 AM
ATAC-seq QC using gel Fasteno Sample Prep / Library Generation 0 03-12-2015 04:13 AM

Thread Tools
Old 08-26-2015, 08:00 AM   #1
Junior Member
Location: CT

Join Date: Aug 2015
Posts: 3
Default ATAC-seq cutting frequency analysis

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!!
Attached Images
File Type: png Screen Shot 2015-08-26 at 10.55.57 AM.png (47.5 KB, 60 views)
cloverliang is offline   Reply With Quote


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 07:05 PM.

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