How To Calculate Bwa Mappability Score?

Question

[lyz10302012]

Hello,

Can anyone tell me how to calculate bwa mappability score with specific parameters? Is there such existing tools?

Thanks

Eric

[Pierre Lindenbaum]

how do you define mapability ? I only know http://genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=141011952&c=chrX&g=wgEncodeMapability

[Devon Ryan]

Do you maybe mean a read's MAPQ score rather than mapability?

Answer 1

[Chris Miller]

Calculating mapability is really straightforward:

1. Take a reference genome's fasta file and generate all possible reads of a given length.
2. Map those reads back to the reference genome, and identify which of them map uniquely in the expected location.

You'll find some very simple bash/perl code that wraps BWA and then calculates mapability and gc-content scores in 100bp windows here: https://code.google.com/p/readdepth/downloads/detail?name=createCustomAnnotations.v1.tar.gz

Keep in mind that this code only generates single-end reads, but by using paired-end reads, you up your mapability significantly. The code could easily be tweaked to do that. Be sure to think about what insert size is appropriate.

There are also some precomputed mappability tracks available for download through the UCSC table browser.