samtools vs bedtools
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.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Count strand-specific 5' / 3' coverage across whole genome for paired-end fragments
written by acvill •
-
How to remove reads mapped to both plasmid and chromosome from the .sam file?
written by Vasiliy Krestov •
-
Copy number (repeat count) of de nove assembly contigs
written by Caronkey •
-
Set the max coverage wanted with Hisat2 ?
written by Chvatil
-
How can I get the read depth of my polymorphism?
written by apl00028
-
How to get regions with mapped reads
written by nanoide
-
bowtie --fr/--rf/--ff reporting
written by ThePresident
-
mapping the filtered reads (unmapped reeads) on the genome using bowtie2
written by zizigolu
-
How do I obtain uniquely mappable reads from bowtie on galaxy?
written by Jason2 •
-
Bwa And Splicing Junctions
written by GPR
No, you cannot use bedtools as an aligner in any meaningful way.
Your question is not silly, but is poorly formulated, and suggests you made little effort to find the answers. For example,
samtools depthvsbedtools genomecovhas 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.