Hi there,
Does anyone have any suggestions for a tool to calculate the mapping coverage (e.g. for each contig in the reference genome, how many reads are mapped against this contig)?
Also if any one can provide me with references covering this topic, I will be gratefully thankful.
thanks in advance.
Shazly
•
written
by
Hatem Elshazly •
2 answers
Using the BBMap package, you can do this with the Pileup tool:
pileup.sh in=mapped.bam out=stats.txt
That will give the number of reads mapped to each reference sequence, average coverage, percent of bases covered, standard deviation, etc.
Hello Shazly,
bcbio is a "python toolkit providing best-practice pipelines for fully automated high throughput sequencing analysis".
We want to add coverage to the pipeline and are considering a new project called chanjo. I haven't done due diligence yet but you might want to check it out.
Regards,
John
•
written
by
John Kern •
Log in to answer this question.
-
Different counts of mapped bases between Samtools coverage and Qualimap bamqc
written by pauline10albert •
-
Tools for Subsampling Paired-End Reads to 25X Coverage ?
written by Sony Nguyen •
-
VCF file transfer to a better genome
written by Nilo •
-
Bam file depth of coverage and statstics
written by daewowo •
-
Get per-nucleotide mapping position for each mapped reads from bam file
written by Changhao Li •
-
Mapping to duplicated region
written by eadele •
-
Poor Output from Spades with High Coverage Input Data
written by callamartyn •
-
Calculating coverage of genes which have accounted for multiple best hits from bam/sam files
written by cwwong13 •
-
Whole Genome Sequencing (WGS) Alignment Quality check softwares/tools
written by ravi.uhdnis
-
Split BAM into 1 BAM per contig
written by Marvin
See this link: Tools To Calculate Average Coverage For A Bam File?
Though this link is pretty old but you may still find some tools that may still exist. Search for coverage in the search window and you should be able to find lots of posts dealing with the same problem.
Hi,
I would share my two cents of experience. Its straightforward to calculate mapping coverage using Bedtools but what I missed initially is that its not feasible to apply mapping quality, base quality, number of covering reads and similar filters using a simple approach like Bedtools. Its much better to use GATK DepthofCoverage and CallableLoci which can apply all those filters and give a "true" representation of coverage & depth. Look into their documentation plz.