aboutsummaryrefslogtreecommitdiff
path: root/workflows/yamlfa2ttl/check_sequence.py
diff options
context:
space:
mode:
authorAndreaGuarracino2021-01-08 23:44:47 +0100
committerAndreaGuarracino2021-01-08 23:44:47 +0100
commit53a47a18a45270c3f037fabe9cc973c66bfe50bf (patch)
tree1fb71fae0d629127a5da6f9ad879fff259016b14 /workflows/yamlfa2ttl/check_sequence.py
parent2ddf72a4028475eb65dfeab153f3565e698d5338 (diff)
downloadbh20-seq-resource-53a47a18a45270c3f037fabe9cc973c66bfe50bf.tar.gz
bh20-seq-resource-53a47a18a45270c3f037fabe9cc973c66bfe50bf.tar.lz
bh20-seq-resource-53a47a18a45270c3f037fabe9cc973c66bfe50bf.zip
sample_id in the FASTA has to match the sample_id in the YAMLyamlfa2ttl
Diffstat (limited to 'workflows/yamlfa2ttl/check_sequence.py')
-rw-r--r--workflows/yamlfa2ttl/check_sequence.py4
1 files changed, 2 insertions, 2 deletions
diff --git a/workflows/yamlfa2ttl/check_sequence.py b/workflows/yamlfa2ttl/check_sequence.py
index f92bf6d..58a65b9 100644
--- a/workflows/yamlfa2ttl/check_sequence.py
+++ b/workflows/yamlfa2ttl/check_sequence.py
@@ -26,7 +26,7 @@ def read_single_fasta(path_fasta):
return header, ''.join(sequence)
-print("FASTA QC: checking similarity to the reference")
+print("FASTA QC: checking similarity to the reference", file=sys.stderr)
header, sequence = read_single_fasta(path_fasta)
@@ -46,7 +46,7 @@ with tempfile.NamedTemporaryFile() as tmp_fasta:
"minimap2", "-c", "-x", "asm20",
tmp_sars_cov_2_reference_fasta.name, tmp_fasta.name
]
- print(" ".join(cmd))
+ print(" ".join(cmd), file=sys.stderr)
result = subprocess.run(cmd, stdout=subprocess.PIPE)
result.check_returncode()