This is a test version of Biostars. For the public version, visit https://www.biostars.org.
obstain 100-way alignment of 100 species for a certain gene

I want to get the alignment file (Amino acid for protein, and nucleotide for CDS) for a particular gene from 100-way alignment of 100 species (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/multiz100way). Is there an easy way to get it?

Thanks

multiple-sequence-alignment 100-way-alignment
• 2,503 views • link
updated by Ram Ram • written by qwzhang0601 qwzhang0601

Thank you for your suggestion. I tried to get the alignment, but it seem it generate the alignment on the gene DNA sequence level (including noncoding part), right? I see many long gaps in the alignment. Is there a way to only include CDS sequence or protein sequence? Thanks

• 0 views • link
written by qwzhang0601 qwzhang0601

1 answer

You can find those alignments here. You will need to download the entire set and then find the gene/exons you need in that file.

It may also be worth looking at NCBI's homologene project where you would be able to search for one specific gene.

• 0 views • link
written by GenoMax GenoMax

Thanks. I got the files (like refGene.exonAA.fa.gz) from the link you provided. But it seems, it only includes alignment on exon level. We need to merge all exons by myself. Am I right? Not sure whether there is a direct way to do this, or any tool can help to merge alignment of exons into alignment of a gene.

I also looked at NCBI homologenes, but only found alignment for 10 species.

Thank you again for your suggestions.

• 0 views • link
written by qwzhang0601 qwzhang0601

Log in to answer this question.

More posts like this