./buildali_uniref50_dssp.pl  /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26/query.fasta  /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26 
$ ./reformat.pl -M first -noss a3m a3m /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26/query.fasta /tmp/phyre2/23838/query.in.a3m -> 1
$ ./reformat.pl -r a3m psi /tmp/phyre2/23838/query.in.a3m /tmp/phyre2/23838/query.psi -> 1

************************************************
 Building alignment for /tmp/phyre2/23838/query (batch____)
************************************************
Building alignment for query with PSI-BLAST ...
$ /bmm/www/servers/phyre2/bin/blast/blastpgp -I T -s T -b 20000 -v 1 -e 0.001 -d /data/data/blast/uniref50.fasta -i /tmp/phyre2/23838/query.seq &> /tmp/phyre2/23838/query.bla
$ ./alignhits.pl -cov 80 -e 0.001  -bl 0.33 -bs 0.67  -best -psi -q /tmp/phyre2/23838/query.seq /tmp/phyre2/23838/query.bla /tmp/phyre2/23838/query.core.psi -> 4
$ ./alignhits.pl -cov 80 -e 0.001  -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/23838/query.core.psi  -best -psi -q /tmp/phyre2/23838/query.seq /tmp/phyre2/23838/query.bla /tmp/phyre2/23838/query.psi -> 6
$ ./alignhits.pl -cov 80 -e 0.001  -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/23838/query.core.psi  -best -a3m -q /tmp/phyre2/23838/query.seq /tmp/phyre2/23838/query.bla /tmp/phyre2/23838/query.a3m -> 6
Found 6 sequences in PSI-BLAST round 1 (E-value<1E-03, coverage>20%, db=uniref50.fasta)

$ /bmm/www/servers/phyre2/bin/blast/blastpgp -I T -s T -b 20000 -v 1 -e 0.001 -d /data/data/blast/uniref50.fasta -i /tmp/phyre2/23838/query.seq -B /tmp/phyre2/23838/query.psi &> /tmp/phyre2/23838/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/23838/query.core.psi   -a3m -q /tmp/phyre2/23838/query.seq /tmp/phyre2/23838/query.bla /tmp/phyre2/23838/query.tmp -> 6
$ cat /tmp/phyre2/23838/query.tmp >> /tmp/phyre2/23838/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/23838/query.core.psi   -psi -q /tmp/phyre2/23838/query.seq /tmp/phyre2/23838/query.bla /tmp/phyre2/23838/query.tmp -> 6
$ cat /tmp/phyre2/23838/query.tmp >> /tmp/phyre2/23838/query.psi
Found 6 sequences in PSI-BLAST round 2 (E-value<1E-03, coverage>20%, db=uniref50.fasta)


Extending query sequence ...
$ /bmm/www/servers/phyre2/bin/blast/blastpgp -I T -s T -I T -b 10 -v 10 -e 1e-6 -A 10 -f 15 -d /data/data/blast/uniref50.fasta -i /tmp/phyre2/23838/query.seq &> /tmp/phyre2/23838/query.bla
/bmm/www/servers/phyre2/bin/blast//fastacmd -d /data/data/blast/uniref50.fasta -s 'UniRef50_Q31US3'
Extended residues=0  original residues=97
$ ./hhfilter -M a3m -id 90  -cov 20 -diff 0 -i /tmp/phyre2/23838/query.a3m -o /tmp/phyre2/23838/query.fil
Input file = /tmp/phyre2/23838/query.a3m
Output file = /tmp/phyre2/23838/query.fil
Read /tmp/phyre2/23838/query.a3m with 12 sequences
Alignment in /tmp/phyre2/23838/query.a3m contains 97 match states
5 out of 12 sequences passed filter (20% min coverage, 90% max pairwise sequence identity)
$ ./reformat.pl -r a3m psi /tmp/phyre2/23838/query.fil /tmp/phyre2/23838/query.psi -> 5
Read DSSP state sequence

Searching DSSP state assignments...
name=batch____  range=
WARNING: Cannot open /bmm/data/dssp//batch____.dssp: No such file or directory
WARNING Cannot open /pdbbatch____.ent: No such file or directory
Predicting secondary structure with PSIPRED ...
$ /bmm/www/servers/phyre2/bin/blast/blastpgp -I T -s T -b 1 -j 1 -h 0.001 -d /bmm/www/servers/phyre/public_html/Version_1.1/data/dummydb -i /tmp/phyre2/23838/query.seq -B /tmp/phyre2/23838/query.psi -C /tmp/phyre2/23838/query.chk &> /tmp/phyre2/23838/query.bla
$ /bmm/home/kelley/psipred25/bin/psipred /tmp/phyre2/23838/query.mtx /bmm/home/kelley/psipred25/data/weights.dat /bmm/home/kelley/psipred25/data/weights.dat2 /bmm/home/kelley/psipred25/data/weights.dat3 /bmm/home/kelley/psipred25/data/weights.dat4 > /tmp/phyre2/23838/query.ss
$ /bmm/home/kelley/psipred25/bin/psipass2 /bmm/home/kelley/psipred25/data/weights_p2.dat 1 0.98 1.09 /tmp/phyre2/23838/query.ss2 /tmp/phyre2/23838/query.ss > /tmp/phyre2/23838/query.horiz
$ cat /tmp/phyre2/23838/query.a3m /tmp/phyre2/23838/query.in.a3m /tmp/phyre2/23838/query.fil > /tmp/phyre2/23838/query.tmp.a3m
$ ./hhfilter -M a3m -id 90  -cov 20 -diff 0 -i /tmp/phyre2/23838/query.tmp.a3m -o /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26/query.a3m
Input file = /tmp/phyre2/23838/query.tmp.a3m
Output file = /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26/query.a3m
Read /tmp/phyre2/23838/query.tmp.a3m with 8 sequences
Alignment in /tmp/phyre2/23838/query.tmp.a3m contains 97 match states
5 out of 6 sequences passed filter (20% min coverage, 90% max pairwise sequence identity)
Written alignment with 5 sequences to /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26/query.a3m

Finished  buildali.pl /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26/query.fasta /bmm/www/servers/phyre2/public_html/phyre2_output/50b0e543e40ddc26
