hi all,
I need to pull out a specific protein from a transcriptome data deposited in a local server (using linux platform). I have the protein sequence in fasta. I get to know from few papers that i need to use local blast and scan against all protein set.
Can anyone tell me how to use unix bash scripting to pull out the protein from the transcriptomic data by providing a sample scripts how to perform the task? If my Idea is wrong please advice me on a preferable steps. Your favor is much appreciated. thank you
1 answer
What you probably wanna do is first translate your transcriptome data into proteins data and put it in one big fasta file. This fasta file will be used to make a local db on your computer.
Then, use blast with your specific protein fasta file against this db.
There is a nice manual about how to use blast on unix/linux.
Log in to answer this question.
-
Tool: Neurosnap: Platform for making access to bioinformatic tools as simple as possible
written by keaun •
-
How to copy a file from a two layer ssh linux
written by Alex S •
-
Transfer annotation from fasta headers to associated gff
written by EJB •
-
RNA-Seq Differential Expression Bash Scripting Error
written by timothyfisher101 •
-
Flat file for functional annotation of draft fungal genome
written by minions-b •
-
aligning with STAR
written by snp87 •
-
Stand alone Blast annotation with multiple databases
written by Biogeek
-
SNP microarray dataset repositories
written by piscean1987 •
-
Sequencing analysis server
written by munnadt2 •
-
Variant Discovery Paper To Repeat?
written by Alice
