@@ -149,11 +149,11 @@ class CorrectUmisTest extends UnitSpec {
149149 readBamRecs(rejects).map(_.name) should contain theSameElementsAs Seq (" q9" , " q10" )
150150
151151 val metricsByUmi = Metric .read[UmiCorrectionMetrics ](metrics).map(m => m.umi -> m).toMap
152- metricsByUmi(" AAAAAA" ) shouldBe UmiCorrectionMetrics (umi= " AAAAAA" , total_matches= 5 , perfect_matches= 2 , one_mismatch_matches= 1 , two_mismatch_matches= 1 , other_matches = 1 , fraction_of_matches = 5 / 10.0 , representation = (5 / 8.0 ) / (1 / 4d ))
153- metricsByUmi(" CCCCCC" ) shouldBe UmiCorrectionMetrics (umi= " CCCCCC" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
154- metricsByUmi(" GGGGGG" ) shouldBe UmiCorrectionMetrics (umi= " GGGGGG" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
155- metricsByUmi(" TTTTTT" ) shouldBe UmiCorrectionMetrics (umi= " TTTTTT" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
156- metricsByUmi(" NNNNNN" ) shouldBe UmiCorrectionMetrics (umi= " NNNNNN" , total_matches= 2 , perfect_matches= 0 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
152+ metricsByUmi(" AAAAAA" ) shouldBe UmiCorrectionMetrics (umi= " AAAAAA" , total_matches= 5 , perfect_matches= 2 , one_mismatch_matches= 1 , two_mismatch_matches= 1 , fraction_of_matches = 5 / 10.0 , representation = (5 / 8.0 ) / (1 / 4d ))
153+ metricsByUmi(" CCCCCC" ) shouldBe UmiCorrectionMetrics (umi= " CCCCCC" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
154+ metricsByUmi(" GGGGGG" ) shouldBe UmiCorrectionMetrics (umi= " GGGGGG" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
155+ metricsByUmi(" TTTTTT" ) shouldBe UmiCorrectionMetrics (umi= " TTTTTT" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
156+ metricsByUmi(" NNNNNN" ) shouldBe UmiCorrectionMetrics (umi= " NNNNNN" , total_matches= 2 , perfect_matches= 0 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
157157 }
158158
159159 it should " run end to end and do the right thing for a handful of fragments with duplex UMIs" in {
@@ -177,11 +177,11 @@ class CorrectUmisTest extends UnitSpec {
177177 readBamRecs(rejects).map(_.name) should contain theSameElementsAs Seq (" q3" , " q5" )
178178
179179 val metricsByUmi = Metric .read[UmiCorrectionMetrics ](metrics).map(m => m.umi -> m).toMap
180- metricsByUmi(" AAAAAA" ) shouldBe UmiCorrectionMetrics (umi= " AAAAAA" , total_matches= 3 , perfect_matches= 2 , one_mismatch_matches= 1 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 3 / 10.0 , representation = (3 / 8.0 ) / (1 / 4d ))
181- metricsByUmi(" CCCCCC" ) shouldBe UmiCorrectionMetrics (umi= " CCCCCC" , total_matches= 2 , perfect_matches= 1 , one_mismatch_matches= 1 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
182- metricsByUmi(" GGGGGG" ) shouldBe UmiCorrectionMetrics (umi= " GGGGGG" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
183- metricsByUmi(" TTTTTT" ) shouldBe UmiCorrectionMetrics (umi= " TTTTTT" , total_matches= 2 , perfect_matches= 2 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
184- metricsByUmi(" NNNNNN" ) shouldBe UmiCorrectionMetrics (umi= " NNNNNN" , total_matches= 2 , perfect_matches= 0 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , other_matches = 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
180+ metricsByUmi(" AAAAAA" ) shouldBe UmiCorrectionMetrics (umi= " AAAAAA" , total_matches= 3 , perfect_matches= 2 , one_mismatch_matches= 1 , two_mismatch_matches= 0 , fraction_of_matches = 3 / 10.0 , representation = (3 / 8.0 ) / (1 / 4d ))
181+ metricsByUmi(" CCCCCC" ) shouldBe UmiCorrectionMetrics (umi= " CCCCCC" , total_matches= 2 , perfect_matches= 1 , one_mismatch_matches= 1 , two_mismatch_matches= 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
182+ metricsByUmi(" GGGGGG" ) shouldBe UmiCorrectionMetrics (umi= " GGGGGG" , total_matches= 1 , perfect_matches= 1 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 1 / 10.0 , representation = (1 / 8.0 ) / (1 / 4d ))
183+ metricsByUmi(" TTTTTT" ) shouldBe UmiCorrectionMetrics (umi= " TTTTTT" , total_matches= 2 , perfect_matches= 2 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
184+ metricsByUmi(" NNNNNN" ) shouldBe UmiCorrectionMetrics (umi= " NNNNNN" , total_matches= 2 , perfect_matches= 0 , one_mismatch_matches= 0 , two_mismatch_matches= 0 , fraction_of_matches = 2 / 10.0 , representation = (2 / 8.0 ) / (1 / 4d ))
185185 }
186186
187187 it should " do the same thing whether the UMIs are on the command line or in a file" in {
0 commit comments