Entering edit mode
2.3 years ago
Xylanaser
▴
80
hey i wrote simple script that checks number of hommologs and organism in tax but i have problem with orgnism number... its often 0 where number of homologues are for example 1000
part of my code:
os.system('blastp -query %s -db /home/xyl/nr/nr -seg no -outfmt 7
-taxidlist /home/xyl/blast_tax/taxids/eukarya.ids -evalue 1e-3 -out /home/xyl/blast_tax/results/%s -max_target_seqs 20000 -num_threads 8' % (file, fileout3))
os.system('blastp -query %s -db /home/xyl/nr/nr -seg no -outfmt 18 -taxidlist /home/xyl/blast_tax/taxids/eukarya.ids -evalue 1e-3 -out /home/xyl/blast_tax/results/%s -max_target_seqs 20000 -num_threads 8' % (file, fileouttax3))
i use taxidlist for example: https://mega.nz/file/atdRXa6B#77ntvE3wObEnzGi6peLrUW5HnU611K74ufsgvSb3Eto
Does
nr
refer to preformatted blast database downloaded from NCBI?yes, i use newest blast and db
Have you downloaded
taxdb
files from NCBI and are they in the same folder as thenr
db?