Hi,
I'm using python and I have a wheat genome portion in .FASTA format (nucleotides).
I want to run BLAST against ENSEMBL and obtain genes annotations in .GFF3 programatically in order to use those items in a custom pipeline.
Any ideas?
Thank you in advance
EDITED:
My vision is:
•
written
by
juanma_lace •
Identifying genes in eukaryotes is a bit more complicated than simply aligning proteins with BLAST. If you have a particular set of reference proteins or ESTs, you can splice align these (with programs like GeneSeqer and GenomeThreader) to get a good first approximation of gene structure.
You can also use gene predictors like SNAP or AUGUSTUS to predict genes. These do not require and protein or transcript sequences, but they typically have a substantial number of false positives.
Tools like Maker and EVM combine these two approaches (spliced alignment and ab initio gene prediction), and produce much more reliable annotations. However, they are a pain to set up and more complicated to run.
Sorry there is no easy answer to this question!
•
written
by
Daniel Standage
I see, thank you for your answer. Upvoted
•
written
by
juanma_lace •
Log in to answer this question.
¡Bienvenidos a BioStar! We'd be happy to help, but it looks like you need to take a bit more time to clarify what you're doing. GFF3 is a format that can encode a wide variety of genomic features. From your question it is not clear what features you are looking for or how BLAST/ENSEMBL will help you identify those features.
Thanks, I've edited the question