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

************************************************
 Building alignment for /tmp/phyre2/3570/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/3570/query.seq &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 80 -e 0.001  -bl 0.33 -bs 0.67  -best -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.core.psi -> 62
$ ./alignhits.pl -cov 80 -e 0.001  -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi  -best -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.psi -> 155
$ ./alignhits.pl -cov 80 -e 0.001  -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi  -best -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.a3m -> 155
Found 155 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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 3078
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 3078
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
Found 3078 sequences in PSI-BLAST round 2 (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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 2336
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 2336
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
Found 2336 sequences in PSI-BLAST round 3 (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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1905
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1905
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
Found 1905 sequences in PSI-BLAST round 4 (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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1728
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1728
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
Found 1728 sequences in PSI-BLAST round 5 (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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1639
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1639
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
Found 1639 sequences in PSI-BLAST round 6 (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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1574
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1574
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
Found 1574 sequences in PSI-BLAST round 7 (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/3570/query.seq -B /tmp/phyre2/3570/query.psi &> /tmp/phyre2/3570/query.bla
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -a3m -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1468
Found 1468 sequences in PSI-BLAST round 8 (E-value<1E-03, coverage>20%, db=uniref50.fasta)

$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.a3m
$ ./alignhits.pl -cov 20 -e 0.001   -bl 0 -bs 0.167 -bg 20 -B /tmp/phyre2/3570/query.core.psi   -psi -q /tmp/phyre2/3570/query.seq /tmp/phyre2/3570/query.bla /tmp/phyre2/3570/query.tmp -> 1468
$ cat /tmp/phyre2/3570/query.tmp >> /tmp/phyre2/3570/query.psi
$ ./hhfilter -M a3m -id 90  -cov 20 -diff 0 -i /tmp/phyre2/3570/query.a3m -o /tmp/phyre2/3570/query.fil
Input file = /tmp/phyre2/3570/query.a3m
Output file = /tmp/phyre2/3570/query.fil
Read /tmp/phyre2/3570/query.a3m with 13883 sequences
Alignment in /tmp/phyre2/3570/query.a3m contains 302 match states
6519 out of 13883 sequences passed filter (20% min coverage, 90% max pairwise sequence identity)
$ ./reformat.pl -r a3m psi /tmp/phyre2/3570/query.fil /tmp/phyre2/3570/query.psi -> 6519
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/3570/query.seq -B /tmp/phyre2/3570/query.psi -C /tmp/phyre2/3570/query.chk &> /tmp/phyre2/3570/query.bla
$ /bmm/home/kelley/psipred25/bin/psipred /tmp/phyre2/3570/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/3570/query.ss
$ /bmm/home/kelley/psipred25/bin/psipass2 /bmm/home/kelley/psipred25/data/weights_p2.dat 1 0.98 1.09 /tmp/phyre2/3570/query.ss2 /tmp/phyre2/3570/query.ss > /tmp/phyre2/3570/query.horiz
$ cat /tmp/phyre2/3570/query.a3m /tmp/phyre2/3570/query.in.a3m /tmp/phyre2/3570/query.fil > /tmp/phyre2/3570/query.tmp.a3m
$ ./hhfilter -M a3m -id 90  -cov 20 -diff 0 -i /tmp/phyre2/3570/query.tmp.a3m -o /bmm/www/servers/phyre2/public_html/phyre2_output/b5b984b7bb33ab28/query.a3m
Input file = /tmp/phyre2/3570/query.tmp.a3m
Output file = /bmm/www/servers/phyre2/public_html/phyre2_output/b5b984b7bb33ab28/query.a3m
Read /tmp/phyre2/3570/query.tmp.a3m with 6522 sequences
Alignment in /tmp/phyre2/3570/query.tmp.a3m contains 302 match states
6519 out of 6520 sequences passed filter (20% min coverage, 90% max pairwise sequence identity)
Written alignment with 6519 sequences to /bmm/www/servers/phyre2/public_html/phyre2_output/b5b984b7bb33ab28/query.a3m

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