Hi! Is there someone who knows what the correct code is if I want to use megablast instead of blastn in: b_results = NCBIWWW.qblast('blastn', 'bta_genome/all_contig', f_record).read() Thank you! Cheers, Eirik