# HG changeset patch # User artbio # Date 1553466180 14400 # Node ID 279fdd92a6150854ed1fd0e7e7b4593ab15702e0 # Parent e11f91575af63eb070d9cc941218dcf03576efbf planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sr_bowtie_dataset_annotation commit 1d30d780dc66069ae9b52c909f65d713156442b3 diff -r e11f91575af6 -r 279fdd92a615 sr_bowtie_dataset_annotation.xml --- a/sr_bowtie_dataset_annotation.xml Wed Mar 20 07:12:53 2019 -0400 +++ b/sr_bowtie_dataset_annotation.xml Sun Mar 24 18:23:00 2019 -0400 @@ -1,4 +1,4 @@ - + by iterative alignments with sRbowtie bowtie @@ -13,6 +13,11 @@ #else: #set index_path = $refGenomeSource1.index.fields.path #end if + + #for $i in $AdditionalQueries: + bowtie-build -f $i.ownFile $i.ownFile.name 1>/dev/null && + #end for + #set method_prefix = "-v %s -k 1 --best" % str($mismatches) #if $input[0].is_of_type('fasta'): #set format = "-f" @@ -38,14 +43,13 @@ #else: #set to_align = "matched.fa" #end if - bowtie-build -f $i.ownFile subgenome 1>/dev/null && touch tmp_class_matched.fa tmp_class_unmatched.fa && bowtie -p \${GALAXY_SLOTS:-4} $method_prefix --al tmp_class_matched.fa --un tmp_class_unmatched.fa --suppress 6,7,8 - subgenome $format '$to_align' > tabular_bowtie_output.tab && + $i.ownFile.name $format '$to_align' > tabular_bowtie_output.tab && class_aligned=\$(( \$(wc -l < tmp_class_matched.fa)/2)) && class_unaligned=\$(( \$(wc -l < tmp_class_unmatched.fa)/2)) && echo -e "$sample\t$i.ownFile.name\t\$class_aligned\t\${genome_aligned}" >> $output &&