How do I find the same gene in two different assemblies?
0
0
Entering edit mode
5.6 years ago

I have two genome assemblies from two phenotypes of a non-model ant species. I have performed gene prediction on them both and have found a difference in the number of predicted genes. How do I find these differing genes? I can't compare gene IDs as they will be different in both assemblies and I can't pull their positions and compare them as the same gene has different coordinates in both assemblies simply due to the assemblies being different.

Any ideas on how I can compare these putative genes between the two genomes?

RNA-Seq Assembly alignment genome • 1.3k views
ADD COMMENT
4
Entering edit mode

Maybe blast can help?

ADD REPLY
1
Entering edit mode

For non coding genes, UTRs etc it would also be good to generate a chain file for liftOver (provided only assemblies are different and species are same/very_similar) . I have used this pipe line https://github.com/wurmlab/flo with good results to generate chain files which then can be used by crossmap/liftOver to convert coordinates.

ADD REPLY
0
Entering edit mode

I also think generating your own chain files and use liftOver is a good shoot. Otherwise, extract the genomic sequence of each gene and then use blat or gmap to map the to the other assembly.

ADD REPLY

Login before adding your answer.

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