I realise this has been asked a lot, but I havent found a tool that fits my simple needs.
I would like to calculate coverage % for an alignment file (sam ot bam) against a reference fasta genome. (NB there may be mutiple aligned sequences in the sam/bam file).
bedtools can calculate coverage - but splits this up and I cant work out how to get a single % value for coverage - other than writing a script to parse the output. Ditto for samtools mpileup - this generates a huge ammount of information, whereas I am after a single percentage.
One basic way to do this is plot in IGV, extract consensus, then take the number of N's divided by reference genome length. There maybe existing tools that can do this in one command.
1 answer
Unfortunately, nobody ever defines what they mean by 'coverage'. Can you please explain what you need such that another person's interpretation of what you mean by coverage can be made beyond a reasonable doubt?
Do you mean the percent of the reference genome that has at least 1 aligned read? - if 'yes', see my code, here: Determine % of reference genome covered by aligned SAM/BAM
Kevin
Log in to answer this question.
-
Bam file depth of coverage and statstics
written by daewowo •
-
Convert sam file with no header to fastq
written by daewowo •
-
Post Alignment basic QC of RNA-Seq data
written by daewowo •
-
RSEM generates zero byte size genomic bam file
written by arozhok •
-
To calculate the length of aligned genome
written by sambioinfo2018 •
-
metagenomic assembly coverage, for multiple samples
written by willnotburn •
-
IGV does not show BAM file information
written by Hans
-
Coverage of reference from the BAM file
written by brs1111 •
-
Tuning Consensus Calling with bcftools/samtools
written by nickp60 •
-
genome coverage by samtools
written by zizigolu
