Skip to content
Snippets Groups Projects
Commit 61a668fd authored by christine.tranchant_ird.fr's avatar christine.tranchant_ird.fr
Browse files

modify vecsceen rule

parent b78264f2
No related branches found
No related tags found
No related merge requests found
...@@ -30,7 +30,7 @@ rule vecscreen: ...@@ -30,7 +30,7 @@ rule vecscreen:
from Bio import SeqIO from Bio import SeqIO
# Load NCBI-TOOLS module # Load NCBI-TOOLS module
cmd_load=f"module load {params.module}" #cmd_load=f"module load {params.module}"
nb_seq, nb_seq_no_hits, nb_seq_suspicious = 0, 0, 0 nb_seq, nb_seq_no_hits, nb_seq_suspicious = 0, 0, 0
...@@ -43,8 +43,9 @@ rule vecscreen: ...@@ -43,8 +43,9 @@ rule vecscreen:
tmp.write('>' + str(seq.description) + '\n') tmp.write('>' + str(seq.description) + '\n')
tmp.write(gs.format_60(str(seq.seq)) + '\n') tmp.write(gs.format_60(str(seq.seq)) + '\n')
cmd = f'module load {params.module} && vecscreen -db {input.univec_file} -query {tmp_file} -outfmt 1 -text_output' cmd = f'vecscreen -db {input.univec_file} -query {tmp_file} -outfmt 1 -text_output'
print(f"\t\t\ {seq.description} : vecscreen cmd : {cmd}") #cmd = f'module load {params.module} && vecscreen -db {input.univec_file} -query {tmp_file} -outfmt 1 -text_output'
#print(f"\t\t\ {seq.description} : vecscreen cmd : {cmd}")
try: try:
subprocess.run(cmd, shell=True, check=True) subprocess.run(cmd, shell=True, check=True)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment