Hello Brian! Thank you for suggesting BBTools. I will try this!
• 0 views
•
link
written
by
Shakunthala Natarajan •
Understanding bedtools coverage for datasets with in-dels.
Hello. I have long reads (ONT) and I have mapped them to a reference genome. I have the sorted bam file in hand. As a next step, I wanted to calculate the read coverage depth per position using bedtools with the command:
bedtools coverage -d -split -ibam sorted_bam_file > output_file
I want to further understand how bedtools deals with insertions and deletions, since my reads are long reads (ONT).
Can someone please help me with how to proceed further to understand how in-dels are dealt with during coverage calculation with bedtools as mentioned above?
Thank you!
• 2,094 views
•
link
written
by
Shakunthala Natarajan •
1 answer
You can calculate positional coverage, including indels, using BBTools, like this:
pileup.sh in=mapped.bam basecov=basecov.txt
It has various options that you can see by running the shellscript with no arguments, but the most salient is "delcoverage=t" which is the default. That means bases covered by a deletion event will have their coverage incremented. You can set it to false, and then deletion events will not have their coverage incremented. There's no similar option for insertions, because those bases don't exist in the reference, so obviously it makes no sense to classify them as covered based on an insertion.
This does not require a sorted bam file, it works fine with an unsorted sam file too.
• 0 views
•
link
written
by
Brian Bushnell
Log in to answer this question.
More posts like this
-
Differences between SNPs and Allele coverage for Prediction
written by Loïc •
-
calculation of gene coverage
written by ameya4microbiology •
-
Count strand-specific 5' / 3' coverage across whole genome for paired-end fragments
written by acvill •
-
Merge regions in bedtools genomecov/bedgraph file
written by eoneill627 •
-
genome coverage visualization
written by gubrins
-
Find genetic duplicates/deletions via coverage
written by AlBaars •
-
How to reassembly using .bed file
written by p_puppo •
-
Trying to check plasmid alignment
written by gracie •
-
Per base statistics on an RNAseq bam file
written by tassa.saldi •
-
coverage on exons with bedtools genomecov
written by guillaume.rbt
You may want to try
pandepthhttps://github.com/HuiyangYu/PanDepth for this. You probably don't need per-base coverage since the files would become gigantic.mosdepthcan generate per-base coverages and has similar similar functionality: https://github.com/brentp/mosdepthBoth these programs would likely be much faster than
bedtools.Hi thank you for bringing this to my notice! I would try this. But presently, I am using bedtools as a part of a script that I am running. To understand the output of the script's working, I need to understand how coverage is being calculated by bedtools when in-dels are present. Essentially, I want to back track from my results, get to the details of how coverage calculation is done by bedtools in the presence of in-dels. Is there some way in which I could do this?
or you can use
samtools depth in.bamHi thank you for suggesting samtools. But I would like to understand the coverage calculation by bedtools in the presence of in-dels as mentioned above.