Is there a way to extract sites that are perfectly aligned to the reference genome in maf format?
0
0
Entering edit mode
8.3 years ago
Cacau ▴ 520

Perfectly aligned to means that the aligned site in the target genome is the same nucleotide as the reference genome.

alignment • 1.7k views
ADD COMMENT
0
Entering edit mode

What is the target genome and reference genome ? what is MAF format ?

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Hi , I had a similar problem and have a very imperfect solution. PHAST (http://compgen.cshl.edu/phast/) provides a utility called maf_parse which can get an maf column for a bed file entry. I created a bed file for entire genome with each entry being 1 base pair. Then I used maf_parse to get maf columns for this bed. The entries where the nucleotides match should be the ones you want.

ADD REPLY
0
Entering edit mode

That is awesome. I happen to think of using maf_parse too. Thank you a lot.

ADD REPLY

Login before adding your answer.

Traffic: 2724 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