Mercurial > repos > ecology > ab1_fastq_converter

annotate ab1_fastq.R @ 0:fef3f3584168 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
author ecology
date Wed, 12 Jan 2022 15:12:28 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
1 #Rscript
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
2
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
3 ############################################
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
4 ## Convert ab1 files into fastq ##
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
5 ############################################
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
6
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
7 #####Packages
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
8
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
9 library(CrispRVariants, quietly = TRUE)
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
10 library(sangerseqR, quietly = TRUE)
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
11
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
12 #####Load arguments
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
13
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
14 args <- commandArgs(trailingOnly = TRUE)
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
15
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
16 if (length(args) == 0) {
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
17 stop("This tool needs at least one argument")
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
18 } else {
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
19 file <- args[1]
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
20 filename <- args[2]
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
21 tr <- as.logical(args[3])
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
22 co <- as.numeric(args[4])
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
23 min_seq <- as.integer(args[5])
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
24 os <- as.numeric(args[6])
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
25 }
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
26
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
27 ##### Conversion
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
28
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
29 if (grepl("^.+\\.[aA][bB][1i]$", filename)) {
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
30 nfile <- sub("^(.+)\\.[aA][bB][1i]$", "\\1", filename)
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
31 } else {
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
32 nfile <- filename
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
33 }
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
34
fef3f3584168 "planemo upload for repository https://github.com/ColineRoyaux/Galaxy_tool_projects/tree/main/ab1_fastq commit dbecaa89a5afa0cc73ae00a716c98ae46fa97b58"
ecology
parents:
diff changeset
35 CrispRVariants::abifToFastq(nfile, file, "output.fastq", trim = tr, cutoff = co, min_seq_len = min_seq, offset = os)