Question: samtools vs bedtools
0
gravatar for xupbuy
2.2 years ago by
xupbuy10
xupbuy10 wrote:

Hi I have a silly question: What is the difference between samtools depth and bedtools genomecov? For both tools, the input file needs to be *.sorted.bam, which was generated by sam->bam->sorted.bam.

I tried both tools using the same *.sorted.bam input file and they generate the same coverage file.

Another question:

Can I use bedtools to map reads file (fastq or fasta ) to a reference genome? Like bowtie or bwa did? Because I like the function of bedtools that gives you how many reads map to the genome, and also the percentage of genome that is mapped by the reads.

Thank you.

sequencing alignment • 1.3k views
ADD COMMENTlink modified 2.2 years ago by genomax78k • written 2.2 years ago by xupbuy10

No, you cannot use bedtools as an aligner in any meaningful way.

ADD REPLYlink written 2.2 years ago by jared.andrews074.8k

Your question is not silly, but is poorly formulated, and suggests you made little effort to find the answers. For example, samtools depth vs bedtools genomecov has been discussed many times (e.g. avrilomics blog, on biostars or bioinformatics stackexchange). Please read Tutorial: How To Ask Good Questions On Technical And Scientific Forums.

ADD REPLYlink written 2.2 years ago by h.mon29k
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: 1983 users visited in the last hour