# HG changeset patch
# User iuc
# Date 1727994949 0
# Node ID d3b0390f325f8e1fea410ae62b8624a8d3c0557c
# Parent 05f8e9d1e8b8a74c3e172a8689d3fc6f1ef07b14
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/brew3r_r commit d287d5688e57f4154d5875789e0cd4d5c87f11ce
diff -r 05f8e9d1e8b8 -r d3b0390f325f brew3r.r_script.R
--- a/brew3r.r_script.R Thu Jun 27 07:31:23 2024 +0000
+++ b/brew3r.r_script.R Thu Oct 03 22:35:49 2024 +0000
@@ -97,6 +97,18 @@
}
}
+if (is.null(input_gr_to_extend$exon_id)) {
+ is.exon <- which(input_gr_to_extend$type == "exon")
+ input_gr_to_extend$exon_id <- NA
+ input_gr_to_extend$exon_id[is.exon] <- paste0(
+ "EXON",
+ sprintf(
+ "%010d",
+ 1:length(is.exon)
+ )
+ )
+}
+
# Run BREW3R.r main function
if (length(input_gr_template) > 0) {
new_gr_exons <- extend_granges(
diff -r 05f8e9d1e8b8 -r d3b0390f325f brew3r_r.xml
--- a/brew3r_r.xml Thu Jun 27 07:31:23 2024 +0000
+++ b/brew3r_r.xml Thu Oct 03 22:35:49 2024 +0000
@@ -2,7 +2,7 @@
Extend GTF
1.0.2
- 0
+ 1
brew3r.r
@@ -67,6 +67,16 @@
+
+
+
+
+
diff -r 05f8e9d1e8b8 -r d3b0390f325f test-data/generate_test.R
--- a/test-data/generate_test.R Thu Jun 27 07:31:23 2024 +0000
+++ b/test-data/generate_test.R Thu Oct 03 22:35:49 2024 +0000
@@ -258,5 +258,7 @@
new.gr <- extend_granges(input_gr, big_gr)
library("rtracklayer")
export.gff(input_gr, "input.gtf")
+input_gr$exon_id <- NULL
+export.gff(input_gr, "input_noexonid.gtf")
export.gff(big_gr, "second_input.gtf")
export.gff(sort(new.gr, ignore.strand = TRUE), "output.gtf")
diff -r 05f8e9d1e8b8 -r d3b0390f325f test-data/input_noexonid.gtf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/input_noexonid.gtf Thu Oct 03 22:35:49 2024 +0000
@@ -0,0 +1,31 @@
+##gff-version 2
+##source-version rtracklayer 1.64.0
+##date 2024-10-02
+chr1 rtracklayer exon 5 10 . + . gene_id "gene11"; transcript_id "transcript11"; gene_name "gene11"
+chr1 rtracklayer exon 20 30 . + . gene_id "gene12"; transcript_id "transcript12"; gene_name "gene12"
+chr1 rtracklayer exon 105 110 . + . gene_id "gene21"; transcript_id "transcript21"; gene_name "gene21"
+chr1 rtracklayer exon 120 125 . + . gene_id "gene22"; transcript_id "transcript22"; gene_name "gene22"
+chr1 rtracklayer exon 205 210 . + . gene_id "gene31"; transcript_id "transcript31"; gene_name "gene31"
+chr1 rtracklayer exon 220 222 . + . gene_id "gene32"; transcript_id "transcript32"; gene_name "gene32"
+chr1 rtracklayer exon 305 310 . + . gene_id "gene41"; transcript_id "transcript41"; gene_name "gene41"
+chr1 rtracklayer exon 305 322 . + . gene_id "gene42"; transcript_id "transcript42"; gene_name "gene42"
+chr1 rtracklayer exon 405 410 . + . gene_id "gene51"; transcript_id "transcript51"; gene_name "gene51"
+chr1 rtracklayer exon 405 425 . + . gene_id "gene52"; transcript_id "transcript52"; gene_name "gene52"
+chr1 rtracklayer exon 505 510 . + . gene_id "gene61"; transcript_id "transcript61"; gene_name "gene61"
+chr1 rtracklayer exon 520 522 . + . gene_id "gene62"; transcript_id "transcript62"; gene_name "gene62"
+chr1 rtracklayer exon 601 610 . + . gene_id "gene71"; transcript_id "transcript71"; gene_name "gene71"
+chr1 rtracklayer exon 601 610 . + . gene_id "gene71"; transcript_id "transcript72"; gene_name "gene71"
+chr1 rtracklayer exon 630 640 . + . gene_id "gene71"; transcript_id "transcript72"; gene_name "gene71"
+chr1 rtracklayer exon 701 710 . + . gene_id "gene81"; transcript_id "transcript81"; gene_name "gene81"
+chr1 rtracklayer exon 701 710 . + . gene_id "gene82"; transcript_id "transcript82"; gene_name "gene82"
+chr1 rtracklayer exon 730 740 . + . gene_id "gene82"; transcript_id "transcript82"; gene_name "gene82"
+chr1 rtracklayer exon 801 808 . + . gene_id "gene1"; transcript_id "transcript91"; gene_name "gene1"
+chr1 rtracklayer exon 801 810 . + . gene_id "gene1"; transcript_id "transcript92"; gene_name "gene1"
+chr1 rtracklayer exon 830 840 . + . gene_id "gene1"; transcript_id "transcript92"; gene_name "gene1"
+chr1 rtracklayer exon 901 908 . + . gene_id "gene101"; transcript_id "transcript101"; gene_name "gene101"
+chr1 rtracklayer exon 901 910 . + . gene_id "gene102"; transcript_id "transcript102"; gene_name "gene102"
+chr1 rtracklayer exon 930 940 . + . gene_id "gene102"; transcript_id "transcript102"; gene_name "gene102"
+chr1 rtracklayer exon 1005 1010 . + . gene_id "gene111"; transcript_id "transcript111"; gene_name "Gm001"
+chr1 rtracklayer exon 1055 1070 . + . gene_id "gene112"; transcript_id "transcript112"; gene_name "gene112"
+chr1 rtracklayer exon 1100 1105 . + . gene_id "gene121"; transcript_id "transcript121"; gene_name "gene121"
+chr1 rtracklayer exon 1110 1120 . - . gene_id "gene122"; transcript_id "transcript122"; gene_name "gene122"