Skip to content
Snippets Groups Projects
Commit b2226d71 authored by Sander Bollen's avatar Sander Bollen
Browse files

CheckAllelesVcfInBam testing.

parent 31a7e23d
No related branches found
No related tags found
No related merge requests found
...@@ -15,8 +15,11 @@ ...@@ -15,8 +15,11 @@
*/ */
package nl.lumc.sasc.biopet.tools package nl.lumc.sasc.biopet.tools
import java.io.File
import java.nio.file.Paths import java.nio.file.Paths
import htsjdk.samtools.{SamReaderFactory, SamReader}
import htsjdk.variant.vcf.VCFFileReader
import org.scalatest.Matchers import org.scalatest.Matchers
import org.scalatest.mock.MockitoSugar import org.scalatest.mock.MockitoSugar
import org.scalatest.testng.TestNGSuite import org.scalatest.testng.TestNGSuite
...@@ -38,6 +41,7 @@ class CheckAllelesVcfInBamTest extends TestNGSuite with MockitoSugar with Matche ...@@ -38,6 +41,7 @@ class CheckAllelesVcfInBamTest extends TestNGSuite with MockitoSugar with Matche
val vcf = resourcePath("/chrQ.vcf") val vcf = resourcePath("/chrQ.vcf")
val bam = resourcePath("/single01.bam") val bam = resourcePath("/single01.bam")
val vcf2 = new File(resourcePath("/chrQ2.vcf.gz"))
val rand = new Random() val rand = new Random()
@Test def testOutputTypeVcf() = { @Test def testOutputTypeVcf() = {
...@@ -58,4 +62,19 @@ class CheckAllelesVcfInBamTest extends TestNGSuite with MockitoSugar with Matche ...@@ -58,4 +62,19 @@ class CheckAllelesVcfInBamTest extends TestNGSuite with MockitoSugar with Matche
main(arguments) main(arguments)
} }
@Test
def testCheckAllelesNone() = {
val variant = new File(vcf)
val samRecord = SamReaderFactory.makeDefault().open(new File(bam)).iterator().next()
val varRecord = new VCFFileReader(variant, false).iterator().next()
checkAlleles(samRecord, varRecord) shouldBe None
}
@Test
def testCheckAlleles() = {
val samRecord = SamReaderFactory.makeDefault().open(new File(bam)).iterator().next()
val varRecord = new VCFFileReader(vcf2).iterator().next()
checkAlleles(samRecord, varRecord) shouldBe Some("T")
}
} }
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