Commit b49b68db authored by Mikael Salson's avatar Mikael Salson

repseq_vdj.py: Change header separator

The # is a comment in should-vdj.fa and it may be useful to see it in the output
(because we may put there a sequence ID)
parent c6492a75
Pipeline #100423 failed with stages
in 3 minutes and 7 seconds
......@@ -301,8 +301,8 @@ def header_igblast_results(ff_fasta, ff_igblast):
### Vidjil
VIDJIL_FINE = '{directory}/vidjil-algo --header-sep "#" -c designations -1 -2 -3 -d -g {directory}/germline/homo-sapiens.g %s >> %s'
VIDJIL_KMER = '{directory}/vidjil-algo -w 20 --header-sep "#" -b out -c windows -uuuU -1 -2 -g {directory}/germline/homo-sapiens.g %s > /dev/null ; cat out/out.segmented.vdj.fa out/out.unsegmented.vdj.fa >> %s'
VIDJIL_FINE = '{directory}/vidjil-algo --header-sep "|" -c designations -1 -2 -3 -d -g {directory}/germline/homo-sapiens.g %s >> %s'
VIDJIL_KMER = '{directory}/vidjil-algo -w 20 --header-sep "|" -b out -c windows -uuuU -1 -2 -g {directory}/germline/homo-sapiens.g %s > /dev/null ; cat out/out.segmented.vdj.fa out/out.unsegmented.vdj.fa >> %s'
def should_results_from_vidjil_output(f_log):
'''
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment