Entering edit mode
2.4 years ago
Xylanaser
▴
80
hey i need to run BLAST+ search limited to specific superkingdoms (Bacteria, eukarya, viruses, archaea) but tax ids does not work ;( Do you know what the problem?
os.system('blastp -query input -db /home/xyl/nr/nr -seg no -outfmt 18 -taxids 2 -evalue 1e-3 -out output -max_target_seqs 50000 -num_threads 8' % (file, fileouttax1))