Mercurial > repos > iuc > stacks2_shortreads

changeset 1:1974fee35ca7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/stacks2 commit c300b84615660459bb0020fa74ccd3b874d329a4"
author iuc
date Mon, 30 Sep 2019 14:16:13 -0400
parents ad7a60726fc3
children 43e3eeb2e0ec
files macros.xml macros_process.xml stacks_shortreads.xml test-data/procrad/process_radtags2.out
diffstat 4 files changed, 100 insertions(+), 10 deletions(-) [+]
line wrap: on
line diff
--- a/macros.xml	Mon Jul 01 11:02:39 2019 -0400
+++ b/macros.xml	Mon Sep 30 14:16:13 2019 -0400
@@ -8,7 +8,7 @@
     </xml>
 
     <token name="@STACKS_VERSION@">2.4</token>
-    <token name="@WRAPPER_VERSION@">0</token>
+    <token name="@WRAPPER_VERSION@">1</token>
     <!-- fix to 18.01 since https://github.com/galaxyproject/galaxy/pull/7032 -->
     <token name="@PROFILE@">18.01</token>
 
@@ -216,8 +216,8 @@
         </expand>
     </xml>
 
-    <!-- fastq input (used in denovomap, tsv2bam, ustacks) 
-         - fastq_optional: makes fastq input optional (true/false) 
+    <!-- fastq input (used in denovomap, tsv2bam, ustacks)
+         - fastq_optional: makes fastq input optional (true/false)
          - se_option: wording for "single end" option (for tsv2bam this is the
               reverse reads for the others its the forward reads)
          - help: help text -->
@@ -568,10 +568,10 @@
         </expand>
     </xml>
 
-	<!-- variant calling option for use in gstacks and denovomap 
+	<!-- variant calling option for use in gstacks and denovomap
          default for var_alpha is 0.01 if model == marukilow (which is the case in denovomap and refmap)
          otherwise no default is is available and gstacks will output and error
-         "Error: No value was provided for \-\-var-alpha and there is no default for this model)" 
+         "Error: No value was provided for \-\-var-alpha and there is no default for this model)"
 	-->
     <xml name="variant_calling_options_vg" token_varalpha_default="">
         <param argument="--var-alpha" name="var_alpha" type="float" value="@VARALPHA_DEFAULT@" min="0" label="Alpha threshold for discovering SNPs" help="Default is 0.01 if the marukilow model is used (which is the case in refmap and denovomap), otherwise no default value is available." />
--- a/macros_process.xml	Mon Jul 01 11:02:39 2019 -0400
+++ b/macros_process.xml	Mon Sep 30 14:16:13 2019 -0400
@@ -75,7 +75,7 @@
                 <param name="score" type="integer" value="10" min="0" max="40" argument="-s" label="Set the score limit. If the average score within the sliding window drops below this value, the read is discarded" />
                 <param name="remove" type="boolean" checked="false" truevalue="-c" falsevalue="" argument="-c" label="Clean data, remove any read with an uncalled base" />
                 <param name="discard" type="boolean" checked="false" truevalue="-q" falsevalue="" argument="-q" label="Discard reads with low quality scores"/>
-                <param argument="--filter_illumina" type="boolean" checked="false" truevalue="--filter_illumina" falsevalue="" label="Discard reads that have been marked by Illumina's chastity/purity filter as failing" />
+                <param argument="--filter-illumina" name="filter_illumina" type="boolean" checked="false" truevalue="--filter-illumina" falsevalue="" label="Discard reads that have been marked by Illumina's chastity/purity filter as failing" />
             </when>
             <when value="no">
                 <param argument="--len_limit" type="integer" value="" optional="true" label="Minimum sequence length" help="useful if your data has already been trimmed"/>
--- a/stacks_shortreads.xml	Mon Jul 01 11:02:39 2019 -0400
+++ b/stacks_shortreads.xml	Mon Sep 30 14:16:13 2019 -0400
@@ -12,6 +12,10 @@
 
 mkdir stacks_inputs stacks_outputs &&
 
+#if $output_log
+    ln -s '$output_log' stacks_outputs/process_shortreads.log &&
+#end if
+
 #set ($link_command, $inputype) = $fastq_input_nonbatch( $input_type.fqinputs, $input_type.input_type_select, "_R%d_0" )
 $link_command
 
@@ -29,9 +33,6 @@
 $options_advanced.mate_pair
 $options_advanced.no_overhang
 
-#if $output_log
-    && mv stacks_outputs/process_shortreads.log $output_log
-#end if
 @PROCESS_FASTQ_POSTPROC@
     ]]></command>
 
@@ -82,13 +83,14 @@
             </param>
             <param name="input_type|barcode_encoding" value="--inline_null"/>
             <param name="barcode" value="procrad/barcodes"/>
-            <param name="capture" value="true"/>
+            <param name="capture" value="-D"/>
             <param name="no_read_trimming" value="--no_read_trimming" />
             <param name="mate_pair" value="--mate-pair" />
             <param name="no_overhang" value="--no_overhang" />
             <param name="outype" value="gzfastq"/>
             <param name="add_log" value="yes" />
             <assert_command>
+                <has_text text="-D" />
                 <has_text text="--no_read_trimming" />
                 <has_text text="--mate-pair" />
                 <has_text text="--no_overhang" />
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/procrad/process_radtags2.out	Mon Sep 30 14:16:13 2019 -0400
@@ -0,0 +1,88 @@
+process_radtags v2.4, executed 2019-09-30 14:57:00 (zlib-1.2.11)
+process_radtags -p stacks_inputs/ -i gzfastq -b /tmp/tmpe_shit/files/d/e/d/dataset_ded45538-15ef-471f-9585-07ee09b6e6d4.dat --inline_null -y gzfastq -o stacks_outputs -w 0.1 -s 11 -c -q -D -e ecoRI
+File	Retained Reads	Low Quality	Barcode Not Found	RAD cutsite Not Found	Total
+R1.fastq.gz	2222	0	0	4778	7000
+R2.fastq.gz	0	0	7000	0	7000
+
+Total Sequences	14000
+Barcode Not Found	7000
+Low Quality	0
+RAD Cutsite Not Found	4778
+Retained Reads	2222
+
+Barcode	Filename	Total	NoRadTag	LowQuality	Retained
+ATGGGG	PopA_01	160	94	0	66
+GGGTAA	PopA_02	146	86	0	60
+AGGAAA	PopA_03	224	132	0	92
+TTTAAG	PopA_04	168	146	0	22
+GGTGTG	PopA_05	220	154	0	66
+TGATGT	PopA_06	214	147	0	67
+GGTTGT	PopA_07	182	137	0	45
+ATAAGT	PopA_08	232	164	0	68
+AAGATA	PopA_09	182	134	0	48
+TGTGAG	PopA_10	168	128	0	40
+ATAGTT	PopA_11	200	132	0	68
+GGAAGG	PopA_12	178	136	0	42
+TTTGTG	PopA_13	190	152	0	38
+TTAAAT	PopA_14	166	96	0	70
+AATAAG	PopA_15	152	74	0	78
+AAGAGG	PopA_16	174	112	0	62
+TAGTGT	PopA_17	242	146	0	96
+TGGAAG	PopA_18	148	98	0	50
+GGGTTG	PopA_19	110	64	0	46
+CATCAT	PopA_20	178	144	0	34
+GGAGAG	PopB_20	210	142	0	68
+GTTTTA	PopB_01	136	90	0	46
+TGATAA	PopB_02	166	108	0	58
+GTTGAT	PopB_03	148	74	0	74
+AGATTA	PopB_04	180	156	0	24
+GTATAG	PopB_05	136	106	0	30
+TTGGGA	PopB_06	166	112	0	54
+ATATAT	PopB_07	242	156	0	86
+GATGAG	PopB_08	208	152	0	56
+GGGAAT	PopB_09	198	126	0	72
+AGTAAT	PopB_10	164	120	0	44
+GGGATA	PopB_11	146	90	0	56
+GAGAAG	PopB_12	162	110	0	52
+AGTAGA	PopB_13	204	146	0	58
+AAGGAT	PopB_14	180	128	0	52
+AGGGTA	PopB_15	134	72	0	62
+TGTTTT	PopB_16	138	98	0	40
+ATGATG	PopB_17	142	100	0	42
+GAGTTA	PopB_18	160	112	0	48
+ATGTAG	PopB_19	146	104	0	42
+
+Sequences not recorded
+Barcode	Total
+CCGGTG	772
+CCGCGC	750
+CCGACA	743
+CCGGTC	723
+CCGCCT	695
+CCGTGG	675
+CCGGTT	667
+CCGATC	579
+CCGTGT	417
+CCGTAT	252
+CCGCGT	198
+CCGCGG	150
+CCGCAT	117
+CCGTTC	94
+CCGTAC	38
+CCGGGA	32
+CCGGAG	26
+CCGAAC	24
+CCGAGT	18
+CCGTGA	10
+CCGGGT	7
+CCGCTT	2
+CCGACG	2
+CCGATT	1
+CCGCCC	1
+CCGGCA	1
+CCGGCT	1
+CCGGGG	1
+CCGCCG	1
+CCGTGC	1
+CCGTAG	1
+CCGACT	1