Question: what does the bed file in bedtools multicov command refer to ?
0
gravatar for rania.hamdy1
12 months 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 12 months ago by ATpoint34k • written 12 months 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 12 months ago • written 12 months ago by geek_y11k

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

ADD REPLYlink written 12 months ago by rania.hamdy10

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

ADD REPLYlink written 12 months ago by geek_y11k
0
gravatar for ATpoint
12 months ago by
ATpoint34k
Germany
ATpoint34k wrote:

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

ADD COMMENTlink modified 12 months ago • written 12 months ago by ATpoint34k

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

ADD REPLYlink written 12 months 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 12 months ago by ATpoint34k
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: 1018 users visited in the last hour