[Biopython] NCBIWWW and megablast
Bart
bratdaking at gmail.com
Fri Oct 8 12:00:53 UTC 2010
Hey,
I was wondering why the megablast option (the greedy extension) in the
qblast is left out in the NCBIWWW.py?
I want to map a sequence to the human genome, and to mimic the NCBI website
I need a gapcost setting of "0 0", with the megablast option set to True.
The fix was to add the following line
('LCASE_MASK',lcase_mask),
('MEGABLAST',megablast),
('MATRIX_NAME',matrix_name),
to the parameters list of the qblast def and add:
megablast=None,
to the arguments.
But is there a reason this setting has been left out (it is as far as I can
see the only setting from the NCBI api missing)?
Cheers,
Bart
More information about the Biopython
mailing list