Calculate Average Genome Coverage Accurately for hg38 Reference Using Alternate Contigs and Decoys
I have whole genome sequencing data aligned to hg38 and alternate contigs and decoys (the one from bwakit). What command could I use to accurately calculate the average genome coverage of each BAM file?
• 1,330 views
•
link
written
by
dario.garvan
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Gene performance/coverage calculation methodology in .bam files using samtools
written by slzr_ •
-
Extract human, mouse and pig aligned reads from .bam file alignment to conactenated genomes
written by daewowo •
-
genome coverage visualization
written by gubrins
-
Comparative analysis of experimental and control data from BS seq run
written by kspata
-
Efficent Estimation of Sequencing Coverage
written by dario.garvan
-
How to calculate average depth for each Amplicon in targeted amplicon sequencing?
written by bioinfo89 •
-
What tools can I use to calculate sequencing depth and sequencing coverage for ddRADseq data
written by jingjin2203 •
-
Coverage and fragment length
written by aaron.agri •
-
Best practice to identify Copy number alteration from mouse whole genome low coverage sequencing st…
written by ifudontmind_plzz
-
Split BAM by average coverage
written by h.mon
Tools To Calculate Average Coverage For A Bam File? and the search function please. Has been asked many times before.
Yes, but I don't think any of those tools explicitly handle the case of alternate contigs. Somehow, the main chromosome and its alternate contigs need to be merged before calculation happens. Otherwise, the coverage of the main chromosome is artificially reduced if the reads map to alternate contigs.