1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54

use blul_core::domain::{
    dtos::blast_builder::BlastBuilder,
    entities::execute_blastn::{ExecuteBlastn, ExecutionResponse},
};
use mycelium_base::utils::errors::{execution_err, MappedErrors};
use subprocess::{Exec, Redirection};

#[derive(Debug)]
pub struct ExecuteBlastnProcRepository {}

impl ExecuteBlastn for ExecuteBlastnProcRepository {
    fn run(
        &self,
        query_sequences: String,
        blast_config: BlastBuilder,
    ) -> Result<ExecutionResponse, MappedErrors> {
        let blast_response = match Exec::cmd("blastn")
            .stdin(&*query_sequences)
            .arg("-subject")
            .arg(&blast_config.subject_reads)
            .arg("-outfmt")
            .arg(&blast_config.out_format.to_string())
            .arg("-max_target_seqs")
            .arg(&blast_config.max_target_seqs.to_string())
            .arg("-perc_identity")
            .arg(&blast_config.perc_identity.to_string())
            .arg("-qcov_hsp_perc")
            .arg(&blast_config.query_cov.to_string())
            .arg("-strand")
            .arg(&blast_config.strand.to_string())
            .arg("-evalue")
            .arg(&blast_config.e_value.to_string())
            .arg("-word_size")
            .arg(&blast_config.word_size.to_string())
            .stdout(Redirection::Pipe)
            .stderr(Redirection::Pipe)
            .capture()
        {
            Err(err) => {
                return execution_err(format!(
                    "Unexpected error detected on execute blast: {err}"
                ))
                .as_error()
            }
            Ok(res) => res,
        };

        if !blast_response.success() {
            return Ok(ExecutionResponse::Fail(blast_response.stderr_str()));
        }

        Ok(ExecutionResponse::Success(blast_response.stdout_str()))
    }
}