Question: How to run PSSPred on multple sequences
gravatar for Tom
2.8 years ago by
Tom30 wrote:

Hi all, I'm using PSSpred link here to predict secondary structures for a set of peptides.

I am running the command on linux. When I run the command on a fasta file with a single sequence in it like this:

 perl test.fasta

it works (i.e, it correctly predicts the secondary structure for this one peptide and produces a file with the output).

However, I want the input fasta file to have multiple peptide sequences, and ideally I would like the output to be one predicted secondary structure per peptide.

When I add in three sequences to the test.fasta file, there is still one output file, and all of the sequences have been concatenated to one long string, instead of per peptide.

For example, if the input file is like this:


The output file shows the secondary structure like: AAABBBCCC.

I am trying to somehow get a predicted secondary structure per peptide; either one separate output file per peptide, or all to one file, but with a clear delimiter that I can parse, to show where the secondary structure for one peptide ends and the next begins.

I tried:

 for i in 'cat input_file'; do perl $i >> $i.out ; done

but that obviously doesn't work. Does anyone know what the specific command is to run PSSpred on an input fasta file with multiple peptides in it. Thanks

peptide sequence psspred • 773 views
ADD COMMENTlink modified 2.8 years ago by Pierre Lindenbaum133k • written 2.8 years ago by Tom30
gravatar for Pierre Lindenbaum
2.8 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum133k wrote:

linearize and loop

awk '/^>/ {printf("%s%s\t",(N>0?"\n":""),$0);N++;next;} {printf("%s",$0);} END {printf("\n");}' jeter.fa | cut -c2- | while read -a T; do echo -e ">${T[0]}\n${T[1]}" > ${T[0]}.fa && perl ${T[0]}.fa > ${T[0]}.out ; done
ADD COMMENTlink written 2.8 years ago by Pierre Lindenbaum133k
Please log in to add an answer.


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