# HG changeset patch
# User bcclaywell
# Date 1427315950 14400
# Node ID 221a97bf686d5f512b20a20903689a55c725b007
# Parent ce6db18f5fd321fa6dbf887add82c7030bb7af92
Uploaded
diff -r ce6db18f5fd3 -r 221a97bf686d alignment-wrapper.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/alignment-wrapper.sh Wed Mar 25 16:39:10 2015 -0400
@@ -0,0 +1,14 @@
+#!/bin/bash
+
+source $1
+
+# adapted from yapp/bin/refpkg_align
+ref_sto=$(taxit rp ${REFPKG} aln_sto)
+profile=$(taxit rp ${REFPKG} profile)
+
+sto=$(mktemp -u).sto
+
+cmalign --cpu ${GALAXY_SLOTS:-4} -o "$sto" --sfile "${ALIGNED_SCORES}" --noprob --dnaout "$profile" "${INPUT_SEQS}" | grep -E '^#'
+
+esl-alimerge --dna --outformat afa "$ref_sto" "$sto" | \
+ seqmagick convert --output-format fasta --dash-gap - "${ALIGNED_SEQS}"
diff -r ce6db18f5fd3 -r 221a97bf686d alignment.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/alignment.xml Wed Mar 25 16:39:10 2015 -0400
@@ -0,0 +1,47 @@
+
+ with sequences in a reference package
+
+ yapp_env
+
+
+ macros.xml
+
+ echo "guppy $(guppy --version)"
+
+ alignment-wrapper.sh ${config}
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+REFPKG="${refpkg.extra_files_path}"
+INPUT_SEQS="${input_seqs}"
+
+ALIGNED_SEQS="${aligned_seqs}"
+ALIGNED_SCORES="${aligned_scores}"
+
+
+
+
+
+.. class:: infomark
+
+**What it does**
+
+This tool aligns query sequences with the reference sequences used to make the
+reference tree contained in the reference package and then merges the query and
+reference sequences.
+
+
+
diff -r ce6db18f5fd3 -r 221a97bf686d preprocessing-wrapper.sh
--- a/preprocessing-wrapper.sh Thu Feb 26 19:31:20 2015 -0500
+++ b/preprocessing-wrapper.sh Wed Mar 25 16:39:10 2015 -0400
@@ -7,14 +7,3 @@
--deduplicated-sequences-file ${DEDUP_INFO} \
${INPUT_SEQS} \
${DEDUP_SEQS}
-
-# adapted from yapp/bin/refpkg_align
-ref_sto=$(taxit rp ${REFPKG} aln_sto)
-profile=$(taxit rp ${REFPKG} profile)
-
-sto=$(mktemp -u).sto
-
-cmalign --cpu ${GALAXY_SLOTS:-4} -o "$sto" --sfile "${ALIGNED_SCORES}" --noprob --dnaout "$profile" "${DEDUP_SEQS}" | grep -E '^#'
-
-esl-alimerge --dna --outformat afa "$ref_sto" "$sto" | \
- seqmagick convert --output-format fasta --dash-gap - "${ALIGNED_SEQS}"
diff -r ce6db18f5fd3 -r 221a97bf686d preprocessing.xml
--- a/preprocessing.xml Thu Feb 26 19:31:20 2015 -0500
+++ b/preprocessing.xml Wed Mar 25 16:39:10 2015 -0400
@@ -1,5 +1,5 @@
-
- in preparation for phylogenetic placement
+
+ in preparation for alignment and placement
yapp_env
@@ -14,26 +14,20 @@
-
-
-
-REFPKG="${refpkg.extra_files_path}"
INPUT_SEQS="${input_seqs}"
SPLIT_MAP="${split_map}"
DEDUP_SEQS="${dedup_seqs}"
DEDUP_INFO="${dedup_info}"
-ALIGNED_SEQS="${aligned_seqs}"
-ALIGNED_SCORES="${aligned_scores}"