Question: determining length of fasta cds and protein sequences
0
gravatar for vigneshprbh37
3.8 years ago by
INDIA
vigneshprbh3720 wrote:

I have performed blastn and tblastn executions for a list of selected genes. for purpose of analysis i need to compare the best hit alignment length with best percent identity to the actual length of the sequences of cds protein.

is there a bioinformatic tool i can use for accomplishing this

 

is there any algorithms i can use on linux to accomplish the same

fasta protein cds software • 1.2k views
ADD COMMENTlink modified 14 days ago by Biostar ♦♦ 20 • written 3.8 years ago by vigneshprbh3720
1
gravatar for lieven.sterck
4 months ago by
lieven.sterck3.3k
VIB, Ghent, Belgium
lieven.sterck3.3k wrote:

Assuming you are using the most recent version of blast (and you should) , you can ask to add the query (or hit) length to be added in certain output formats, such as the tabular one.

you'll need to add the following to your blast cmdline:

-outfmt "6 std qlen slen"

this will add, to the normal output the query input seq length and the hit seq length in the second-to-last and last column respectively .

See the blast help (blastp -help) for more info on those parameters (and how to add others for instance)

ADD COMMENTlink written 4 months ago by lieven.sterck3.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2205 users visited in the last hour