This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get Fasta Sequences with their gene name for many genes?

Hi all,

I have a list of many gene name (about 5000) and I want to retrieve their corresponding fasta sequence. To this end, I try to use biomart, but it give me just nucleotide sequence. In fact, we cannot select more than one option in the "Attributes" part and we choose "sequence" in this part, so I cannot have their associated gene name for these sequences. Could you explain me how to have both nucleotide sequence and their related gene name. any suggestion would be highly appreciated.

gene sequence
• 6,237 views • link
updated by Matt Shirley Matt Shirley • written by seta seta

2 answers

I think this query is what you are looking for. You just need to select the gene name in the header information for the returned data file.

• 0 views • link
written by Matt Shirley Matt Shirley

you can have your gene in bed file (if not then download genes.gtf and convert it to genes.bed

first coloumn 'Chr' , second coloumn 'Start' third coloumn 'End' and then 'strand'

if you have whole genome fasta e.g. genome.fa

just use bedtools to extract fasta

getFastaFromBed -fi genome.fa -bed genes.bed -fo genes.fasta.out

details are here

http://bedtools.readthedocs.org/en/latest/content/tools/getfasta.html

• 0 views • link
written by Manvendra Singh Manvendra Singh

Log in to answer this question.

More posts like this