Skip to content
Snippets Groups Projects
Commit 9eb01b6b authored by Wai Yi Leung's avatar Wai Yi Leung
Browse files

Fix FastQC test since we now include the reverse complement of a found...

Fix FastQC test since we now include the reverse complement of a found adapter, the test had to be modified to check the change.
parent fe444f51
No related branches found
No related tags found
No related merge requests found
......@@ -185,10 +185,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] = (adapterSet ++ contaminantSet)
.filter(x => foundAdapterNames.exists(_.startsWith(x.name))).map {
x => AdapterSequence(x.name + "_RC", reverseComplement(x.seq))
}
val fromKnownListRC: Set[AdapterSequence] = fromKnownList.map {
x => AdapterSequence(x.name + "_RC", reverseComplement(x.seq))
}
// list all sequences found by FastQC
val fastQCFoundSequences: Seq[AdapterSequence] = if (sensitiveAdapterSearch) {
......
......@@ -71,10 +71,16 @@ class FastqcV0101Test extends TestNGSuite with Matchers {
val fqc = new Fastqc(null)
fqc.output = outputv0101
fqc.contaminants = Option(resourceFile("fqc_contaminants_v0101.txt"))
// found adapters also contain the adapters in reverse complement (done within flexiprep/fastqc only)
val adapters = fqc.foundAdapters
adapters.size shouldBe 1
adapters.head.name should ===("TruSeq Adapter, Index 1")
// from fqc_contaminants_v0101.txt
adapters.head.seq should ===("GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG")
// we find 1 adapter which comes with the Reverse Complement counterpart
adapters.size shouldBe 2
adapters.head.name shouldEqual "TruSeq Adapter, Index 1_RC"
adapters.head.seq shouldEqual "CAAGCAGAAGACGGCATACGAGATCGTGATGTGACTGGAGTTCAGACGTGTGCTCTTCCGATC"
adapters.last.name shouldEqual "TruSeq Adapter, Index 1"
adapters.last.seq shouldEqual "GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG"
}
}
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment