Question: ATAC data analysis
0
gravatar for krushnach80
23 months ago by
krushnach80500
krushnach80500 wrote:

I am dooing atac seq data analysis so far what i have done

  • I aligned the data with bowtie2
  • Removed duplicates
  • made bed files out of those bed file
  • Did peak calling with Homer using it findPeak command -Now i have list of peak files with respect to various condition . In Homer manual i read a function called findMotifsgenome function , i have some confusion here how do i proceed from here do i have to just give the peak file as it is or I have to process the peak files .

Any suggestion or help would be highly appreciated

atac-seq next-gen • 1.6k views
ADD COMMENTlink modified 23 months ago by Kasthuri240 • written 23 months ago by krushnach80500
1
gravatar for Kasthuri
23 months ago by
Kasthuri240
United States
Kasthuri240 wrote:

Typically we use diffbind to find differentially binding sites between the control (closed chromatin) and experiment (open chromatin). I don't use Homer. But one can use bedtools getfasta to get the fasta sequences for those open chromatin peak sites and use MEME to get the consensus motif if that is what one wants. There are other analysis one can do, like using deepTools etc.

ADD COMMENTlink written 23 months ago by Kasthuri240

okay so i have to convert the peak files into fast and then use MEME? thats what you are suggesting

ADD REPLYlink written 23 months ago by krushnach80500

Yes, that is one reasonable approach...

ADD REPLYlink written 23 months ago by Kasthuri240

i have all the peak files but what do i need from the peak file as input to to do MEME becasue it contains field like peak ID chromosome location , fold change etc etc .I dont think i need all those .Do i have to extract some columns from my peak files or i just have to use the whole file?

ADD REPLYlink written 23 months ago by krushnach80500

You need to have a bed file - just use the chromosome, start and end position and get the fasta sequences through bedtools and use MEME.

ADD REPLYlink written 23 months ago by Kasthuri240
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1283 users visited in the last hour