Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

addss.pl calls blastpgp and makemat which are legacy blast commands #374

Open
Beigephage opened this issue May 24, 2024 · 3 comments
Open

Comments

@Beigephage
Copy link

addss.pl runs blastpgp -b 1 -j 1 -h 0.001 -d $dummydb -i $infile -B $tmpfile.in.psi -C $tmpfile.chk 1> $tmpfile.blalog 2> $tmpfile.blalog and &HHPaths::System("$ncbidir/makemat -P $tmpfile"); but both commands are part of legacy blast and not the more recent BLAST+. How do you recommend updating addss.pl to allow this to still run?

Thank you

@milot-mirdita
Copy link
Member

We don't intend it update addss.pl currently. I'd be happy to review a pull request though.

@Beigephage
Copy link
Author

Hi Milot, how do you suggest running addss.pl without the legacy blast, which is not downloadable? Any help would be appreciated.

@Beigephage
Copy link
Author

commented out #"$ncbidir/blastpgp -b 1 -j 1 -h 0.001 -d $dummydb -i $infile -B $tmpfile.in.psi -C $tmpfile.chk 1> $tmpfile.blalog 2> $tmpfile.blalog"

to replace with

"$ncbidir/psiblast -db $dummydb -query $infile -inclusion_ethresh 0.001 -out_pssm $tmpfile.chk -num_iterations 3 -num_alignments 0 1> $tmpfile.blalog 2> $tmpfile.blalog"

and commented out
#&HHPaths::System("$ncbidir/makemat -P $tmpfile");
but unfortunately it is still incapable of adding secondary structure to the alignment.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants