@@ -535,6 +535,60 @@ test_that("linkMultipleVariants works", {
535535 expect_equal(res $ filtSeparate $ V3 [, " nbrRetained" ],
536536 sum(! grepl(" C.*mut|del|V3\\ ." , truth $ truth $ status )))
537537
538+ # --------------------------------------------------------------------------
539+ # Get warning when summing UMIs
540+ expect_warning({
541+ res <- linkMultipleVariants(
542+ combinedDigestParams = list (
543+ fastqForward = system.file(" extdata" , " multipleVariableRegions_R1.fastq.gz" ,
544+ package = " mutscan" ),
545+ fastqReverse = system.file(" extdata" , " multipleVariableRegions_R2.fastq.gz" ,
546+ package = " mutscan" ),
547+ elementsForward = " UVCVCV" , elementLengthsForward = c(6 , 24 , 10 , 30 , 10 , - 1 ),
548+ elementsReverse = " CVCV" , elementLengthsReverse = c(6 , 40 , 10 , - 1 ),
549+ mergeForwardReverse = TRUE , minOverlap = 20 ,
550+ maxOverlap = 30 , maxFracMismatchOverlap = 0 ,
551+ revComplForward = FALSE , revComplReverse = TRUE ,
552+ avePhredMinForward = 20 , avePhredMinReverse = 20 ,
553+ verbose = verbose ),
554+ barcode = list (
555+ fastqForward = system.file(" extdata" , " multipleVariableRegions_R1.fastq.gz" ,
556+ package = " mutscan" ),
557+ elementsForward = " CVCSCS" , elementLengthsForward = c(6 , 24 , 10 , 30 , 10 , - 1 ),
558+ avePhredMinForward = 20 ,
559+ collapseMaxDist = 1 , collapseMinScore = 1 , collapseMinRatio = 1 ,
560+ verbose = verbose ),
561+ V2 = list (
562+ fastqForward = system.file(" extdata" , " multipleVariableRegions_R1.fastq.gz" ,
563+ package = " mutscan" ),
564+ fastqReverse = system.file(" extdata" , " multipleVariableRegions_R2.fastq.gz" ,
565+ package = " mutscan" ),
566+ elementsForward = " CSCVCS" , elementLengthsForward = c(6 , 24 , 10 , 30 , 10 , - 1 ),
567+ elementsReverse = " CSCV" , elementLengthsReverse = c(6 , 40 , 10 , - 1 ),
568+ mergeForwardReverse = TRUE , minOverlap = 10 ,
569+ maxOverlap = 20 , maxFracMismatchOverlap = 0 ,
570+ revComplForward = FALSE , revComplReverse = TRUE ,
571+ avePhredMinForward = 20 , avePhredMinReverse = 20 ,
572+ wildTypeForward = truth $ WTV2 ,
573+ nbrMutatedCodonsMaxForward = - 1 , nbrMutatedBasesMaxForward = 1 ,
574+ forbiddenMutatedCodonsForward = " " , collapseToWTForward = FALSE ,
575+ verbose = verbose , useTreeWTmatch = TRUE ),
576+ V3 = list (
577+ fastqForward = system.file(" extdata" , " multipleVariableRegions_R1.fastq.gz" ,
578+ package = " mutscan" ),
579+ fastqReverse = system.file(" extdata" , " multipleVariableRegions_R2.fastq.gz" ,
580+ package = " mutscan" ),
581+ elementsForward = " CSCSCV" , elementLengthsForward = c(6 , 24 , 10 , 30 , 10 , - 1 ),
582+ elementsReverse = " CVCS" , elementLengthsReverse = c(6 , 40 , 10 , - 1 ),
583+ mergeForwardReverse = TRUE , minOverlap = 5 ,
584+ maxOverlap = 15 , maxFracMismatchOverlap = 0 ,
585+ revComplForward = FALSE , revComplReverse = TRUE ,
586+ avePhredMinForward = 20 , avePhredMinReverse = 20 ,
587+ wildTypeForward = truth $ WTV3 ,
588+ nbrMutatedCodonsMaxForward = - 1 , nbrMutatedBasesMaxForward = 1 ,
589+ forbiddenMutatedCodonsForward = " " , collapseToWTForward = FALSE ,
590+ verbose = verbose , useTreeWTmatch = TRUE )
591+ )}, " Aggregating UMI counts" )
538592})
539593
540594
0 commit comments