Question: How to run PSSPred on multple sequences
gravatar for Tom
9 months ago by
Tom20 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 • 242 views
ADD COMMENTlink modified 9 months ago by Pierre Lindenbaum116k • written 9 months ago by Tom20
gravatar for Pierre Lindenbaum
9 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum116k 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 9 months ago by Pierre Lindenbaum116k
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: 2336 users visited in the last hour