sp10 + t 60 et 80

Snakefile

@@ -31,7 +31,7 @@ READS_ILLUMINA_BMB="ERR1308675.fastq"
 REF_CCN="ccnscafold.fa"
 REF_BMB="bmbscafold.fa"
 ID=["CCN_yeast","BMB_yeast","diploid_yeast","ecoli_hifi","ecoli_illumina","human"]
-CONFIGS=["","_clustal","_tcoffee","_sp100"]
+CONFIGS=["_sp10","_clustal","_tcoffee","_sp100"]
 
 READS_HUMAN="reads.fastq"
 REF_HUMAN="chm13.draft_v1.1.fasta"
@@ -859,23 +859,26 @@ rule configuration_file_simulated_data:
 #-------------------------------------------------------------------------------
 # Configuration files
 #-------------------------------------------------------------------------------
-rule configuration_files:
+rule configuration_files_sp10:
     input :
         os.path.join('before_configuration_files','config_{id}.yaml')
     output :
-        os.path.join('configuration_files','config_{id}.yaml')
+        os.path.join('configuration_files','config_sp10_{id}.yaml')
     message :
-        "Configuration files ({wildcards.id})"
+        "Configuration files for sp10 ({wildcards.id})"
     log:
-        os.path.join('logs','configuration_files','1_configuration_files_{id}.log')
+        os.path.join('logs','configuration_files','1_configuration_files_for_sp10_{id}.log')
     shell:
-        'cp {input} {output};'
+        'set +o pipefail;'
+        'head -1 {input}| sed "s/$/_sp10/" >{output};'
+        'tail -2 {input} >>{output};'
         'echo -e "D: [10,20,30,45,50,60,100,150,200]" >>{output};'
         'echo -e "S: [100,200,500,1000,2000,5000,10000]" >>{output};'
-        'echo -e "T: [50,70]" >>{output};'
+        'echo -e "T: [50,60,70,80]" >>{output};'
         'echo -e "M: [muscle,mafft,poa,spoa,abpoa,kalign,kalign3]" >>{output};'
         'echo -e "O: 10" >>{output}'
         
+
 rule configuration_files_for_clustal:
     input :
         os.path.join('before_configuration_files','config_{id}.yaml')
@@ -891,7 +894,7 @@ rule configuration_files_for_clustal:
         'tail -2 {input} >>{output};'
         'echo -e "D: [10,20,30,45,50,60,100,150,200]" >>{output};'
         'echo -e "S: [100,200,500,1000,2000,5000,10000]" >>{output};'
-        'echo -e "T: [50,70]" >>{output};'
+        'echo -e "T: [50,60,70,80]" >>{output};'
         'echo -e "M: [clustalo]" >>{output};'
         'echo -e "O: 10" >>{output}'
 
@@ -911,13 +914,13 @@ rule configuration_files_for_tcoffee:
         'tail -2 {input.config} >>{output};'
         'echo -e "D: [10,20,30,45,50,60,100]" >>{output};'
         'echo -e "S: [100,200,500,1000]" >>{output};'
-        'echo -e "T: [50,70]" >>{output};'
+        'echo -e "T: [50,60,70,80]" >>{output};'
         'echo -e "M: [tcoffee]" >>{output};'
         'echo -e "B: [`cat {input.start}`]" >>{output}'
         
 rule configuration_files_for_sp100:
     input :
-        config=os.path.join('before_configuration_files','config_{id}.yaml')
+        os.path.join('before_configuration_files','config_{id}.yaml'),
     output :
         os.path.join('configuration_files','config_sp100_{id}.yaml')
     message :
@@ -926,10 +929,10 @@ rule configuration_files_for_sp100:
         os.path.join('logs','configuration_files','3_configuration_files_for_sp100_{id}.log')
     shell:
         'set +o pipefail;'
-        'head -1 {input.config}| sed "s/$/_sp100/" >{output};'
-        'tail -2 {input.config} >>{output};'
+        'head -1 {input}| sed "s/$/_sp100/" >{output};'
+        'tail -2 {input} >>{output};'
         'echo -e "D: [10,20,30,45,50,60,100,150,200]" >>{output};'
         'echo -e "S: [500]" >>{output};'
-        'echo -e "T: [50,70]" >>{output};'
+        'echo -e "T: [50,60,70,80]" >>{output};'
         'echo -e "M: [muscle,mafft,poa,spoa,abpoa,kalign,kalign3]" >>{output};'
         'echo -e "O: 100" >>{output}'