Which Mappability Track To Use?
1
3
Entering edit mode
12.5 years ago
Biomed 5.0k

There are a number of mappability tracks in the UCSC browser many of which have different/similar k-mer versions. I choose the ones that are similar to the k-mer lengths we use in our mapping process but my question is: What are the main issues when deciding which one of these tracks to use? Are they simply interchangeable if they have the same k-mer option? What is your heuristic to choose from among those tracks when you use them?

I would like to see if there is a relationship between the mappability of a region (as defined by the track) and the genotypes I detect using exome sequencing (illumina)

Thanks

ucsc • 5.3k views
ADD COMMENT
1
Entering edit mode

Almost certainly there is a relationship. Actually you do not need a track. The mapper already tells you which region tend to be repetitive.

ADD REPLY
0
Entering edit mode

Can you give more information on what you're using these mapability tracks for?

ADD REPLY
0
Entering edit mode

I would like to see if there is a relationship between the mappability of a region (as defined by the track) and the genotypes I detect using exome sequencing (illumina)

ADD REPLY
0
Entering edit mode

Almost certainly there is a relationship.

ADD REPLY
4
Entering edit mode
12.5 years ago

You might want to look at this paper: http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0030377

I know that there is a group (different from the one in the paper above) who are developing a general tool for generating mappability statistics for your particular genome and reads. I will update here when I hear that it has been released.

ADD COMMENT
0
Entering edit mode

I briefly looked at the paper. Quite nice. I would probably go for that if I had constructed my own mask.

ADD REPLY
0
Entering edit mode

Is this not the GEM mappability track that is available inthe UCSC genome browser? If so I guess the answer is to use the GEM track as opposed to broad, duke or Umass tracks.

ADD REPLY

Login before adding your answer.

Traffic: 2480 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6