# HG changeset patch # User petr-novak # Date 1652948515 0 # Node ID b91ca438a1cb9b7c927f2f8a93dba4b6c62e8b27 # Parent 0c3111ab729bed22754bec9b749599cb5760f8bc "planemo upload commit 9633fb98932151f059ce02a0ce202a4374ef8d68" diff -r 0c3111ab729b -r b91ca438a1cb clean_dante_ltr.xml --- a/clean_dante_ltr.xml Mon May 16 07:50:41 2022 +0000 +++ b/clean_dante_ltr.xml Thu May 19 08:21:55 2022 +0000 @@ -1,4 +1,4 @@ - + r-optparse @@ -24,33 +24,30 @@ ]]> - - + + - - - - + + - + + + - + - + \ No newline at end of file diff -r 0c3111ab729b -r b91ca438a1cb clean_ltr.R --- a/clean_ltr.R Mon May 16 07:50:41 2022 +0000 +++ b/clean_ltr.R Thu May 19 08:21:55 2022 +0000 @@ -177,10 +177,14 @@ gff_te <- gff_out[gff_out$type %in% "transposable_element"] gff_5ltr <- gff_out[gff_out$LTR %in% "5LTR"] gff_3ltr <- gff_out[gff_out$LTR %in% "3LTR"] + full_te <- getSeqNamed(s, gff_te) +names(full_te) <- paste0(gff_te$ID,":",names(full_te)) ltr5 <- getSeqNamed(s, gff_5ltr) +names(ltr5) <- paste0(gff_5ltr$Parent,":",names(ltr5)) ltr3 <- getSeqNamed(s, gff_3ltr) -inc <- gff_te$Rank != "DL" +names(ltr3) <- paste0(gff_3ltr$Parent,":",names(ltr3)) +inc <- gff_te$Rank != "DL" writeXStringSet(seq_representative, paste0(opt$output, "_RM_lib_non_redundant.fasta")) writeXStringSet(full_te, paste0(opt$output, "_RM_lib_full_TE.fasta")) diff -r 0c3111ab729b -r b91ca438a1cb dante_ltr_search.xml --- a/dante_ltr_search.xml Mon May 16 07:50:41 2022 +0000 +++ b/dante_ltr_search.xml Thu May 19 08:21:55 2022 +0000 @@ -1,4 +1,4 @@ - + blast r-optparse @@ -13,12 +13,12 @@ mv output.gff3 $te_ltr_gff ]]> - + - +