Repository 'featurecounts'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/featurecounts

Changeset 13:386220cf6877 (2018-05-19)
Previous changeset 12:b714f4620411 (2018-05-07) Next changeset 14:85aaf50ad9dc (2018-06-01)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/featurecounts commit 2bd06c2b43c295fb4cf172c4f156fed5475855a4
modified:
featurecounts.xml
test-data/output_1_full.tab
test-data/output_1_jcounts.tab
test-data/output_1_medium.tab
test-data/output_1_short.tab
test-data/output_1_summary.tab
test-data/output_feature_lengths.tab
removed:
test-data/output_1_jcounts_with_header.tab
test-data/output_1_short_with_header.tab
test-data/output_1_summary_with_header.tab
b
diff -r b714f4620411 -r 386220cf6877 featurecounts.xml
--- a/featurecounts.xml Mon May 07 06:21:58 2018 -0400
+++ b/featurecounts.xml Sat May 19 03:53:55 2018 -0400
b
@@ -1,4 +1,4 @@
-<tool id="featurecounts" name="featureCounts" version="1.6.0.4" profile="16.04">
+<tool id="featurecounts" name="featureCounts" version="1.6.0.5" profile="16.04">
     <description>Measure gene expression in RNA-Seq experiments from SAM or BAM files.</description>
     <requirements>
         <requirement type="package" version="1.6.0">subread</requirement>
@@ -81,16 +81,9 @@
 
         '${alignment}'
 
-        ## Removal of comment
-        && grep -v "^#" "output"
-
-        #if $format.value != "tabdel_short":
-            ## and remove column-header line
-            | tail -n+2
-        #else
-            ## update header
-            | sed --expression='s|${alignment}|${alignment.element_identifier}|g'
-        #end if
+        ## Remove comment and add sample name to header
+        && grep -v "^#" "output" 
+        | sed -e 's|${alignment}|${alignment.element_identifier}|g'
         > body.txt
         ## Set the right columns for the tabular formats
         #if $format.value == "tabdel_medium":
@@ -101,7 +94,7 @@
             && cut -f 6 body.txt > gene_lengths.txt
             && paste expression_matrix.txt gene_lengths.txt > expression_matrix.txt.bak
             && mv -f expression_matrix.txt.bak '${output_medium}'
-        #elif $format.value == "tabdel_short" or $format.value == "tabdel_short_noheader":
+        #elif $format.value == "tabdel_short":
             && cut -f 1,7 body.txt > '${output_short}'
         #else:
             && cp body.txt '${output_full}'
@@ -112,19 +105,10 @@
         #end if
 
         #if str($extended_parameters.exon_exon_junction_read_counting_enabled.count_exon_exon_junction_reads) == "-J":
-            #if $format.value != "tabdel_short":
-              && tail -n+2 'output.jcounts' > '${output_jcounts}'
-            #else:
-
-              && sed --expression='s|${alignment}|${alignment.element_identifier}|g' 'output.jcounts' > '${output_jcounts}'
-            #end if
+            && sed -e 's|${alignment}|${alignment.element_identifier}|g' 'output.jcounts' > '${output_jcounts}'
         #end if
 
-        #if $format.value != "tabdel_short":
-            && tail -n+2 'output.summary' > '${output_summary}'
-        #else:
-            && sed --expression='s|${alignment}|${alignment.element_identifier}|g' 'output.summary' > '${output_summary}'
-        #end if
+        && sed -e 's|${alignment}|${alignment.element_identifier}|g' 'output.summary' > '${output_summary}'
     ]]></command>
     <inputs>
         <param name="alignment"
@@ -170,8 +154,7 @@
                type="select"
                label="Output format"
                help="The output format will be tabular, select the preferred columns here">
-            <option value="tabdel_short_noheader" selected="true">Gene-ID "\t" read-count (DESeq2 IUC wrapper compatible)</option>
-            <option value="tabdel_short">Gene-ID "\t" read-count (MultiQC/edgeR/limma-voom compatible, includes header in output)</option>
+            <option value="tabdel_short">Gene-ID "\t" read-count (MultiQC/DESeq2/edgeR/limma-voom compatible)</option>
             <option value="tabdel_medium">Gene-ID "\t" read-count "\t" gene-length</option>
             <option value="tabdel_full">featureCounts 1.4.0+ default (includes regions provided by the GTF file)</option>
         </param>
@@ -428,7 +411,7 @@
         <data format="tabular"
               name="output_short"
               label="${tool.name} on ${on_string}">
-            <filter>format == "tabdel_short_noheader" or format == "tabdel_short"</filter>
+            <filter>format == "tabdel_short"</filter>
             <actions>
                 <action name="column_names" type="metadata" default="Geneid,${alignment.element_identifier}" />
             </actions>
@@ -470,23 +453,6 @@
         </data>
     </outputs>
     <tests>
-        <test expect_num_outputs="4">
-            <param name="alignment" value="featureCounts_input1.bam" ftype="bam" dbkey="hg38" />
-            <param name="anno_select" value="history"/>
-            <param name="reference_gene_sets" value="featureCounts_guide.gff" ftype="gff" dbkey="hg38" />
-            <param name="format" value="tabdel_short_noheader" />
-            <param name="include_feature_length_file" value="true"/>
-            <param name="count_exon_exon_junction_reads" value="-J"/>
-            <output name="output_short" file="output_1_short.tab">
-                <metadata name="column_names" value="Geneid,featureCounts_input1.bam"/>
-            </output>
-            <output name="output_summary" file="output_1_summary.tab">
-                <metadata name="column_names" value="Status,featureCounts_input1.bam"/>
-            </output>
-            <output name="output_jcounts" file="output_1_jcounts.tab">
-                <metadata name="column_names" value="PrimaryGene,SecondaryGene,Site1_chr,Site1_location,Site1_strand,Site2_chr,Site2_location,Site2_strand,featureCounts_input1.bam"/>
-            </output>
-        </test>
         <test expect_num_outputs="3">
             <param name="alignment" value="featureCounts_input1.bam" ftype="bam" dbkey="hg38" />
             <param name="anno_select" value="history"/>
@@ -523,13 +489,13 @@
             <param name="format" value="tabdel_short" />
             <param name="include_feature_length_file" value="true"/>
             <param name="count_exon_exon_junction_reads" value="-J"/>
-            <output name="output_short" file="output_1_short_with_header.tab">
+            <output name="output_short" file="output_1_short.tab">
                 <metadata name="column_names" value="Geneid,featureCounts_input1.bam"/>
             </output>
-            <output name="output_summary" file="output_1_summary_with_header.tab">
+            <output name="output_summary" file="output_1_summary.tab">
                 <metadata name="column_names" value="Status,featureCounts_input1.bam"/>
             </output>
-            <output name="output_jcounts" file="output_1_jcounts_with_header.tab">
+            <output name="output_jcounts" file="output_1_jcounts.tab">
                 <metadata name="column_names" value="PrimaryGene,SecondaryGene,Site1_chr,Site1_location,Site1_strand,Site2_chr,Site2_location,Site2_strand,featureCounts_input1.bam"/>
             </output>
         </test>
@@ -582,7 +548,7 @@
 
 Output format
 -------------
-FeatureCounts produces a table containing counted reads, per gene, per row. Optionally the last column can be set to be the effective gene-length. These tables are compatible with the DESeq2 Galaxy wrapper by IUC. Column names are added as metadata object.
+FeatureCounts produces a table containing counted reads, per gene, per row. Optionally the last column can be set to be the effective gene-length. These tables are compatible with the DESeq2, edgeR and limma-voom Galaxy wrappers by IUC.
 
 .. _Subread: http://bioinf.wehi.edu.au/subread-package/SubreadUsersGuide.pdf
     ]]></help>
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_full.tab
--- a/test-data/output_1_full.tab Mon May 07 06:21:58 2018 -0400
+++ b/test-data/output_1_full.tab Sat May 19 03:53:55 2018 -0400
b
b'@@ -1,2 +1,3 @@\n+Geneid\tChr\tStart\tEnd\tStrand\tLength\tfeatureCounts_input1.bam\n left\tX;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X;X\t30948;51312;64215;90012;92326;168410;212059;322007;445010;539893;577936;587701;598742;604149;615232;705855;721190;800289;915767;926940;957344;1407723;1714165;2025638;2293426;2319538;2969222;3072054;3833560;4051480;4524542;4581901;4656840;5374647;5418429;5620971;6228798;6286459;6649313;6917189;7018834;7596495;8186836;8206413;8258438;8815918;8842347;9401726;9469379;10031803;10242387;10418659;10444444;10989586;11134085;11162151;11465040;11499223;11637155;12057077;12061726;12639971;13295077;13323052;13357233;13924339;13998611;14197670;14483263;14612714;14664006;14830448;14916699;14927680;14970356;15220273;15306642;15329260;15375323;15730527;15734786;16435015;16441298;16829597;17125391;17450531;17686539;18250193;18308377;18323082;18448967;18647925;18706010;19030140;19396008;19559764;19588853;19926904;19990109;20004574;20037228;20073287;20081159;20088846;20096706;20129221;20235565;20238464;20247900;20316429;20374707;20397489;20446575;20707308;20872979;20904331;20907486;20970376;20996038;21046284;21061792;21125285;21153965;21171043;21291889;21322696;21339280;21344984;21376449;21455828;21497479;21685998;21860031;21874192;21886041;21887574;21895770;21907764;21914659;21946675;21966395;22003324;22021878;22030491;22034624;22036925;22048433;22057124;22082873;22139935;22195806;22223174;22239374;22243264;22276936;22283375;22289776;22339332;22385377;22397699;22406838;53883;77479;79405;82402;184346;186159;189437;200157;230841;558910;1461250;1521562;1630699;1712287;1824125;1848474;2202381;2503973;2557975;2717376;2738488;2951607;2963803;3112094;3308606;3357795;3386680;3556511;3647008;3679400;3689927;3735908;3880451;4178048;4593999;4625418;4683101;4685985;4721556;4989861;5232315;5283350;5767200;5820682;6065171;6323576;6338891;6456406;6499064;6908350;7374377;7677833;7739314;7925028;8353372;9692111;10160891;10339131;10579639;11077839;11205652;11428629;11527005;11799400;12195551;12275459;12824166;13831569;14051181;14309904;14445232;14476735;14530058;14660004;14671630;14709104;14712353;14941525;15167294;15315507;15326753;15448269;15449888;15468199;15469456;16116366;16234256;16386234;16951522;17007337;17088792;17141025;17155869;17174503;17324955;17787185;18024852;18087657;18213198;18513283;18677415;18850658;18954002;19001351;19015489;19085704;19224162;19250578;19265090;19273743;19330858;19365718;19489665;19567246;19603659;19690936;19703442;19766215;19819597;19829825;19854295;19892332;19986996;19993521;20033920;20214636;20267464;20331167;20351467;20485359;20503414;20676657;20727409;20750652;20765308;20811305;20874984;20876509;20932446;20956729;21019641;21047881;21051404;21086568;21180050;21185984;21304742;21329282;21457944;21474559;21494393;21877985;21889900;21891428;21909599;21912243;21974809;21987500;21993503;21997723;22001309;22039357;22043411;22095242;22131802;22143032;22145369;22158980;22177427;22181659;22185911;22190019;22234403;22250617;22256814;22271776;22326110;22336396;22412656\t31447;51811;64714;90511;92825;168909;212558;322506;445509;540392;578435;588200;599241;604648;615731;706354;721689;800788;916266;927439;957843;1408222;1714664;2026137;2293925;2320037;2969721;3072553;3834059;4051979;4525041;4582400;4657339;5375146;5418928;5621470;6229297;6286958;6649812;6917688;7019333;7596994;8187335;820'..b'813211;20875962;20877360;20933374;20958184;21020880;21049038;21052363;21087529;21181012;21186897;21305490;21335152;21458705;21475408;21495636;21879005;21890892;21892197;21910516;21913162;21976502;21988186;21994450;21998648;22002053;22040352;22044257;22096103;22133042;22144129;22146581;22159740;22178135;22182560;22186943;22191085;22235245;22251644;22257959;22272465;22327885;22337937;22413523\t32671;52817;65834;91417;93737;169862;218087;329133;453514;541300;579251;588991;600099;605652;616640;712428;729923;803022;917038;935515;967079;1408985;1715486;2028083;2302720;2320967;2975470;3080661;3834833;4059457;4525816;4583161;4658270;5379533;5423234;5626691;6232251;6292155;6650736;6925521;7028934;7597873;8190816;8207704;8259768;8817277;8850798;9404164;9470645;10040243;10250896;10426654;10446849;10997630;11135498;11164585;11466225;11500610;11645243;12058339;12068095;12641231;13296368;13325139;13363752;13929867;14008703;14204823;14484448;14614203;14665260;14840646;14918203;14936149;14971898;15221460;15311867;15330635;15376659;15731752;15736152;16436510;16442675;16830936;17126752;17451975;17690779;18251655;18309757;18329335;18451430;18652532;18707408;19038001;19397505;19561379;19591476;19931466;19991289;20006033;20038670;20075997;20083884;20091556;20099431;20130616;20236911;20239707;20250305;20317717;20377169;20398964;20448055;20715852;20874459;20906025;20908816;20971789;20997396;21047529;21063333;21126812;21155238;21172291;21293390;21324123;21340490;21346294;21378082;21457385;21499776;21687843;21862493;21875822;21887403;21889066;21897335;21909292;21915944;21948116;21972417;22004848;22029088;22031874;22036082;22038719;22049778;22058484;22091278;22141384;22197424;22226496;22240888;22245836;22278264;22285879;22291074;22340730;22386633;22405874;22409269;55779;78840;80653;83663;185622;187527;198008;201435;232070;560290;1463976;1529972;1640805;1713541;1832639;1857356;2211907;2512032;2560060;2726865;2746840;2952792;2967873;3122184;3315732;3360140;3396780;3557708;3655505;3680588;3691342;3737093;3889026;4184192;4595265;4631435;4685294;4687334;4729869;4998370;5234487;5287396;5768399;5826828;6073546;6329745;6345591;6465831;6503028;6909598;7376809;7679065;7741746;7930262;8354975;9700105;10170982;10340355;10585717;11083857;11213940;11436041;11535017;11800858;12197791;12280303;12832222;13835774;14057780;14315601;14453835;14477969;14536428;14661479;14677350;14710306;14718498;14943134;15169400;15317962;15328173;15449537;15451143;15469435;15470709;16125657;16244354;16392753;16952707;17013907;17094820;17142218;17157084;17175693;17326222;17790913;18026536;18089130;18214480;18514605;18680321;18860262;18959719;19003520;19016938;19086895;19225376;19252956;19266350;19274928;19339318;19367219;19491202;19568605;19606389;19699191;19704787;19767602;19820837;19831333;19855478;19893616;19988426;19997525;20035191;20215900;20269568;20332421;20352774;20486718;20504666;20677999;20728615;20752145;20766572;20813710;20876461;20877859;20933873;20958683;21021379;21049537;21052862;21088028;21181511;21187396;21305989;21335651;21459204;21475907;21496135;21879504;21891391;21892696;21911015;21913661;21977001;21988685;21994949;21999147;22002552;22040851;22044756;22096602;22133541;22144628;22147080;22160239;22178634;22183059;22187442;22191584;22235744;22252143;22258458;22272964;22328384;22338436;22414022\t+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;+;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-;-\t170000\t66\n'
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_jcounts.tab
--- a/test-data/output_1_jcounts.tab Mon May 07 06:21:58 2018 -0400
+++ b/test-data/output_1_jcounts.tab Sat May 19 03:53:55 2018 -0400
b
@@ -0,0 +1,1 @@
+PrimaryGene SecondaryGenes Site1_chr Site1_location Site1_strand Site2_chr Site2_location Site2_strand featureCounts_input1.bam
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_jcounts_with_header.tab
--- a/test-data/output_1_jcounts_with_header.tab Mon May 07 06:21:58 2018 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,1 +0,0 @@
-PrimaryGene SecondaryGenes Site1_chr Site1_location Site1_strand Site2_chr Site2_location Site2_strand featureCounts_input1.bam
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_medium.tab
--- a/test-data/output_1_medium.tab Mon May 07 06:21:58 2018 -0400
+++ b/test-data/output_1_medium.tab Sat May 19 03:53:55 2018 -0400
b
@@ -1,2 +1,3 @@
+Geneid featureCounts_input1.bam Length
 left 92 170000
 right 66 170000
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_short.tab
--- a/test-data/output_1_short.tab Mon May 07 06:21:58 2018 -0400
+++ b/test-data/output_1_short.tab Sat May 19 03:53:55 2018 -0400
b
@@ -1,2 +1,3 @@
+Geneid featureCounts_input1.bam
 left 92
 right 66
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_short_with_header.tab
--- a/test-data/output_1_short_with_header.tab Mon May 07 06:21:58 2018 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,3 +0,0 @@
-Geneid featureCounts_input1.bam
-left 92
-right 66
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_summary.tab
--- a/test-data/output_1_summary.tab Mon May 07 06:21:58 2018 -0400
+++ b/test-data/output_1_summary.tab Sat May 19 03:53:55 2018 -0400
b
@@ -1,3 +1,4 @@
+Status featureCounts_input1.bam
 Assigned 158
 Unassigned_Unmapped 0
 Unassigned_MappingQuality 0
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_1_summary_with_header.tab
--- a/test-data/output_1_summary_with_header.tab Mon May 07 06:21:58 2018 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,13 +0,0 @@
-Status featureCounts_input1.bam
-Assigned 158
-Unassigned_Unmapped 0
-Unassigned_MappingQuality 0
-Unassigned_Chimera 0
-Unassigned_FragmentLength 0
-Unassigned_Duplicate 0
-Unassigned_MultiMapping 0
-Unassigned_Secondary 0
-Unassigned_Nonjunction 0
-Unassigned_NoFeatures 6078
-Unassigned_Overlapping_Length 0
-Unassigned_Ambiguity 0
b
diff -r b714f4620411 -r 386220cf6877 test-data/output_feature_lengths.tab
--- a/test-data/output_feature_lengths.tab Mon May 07 06:21:58 2018 -0400
+++ b/test-data/output_feature_lengths.tab Sat May 19 03:53:55 2018 -0400
b
@@ -1,2 +1,3 @@
+Geneid Length
 left 170000
 right 170000