示例#1
0
    def run_blast(self, unique_AA_sequences_fasta_path, unique_AA_sequences_names_dict):
        self.run.warning("You elected to use NCBI's blastp for amino acid sequence search. Running blastp will \
                          be significantly slower than DIAMOND (although, anvi'o developers are convinced that \
                          you *are* doing the right thing, so, kudos to you).")
        blast = BLAST(unique_AA_sequences_fasta_path, run=self.run, progress=self.progress,
                          num_threads=self.num_threads, overwrite_output_destinations=self.overwrite_output_destinations)

        blast.names_dict = unique_AA_sequences_names_dict
        blast.log_file_path = self.log_file_path
        blast.search_output_path = self.get_output_file_path('blast-search-results.txt')

        return blast.get_blast_results()
示例#2
0
    def run_blast(self, unique_AA_sequences_fasta_path,
                  unique_AA_sequences_names_dict):
        self.run.warning(
            "You elected to use NCBI's `blastp` for amino acid sequence search. Running blastp will "
            "be significantly slower than DIAMOND, but in some cases, slightly more sensitive. "
            "We are unsure about whether the slight increase in sensitivity may justify significant "
            "increase in run time, but you are the boss.",
            lc="cyan")

        blast = BLAST(
            unique_AA_sequences_fasta_path,
            run=self.run,
            progress=self.progress,
            num_threads=self.num_threads,
            overwrite_output_destinations=self.overwrite_output_destinations)

        blast.names_dict = unique_AA_sequences_names_dict
        blast.log_file_path = self.log_file_path
        blast.search_output_path = self.get_output_file_path(
            'blast-search-results.txt')

        return blast.get_blast_results()