Skip to content
GitLab
Menu
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Mirrors
biopet.biopet
Commits
39772398
Commit
39772398
authored
May 09, 2016
by
Wai Yi Leung
Browse files
enableRCtrimming is used to enable trimming of Reverse Complement versions of adapters
parent
9e3d9653
Changes
1
Hide whitespace changes
Inline
Side-by-side
flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/Fastqc.scala
View file @
39772398
...
...
@@ -35,6 +35,7 @@ class Fastqc(root: Configurable) extends nl.lumc.sasc.biopet.extensions.Fastqc(r
/** Allow reporting of all found (potentially adapter) sequences in the FastQC */
var
sensitiveAdapterSearch
:
Boolean
=
config
(
"sensitiveAdapterSearch"
,
default
=
false
)
var
enableRCtrimming
:
Boolean
=
config
(
"enableRCtrimming"
,
default
=
false
)
/** Class for storing a single FastQC module result */
protected
case
class
FastQCModule
(
name
:
String
,
status
:
String
,
lines
:
Seq
[
String
])
...
...
@@ -188,9 +189,9 @@ class Fastqc(root: Configurable) extends nl.lumc.sasc.biopet.extensions.Fastqc(r
val
fromKnownList
:
Set
[
AdapterSequence
]
=
(
adapterSet
++
contaminantSet
)
.
filter
(
x
=>
foundAdapterNames
.
exists
(
_
.
startsWith
(
x
.
name
)))
val
fromKnownListRC
:
Set
[
AdapterSequence
]
=
fromKnownList
.
map
{
val
fromKnownListRC
:
Set
[
AdapterSequence
]
=
if
(
enableRCtrimming
)
fromKnownList
.
map
{
x
=>
AdapterSequence
(
x
.
name
+
"_RC"
,
reverseComplement
(
x
.
seq
))
}
}
else
Set
.
empty
// list all sequences found by FastQC
val
fastQCFoundSequences
:
Seq
[
AdapterSequence
]
=
if
(
sensitiveAdapterSearch
)
{
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment