Question: what does the bed file in bedtools multicov command refer to ?
0
gravatar for rania.hamdy1
26 days ago by
rania.hamdy10 wrote:

Dear bedtools users,

Hope everything is well with you

I want to use multicov mode to get reads count on my chip seq reads

according to beedtools documentation site I've sorted and then indexed the bam file

bedtools multicov [OPTIONS] -bams BAM1 BAM2 BAM3 ... BAMn -bed  <BED/GFF/VCF>

but i don't understand what does the bed file refer to ??

anybody who can help me, please?

thanks in advance

ADD COMMENTlink modified 26 days ago by ATpoint17k • written 26 days ago by rania.hamdy10

I want to use multicov mode to get reads count on my chip seq reads

Where do you want to measure the read counts ?

ADD REPLYlink modified 26 days ago • written 26 days ago by geek_y9.7k

i want to reads count histone modification in chip seq file REMC Dataset

ADD REPLYlink written 26 days ago by rania.hamdy10

The chip-seq data is usually represented as peaks of genomic coordinates in bed file format.

ADD REPLYlink written 26 days ago by geek_y9.7k
0
gravatar for ATpoint
26 days ago by
ATpoint17k
Germany
ATpoint17k wrote:

It is the intervals to quantify over (= your peaks) in BED format.

ADD COMMENTlink modified 26 days ago • written 26 days ago by ATpoint17k

I am a newbie in bioinformatic researches i can't understand your comment Can you explain more?

ADD REPLYlink written 26 days ago by rania.hamdy10

You really should get a solid background first. Check for example the vignette (the manual) of csaw which is a Bioconductor package to analyze ChIP-seq data. It offers a good introduction to ChIP-seq. Also spend time to get familiar with the common file formats (BED, BAM/SAM, GTF, FASTQ etc.). In this case the BED file is the template to count the reads so it is the parts of the genome that bedtools uses to count the reads. It does not automatically screen the entire genome but only the part you specify, e.g.:

chr1      1       450

would be an interval on chromosome 1 that could be a peak. Peaks come from peak calling. See e.g. the macs2 software package for peak calling.

ADD REPLYlink written 25 days ago by ATpoint17k
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: 1093 users visited in the last hour