Mercurial > repos > greg > cfsan_snp_pipeline_call_consensus_snps

changeset 1:0fd0a757072c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Wed, 22 Nov 2023 22:12:50 +0000
parents 3275d4866650
children d2ac1f5c05ed
files cfsan_snp_pipeline_call_consensus_snps.xml
diffstat 1 files changed, 4 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/cfsan_snp_pipeline_call_consensus_snps.xml	Fri Oct 13 17:57:44 2023 +0000
+++ b/cfsan_snp_pipeline_call_consensus_snps.xml	Wed Nov 22 22:12:50 2023 +0000
@@ -5,7 +5,10 @@
     </macros>
     <expand macro="requirements"/>
     <command detect_errors="exit_code"><![CDATA[
-cfsan_snp_pipeline call_consensus -v 4 -l '$input_merged_snp_sites' --vcfFileName 'consensus.vcf' '$input_call_pileup' &&
+#set $identifier = $input_call_pileup.element_identifier
+mkdir -p ./$identifier/ &&
+cp '$input_call_pileup' ./'$identifier'/reads.all.pileup &&
+cfsan_snp_pipeline call_consensus -v 4 -l '$input_merged_snp_sites' --vcfFileName 'consensus.vcf' ./'$identifier'/reads.all.pileup 1>&2 &&
 mv ./consensus.fasta '$output_fasta' &&
 mv ./consensus.vcf '$output_vcf'
     ]]></command>