# HG changeset patch # User ieguinoa # Date 1570649901 14400 # Node ID 2f5e9c0fe3676074aa78280b7900e07cd41a9b45 "planemo upload for repository https://github.com/ieguinoa/tximport-galaxy-wrapper commit 2bb25471c1320fb1206afa2c4daf536b6d6e275f-dirty" diff -r 000000000000 -r 2f5e9c0fe367 README.md --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/README.md Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,2 @@ +## tximport-galaxy-wrapper +Wrapper for the package [tximport](https://bioconductor.org/packages/release/bioc/html/tximport.html) diff -r 000000000000 -r 2f5e9c0fe367 test-data/Araport11_subset.gff3 --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/Araport11_subset.gff3 Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,106 @@ +Chr1 Araport11 gene 3631 5899 . + . ID=AT1G01010;Alias=ANAC001,NAC domain containing protein 1;symbol=NAC001;tid=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 mRNA 3631 5899 . + . ID=AT1G01010.1;Parent=AT1G01010;Name=AT1G01010.1;gene_id=AT1G01010 +Chr1 Araport11 exon 3631 3913 . + . ID=AT1G01010.1:exon:1;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 five_prime_UTR 3631 3759 . + . ID=AT1G01010.1:five_prime_UTR;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 exon 3996 4276 . + . ID=AT1G01010.1:exon:2;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 exon 4486 4605 . + . ID=AT1G01010.1:exon:3;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 exon 4706 5095 . + . ID=AT1G01010.1:exon:4;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 exon 5174 5326 . + . ID=AT1G01010.1:exon:5;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 exon 5439 5899 . + . ID=AT1G01010.1:exon:6;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 five_prime_UTR 3759 5439 . + . ID=AT1G01010.1:five_prime_UTR;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 three_prime_UTR 5631 5899 . + . ID=AT1G01010.1:three_prime_UTR;Parent=AT1G01010.1;Name=AT1G01010;gene_id=AT1G01010 +Chr1 Araport11 gene 6788 9130 . - . ID=AT1G01020;symbol=ARV1;tid=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 mRNA 6788 9130 . - . ID=AT1G01020.5;Parent=AT1G01020;Name=AT1G01020.5;gene_id=AT1G01020 +Chr1 Araport11 exon 6788 7069 . - . ID=AT1G01020.5:exon:1;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6788 6914 . - . ID=AT1G01020.5:five_prime_UTR;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7157 7232 . - . ID=AT1G01020.5:exon:2;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7384 7450 . - . ID=AT1G01020.5:exon:3;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7564 7649 . - . ID=AT1G01020.5:exon:4;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7762 7835 . - . ID=AT1G01020.5:exon:5;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7942 7987 . - . ID=AT1G01020.5:exon:6;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8236 8325 . - . ID=AT1G01020.5:exon:7;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8417 8464 . - . ID=AT1G01020.5:exon:8;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6914 8417 . - . ID=AT1G01020.5:five_prime_UTR;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8420 8464 . - . ID=AT1G01020.5:three_prime_UTR;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8594 9130 . - . ID=AT1G01020.5:exon:9;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8594 9130 . - . ID=AT1G01020.5:three_prime_UTR;Parent=AT1G01020.5;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 mRNA 6788 9130 . - . ID=AT1G01020.4;Parent=AT1G01020;Name=AT1G01020.4;gene_id=AT1G01020 +Chr1 Araport11 exon 6788 7069 . - . ID=AT1G01020.4:exon:1;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6788 6914 . - . ID=AT1G01020.4:five_prime_UTR;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7157 7232 . - . ID=AT1G01020.4:exon:2;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7384 7450 . - . ID=AT1G01020.4:exon:3;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7564 7649 . - . ID=AT1G01020.4:exon:4;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7762 7835 . - . ID=AT1G01020.4:exon:5;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7942 7987 . - . ID=AT1G01020.4:exon:6;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8236 8464 . - . ID=AT1G01020.4:exon:7;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6914 8236 . - . ID=AT1G01020.4:five_prime_UTR;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8443 8464 . - . ID=AT1G01020.4:three_prime_UTR;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8594 9130 . - . ID=AT1G01020.4:exon:8;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8594 9130 . - . ID=AT1G01020.4:three_prime_UTR;Parent=AT1G01020.4;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 mRNA 6788 9130 . - . ID=AT1G01020.3;Parent=AT1G01020;Name=AT1G01020.3;gene_id=AT1G01020 +Chr1 Araport11 exon 6788 7069 . - . ID=AT1G01020.3:exon:1;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6788 6914 . - . ID=AT1G01020.3:five_prime_UTR;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7157 7232 . - . ID=AT1G01020.3:exon:2;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7384 7450 . - . ID=AT1G01020.3:exon:3;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7564 7649 . - . ID=AT1G01020.3:exon:4;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7762 7835 . - . ID=AT1G01020.3:exon:5;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7942 7987 . - . ID=AT1G01020.3:exon:6;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8236 8464 . - . ID=AT1G01020.3:exon:7;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6914 8236 . - . ID=AT1G01020.3:five_prime_UTR;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8443 8464 . - . ID=AT1G01020.3:three_prime_UTR;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8571 9130 . - . ID=AT1G01020.3:exon:8;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8571 9130 . - . ID=AT1G01020.3:three_prime_UTR;Parent=AT1G01020.3;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 mRNA 6788 9130 . - . ID=AT1G01020.1;Parent=AT1G01020;Name=AT1G01020.1;gene_id=AT1G01020 +Chr1 Araport11 exon 6788 7069 . - . ID=AT1G01020.1:exon:1;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6788 6914 . - . ID=AT1G01020.1:five_prime_UTR;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7157 7232 . - . ID=AT1G01020.1:exon:2;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7384 7450 . - . ID=AT1G01020.1:exon:3;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7564 7649 . - . ID=AT1G01020.1:exon:4;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7762 7835 . - . ID=AT1G01020.1:exon:5;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7942 7987 . - . ID=AT1G01020.1:exon:6;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8236 8325 . - . ID=AT1G01020.1:exon:7;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8417 8464 . - . ID=AT1G01020.1:exon:8;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8571 9130 . - . ID=AT1G01020.1:exon:9;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6914 8571 . - . ID=AT1G01020.1:five_prime_UTR;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8667 9130 . - . ID=AT1G01020.1:three_prime_UTR;Parent=AT1G01020.1;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 mRNA 6788 8737 . - . ID=AT1G01020.2;Parent=AT1G01020;Name=AT1G01020.2;gene_id=AT1G01020 +Chr1 Araport11 exon 6788 7069 . - . ID=AT1G01020.2:exon:1;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6788 7069 . - . ID=AT1G01020.2:five_prime_UTR;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7157 7450 . - . ID=AT1G01020.2:exon:2;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 7157 7314 . - . ID=AT1G01020.2:five_prime_UTR;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7564 7649 . - . ID=AT1G01020.2:exon:3;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7762 7835 . - . ID=AT1G01020.2:exon:4;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7942 7987 . - . ID=AT1G01020.2:exon:5;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8236 8325 . - . ID=AT1G01020.2:exon:6;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8417 8464 . - . ID=AT1G01020.2:exon:7;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8571 8737 . - . ID=AT1G01020.2:exon:8;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 7314 8571 . - . ID=AT1G01020.2:five_prime_UTR;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8667 8737 . - . ID=AT1G01020.2:three_prime_UTR;Parent=AT1G01020.2;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 mRNA 6788 8737 . - . ID=AT1G01020.6;Parent=AT1G01020;Name=AT1G01020.6;gene_id=AT1G01020 +Chr1 Araport11 exon 6788 7069 . - . ID=AT1G01020.6:exon:1;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 6788 7069 . - . ID=AT1G01020.6:five_prime_UTR;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7157 7450 . - . ID=AT1G01020.6:exon:2;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 7157 7314 . - . ID=AT1G01020.6:five_prime_UTR;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 7564 7649 . - . ID=AT1G01020.6:exon:3;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8236 8325 . - . ID=AT1G01020.6:exon:4;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8417 8464 . - . ID=AT1G01020.6:exon:5;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 five_prime_UTR 7314 8417 . - . ID=AT1G01020.6:five_prime_UTR;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8420 8464 . - . ID=AT1G01020.6:three_prime_UTR;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 exon 8594 8737 . - . ID=AT1G01020.6:exon:6;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 three_prime_UTR 8594 8737 . - . ID=AT1G01020.6:three_prime_UTR;Parent=AT1G01020.6;Name=AT1G01020;gene_id=AT1G01020 +Chr1 Araport11 gene 11649 13714 . - . ID=AT1G01030;symbol=NGA3;full_name=NGATHA3;tid=AT1G01030.1;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 mRNA 11649 13714 . - . ID=AT1G01030.1;Parent=AT1G01030;Name=AT1G01030.1;gene_id=AT1G01030 +Chr1 Araport11 exon 11649 13173 . - . ID=AT1G01030.1:exon:1;Parent=AT1G01030.1;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 five_prime_UTR 11649 11863 . - . ID=AT1G01030.1:five_prime_UTR;Parent=AT1G01030.1;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 three_prime_UTR 12941 13173 . - . ID=AT1G01030.1:three_prime_UTR;Parent=AT1G01030.1;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 exon 13335 13714 . - . ID=AT1G01030.1:exon:2;Parent=AT1G01030.1;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 three_prime_UTR 13335 13714 . - . ID=AT1G01030.1:three_prime_UTR;Parent=AT1G01030.1;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 mRNA 11649 13714 . - . ID=AT1G01030.2;Parent=AT1G01030;Name=AT1G01030.2;gene_id=AT1G01030 +Chr1 Araport11 exon 11649 12354 . - . ID=AT1G01030.2:exon:1;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 five_prime_UTR 11649 11863 . - . ID=AT1G01030.2:five_prime_UTR;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 exon 12424 13173 . - . ID=AT1G01030.2:exon:2;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 five_prime_UTR 11863 12424 . - . ID=AT1G01030.2:five_prime_UTR;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 three_prime_UTR 12941 13173 . - . ID=AT1G01030.2:three_prime_UTR;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 exon 13335 13714 . - . ID=AT1G01030.2:exon:3;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 +Chr1 Araport11 three_prime_UTR 13335 13714 . - . ID=AT1G01030.2:three_prime_UTR;Parent=AT1G01030.2;Name=AT1G01030;gene_id=AT1G01030 + diff -r 000000000000 -r 2f5e9c0fe367 test-data/cached_locally/tx2gene.loc --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/cached_locally/tx2gene.loc Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,5 @@ +#The tx2gene.loc file has this format: +# +# + +Ath_Araport11_subset Ath Arabidopsis thaliana - Araport11 - subset ${__HERE__}/../tx2gene.tab diff -r 000000000000 -r 2f5e9c0fe367 test-data/custom_sample.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/custom_sample.tab Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,14 @@ +Transcript_id_here Here_goes_the_length Extra_useless_column_1 Abundance_goes_here Extra_useless_column_2 Here_goes_the_counts +AT1G01010.1 1688 1487.57 4.6817 1 156 +AT1G01020.3 1420 1219.57 1.68167 2 45.9399 +AT1G01020.4 711 510.567 6.65521 3 76.113 +AT1G01020.1 738 537.567 18.1462 4 218.506 +AT1G01020.5 1179 978.567 0.491433 5 10.7721 +AT1G01020.6 617 416.567 0.607591 6 5.66944 +AT1G01020.2 647 446.567 2.63345e-08 7 2.63425e-07 +AT1G03987.1 272 97.3356 0 8 0 +AT1G01030.2 1456 1255.57 0.519509 9 14.6109 +AT1G01030.1 1905 1704.57 0.298284 10 11.3891 +AT1G01040.1 6276 6075.57 1.56643 11 213.178 +AT1G03993.1 788 587.567 0 12 0 +AT1G01040.2 5877 5676.57 4.30834 13 547.822 diff -r 000000000000 -r 2f5e9c0fe367 test-data/salmon_sample1.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/salmon_sample1.tab Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,14 @@ +Name Length EffectiveLength TPM NumReads +AT1G01010.1 1688 1487.57 4.6817 156 +AT1G01020.3 1420 1219.57 1.68167 45.9399 +AT1G01020.4 711 510.567 6.65521 76.113 +AT1G01020.1 738 537.567 18.1462 218.506 +AT1G01020.5 1179 978.567 0.491433 10.7721 +AT1G01020.6 617 416.567 0.607591 5.66944 +AT1G01020.2 647 446.567 2.63345e-08 2.63425e-07 +AT1G03987.1 272 97.3356 0 0 +AT1G01030.2 1456 1255.57 0.519509 14.6109 +AT1G01030.1 1905 1704.57 0.298284 11.3891 +AT1G01040.1 6276 6075.57 1.56643 213.178 +AT1G03993.1 788 587.567 0 0 +AT1G01040.2 5877 5676.57 4.30834 547.822 diff -r 000000000000 -r 2f5e9c0fe367 test-data/salmon_sample2.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/salmon_sample2.tab Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,14 @@ +Name Length EffectiveLength TPM NumReads +AT1G01010.1 1688 1487.57 4.6817 156 +AT1G01020.3 1420 1219.57 1.68167 45.9399 +AT1G01020.4 711 510.567 6.65521 76.113 +AT1G01020.1 738 537.567 18.1462 218.506 +AT1G01020.5 1179 978.567 0.491433 10.7721 +AT1G01020.6 617 416.567 0.607591 5.66944 +AT1G01020.2 647 446.567 2.63345e-08 2.63425e-07 +AT1G03987.1 272 97.3356 0 0 +AT1G01030.2 1456 1255.57 0.519509 14.6109 +AT1G01030.1 1905 1704.57 0.298284 11.3891 +AT1G01040.1 6276 6075.57 1.56643 213.178 +AT1G03993.1 788 587.567 0 0 +AT1G01040.2 5877 5676.57 4.30834 547.822 diff -r 000000000000 -r 2f5e9c0fe367 test-data/tx2gene.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/tx2gene.tab Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,14 @@ +Transcript Gene +AT1G01010.1 AT1G01010 +AT1G01020.3 AT1G01020 +AT1G01020.4 AT1G01020 +AT1G01020.1 AT1G01020 +AT1G01020.5 AT1G01020 +AT1G01020.6 AT1G01020 +AT1G01020.2 AT1G01020 +AT1G03987.1 AT1G03987 +AT1G01030.2 AT1G01030 +AT1G01030.1 AT1G01030 +AT1G01040.1 AT1G01040 +AT1G03993.1 AT1G03993 +AT1G01040.2 AT1G01040 diff -r 000000000000 -r 2f5e9c0fe367 tool-data/tx2gene.loc.sample --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/tx2gene.loc.sample Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,1 @@ +#value, dbkey, name, path diff -r 000000000000 -r 2f5e9c0fe367 tool_data_table_conf.xml.sample --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_data_table_conf.xml.sample Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,7 @@ + + + + value, dbkey, name, path + +
+
diff -r 000000000000 -r 2f5e9c0fe367 tool_data_table_conf.xml.test --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_data_table_conf.xml.test Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,8 @@ + + + + value, dbkey, name, path + +
+
+ diff -r 000000000000 -r 2f5e9c0fe367 tximport-galaxy-wrapper.tar Binary file tximport-galaxy-wrapper.tar has changed diff -r 000000000000 -r 2f5e9c0fe367 tximport.R --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tximport.R Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,119 @@ +# setup R error handling to go to stderr +options( show.error.messages=F, error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } ) + +# we need that to not crash galaxy with an UTF8 error on German LC settings. +loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") + +library("getopt") +#library("tools") +options(stringAsFactors = FALSE, useFancyQuotes = FALSE) +args <- commandArgs(trailingOnly = TRUE) + +# get options, using the spec as defined by the enclosed list. +# we read the options from the default: commandArgs(TRUE). +spec <- matrix(c( + "help", "h", 0, "logical", + "base_dir", "w", 1, "character", + "out_file", "o", 1, "character", + "countsFiles", "n", 1, "character", + "countsFromAbundance", "r", 1, "character", + "format", "v", 1, "character", + "gff_file", "H", 0, "character", + "tx2gene", "f", 0, "character", + "geneIdCol", "l", 0, "character", + "txIdCol" , "p", 1, "character", + "abundanceCol", "i", 0, "character", + "countsCol", "y", 1, "character", + "lengthCol", "x", 1, "character"), + byrow=TRUE, ncol=4) + +opt <- getopt(spec) + + + + +# if help was asked for print a friendly message +# and exit with a non-zero error code +if (!is.null(opt$help)) { + cat(getopt(spec, usage=TRUE)) + q(status=1) +} + +if (is.null(opt$gff_file) & is.null(opt$tx2gene)) { + cat("A GFF/GTF file or a tx2gene table is required\n") + q(status=1) +} + +if (opt$format == 'none'){ #custom format + if (is.null(opt$txIdCol) | is.null(opt$abundanceCol) | is.null(opt$countsCol) | is.null(opt$lengthCol)) { + cat("If you select a custom format for the input files you need to specify the column names\n") + q(status=1) + } +} + +if (is.null(opt$countsFiles)) { + cat("'countsFiles' is required\n") + q(status=1) +} + +## parse counts files +library(rjson) +dat <- fromJSON(opt$countsFiles) +samples_df <- lapply(dat, function(samples) # Loop through each "sample" +{ + # Convert each group to a data frame. + # This assumes you have 6 elements each time + data.frame(matrix(unlist(samples), ncol=2, byrow=T)) +}) +samples_df <- do.call(rbind, samples_df) +colnames(samples_df) <- c("path","id") +rownames(samples_df) <- NULL + +# Prepare char vector with files and sample names +files <- file.path(samples_df[,"path"]) +names(files) <- samples_df[,"id"] +#files +#all(file.exists(files)) + + + +library(tximport) + + + + +### if the input is a gff/gtf file first need to create the tx2gene table +if (!is.null(opt$gff_file)) { + suppressPackageStartupMessages({ + library("GenomicFeatures") + }) + txdb <- makeTxDbFromGFF(opt$gff_file) + k <- keys(txdb, keytype = "TXNAME") + tx2gene <- select(txdb, keys=k, columns="GENEID", keytype="TXNAME") + # Remove 'transcript:' from transcript IDs (when gffFile is a GFF3 from Ensembl and the transcript does not have a Name) + tx2gene$TXNAME <- sub('^transcript:', '', tx2gene$TXNAME) + +} else { + tx2gene <- read.table(opt$tx2gene,header=FALSE) + } + + + +## +if (is.null(opt$geneIdCol)) { ## there is a tx2gene table + if (opt$format == 'none'){ #predefined format + cat("here i am too\n") + txi_out <- tximport(files, type="none",txIdCol=opt$txIdCol,abundanceCol=opt$abundanceCol,countsCol=opt$countsCol,lengthCol=opt$lengthCol,tx2gene=tx2gene,countsFromAbundance=opt$countsFromAbundance) + } else { + txi_out <- tximport(files, type=opt$format, tx2gene=tx2gene,countsFromAbundance=opt$countsFromAbundance) + } +} else { # the gene_ID is a column in the counts table + if (opt$format == 'none'){ #predefined format + txi_out <- tximport(files, type="none",geneIdCol=opt$geneIdCol,txIdCol=opt$txIdCol,abundanceCol=opt$abundanceCol,countsCol=opt$countsCol,lengthCol=opt$lengthCol,tx2gene=tx2gene,countsFromAbundance=opt$countsFromAbundance) + } else { + txi_out <- tximport(files, type=opt$format, geneIdCol=opt$geneIdCol,countsFromAbundance=opt$countsFromAbundance) + } + +} +# write count as table +write.table(txi_out$counts, file=opt$out_file, row.names = TRUE, col.names = TRUE, quote = FALSE, sep = "\t") diff -r 000000000000 -r 2f5e9c0fe367 tximport.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tximport.xml Wed Oct 09 15:38:21 2019 -0400 @@ -0,0 +1,228 @@ + + Summarize transcript-level estimates for gene-level analysis + + bioconductor-tximport + bioconductor-genomicfeatures + r-getopt + r-rjson + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +.. class:: infomark + +Current version only works in 'merge' mode: A single table of gene summarizations is generated with one column for each sample file. +Take into account that DEseq2 package in Galaxy requires one table per sample. + + + + doi:10.18129/B9.bioc.tximport + + + +