Genome mapping and annotation
1
0
Entering edit mode
16 days ago
claudia • 0

Hello,

I need to create a kind of "translation file" between two versions of the same genome. The two genomes use a completely different naming system. I would like to know what ID on genomeA corresponds to what ID on genomeB.

Is there an easy way or a tool for this?

I thought about two strategies but none really convinced me:

1- Blast genes from versionA to the versionB genome, keep the best hits, use a script to extract the genes matching each position from annotationB, considering a kind of sliding window since the stast-stop position of the blast hit, may not exactly correspond to the start position of a gene on the annotationB (pretty hard considering my coding skills).

2- Trying to align sequences of geneA on indexed genomeB with Bowtie (got errors, since I'm not working with reads).

Any ideas?

Thank you very much

annotation genome • 176 views
ADD COMMENT
0
Entering edit mode
16 days ago

Hi claudia

For 1 BLAST Reciprocal Best Hits (RBH)

For 2 liftoff

ADD COMMENT
0
Entering edit mode

Thank you I will try :)

ADD REPLY

Login before adding your answer.

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