From aa9e5d7ec5c666a950bffa3067e07ecb27217444 Mon Sep 17 00:00:00 2001 From: asylvz Date: Fri, 29 Nov 2019 14:07:42 +0300 Subject: [PATCH] Fixed RP-SR support problem --- vh/vh_setcover.c | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/vh/vh_setcover.c b/vh/vh_setcover.c index d2713bc..ad3d2d0 100644 --- a/vh/vh_setcover.c +++ b/vh/vh_setcover.c @@ -965,12 +965,13 @@ float calWeight( bam_info **in_bams, parameters *params, int clusterId, int *cou { if( read_names[ptrReadMapping->readId].readCovered == 0 && listClusterEl[clusterId].indIdCount[ptrReadMapping->indId] > -1) { - supOfIndSeen[ptrReadMapping->indId]++; /* Count the number of split reads supporting this cluster */ if( strcmp( multiLibs[read_names[ptrReadMapping->readId].libId].libName, "SplitRead") == 0) listClusterEl[clusterId].sr_support[ptrReadMapping->indId]++; + supOfIndSeen[ptrReadMapping->indId]++; + //////////////////////////////////////THE HUERISTIC////////////////////////// //Instead of taking the total support as number of reads in each cluster, we use summation of //normalized number of reads. i.e each read contributes as 1/(total number of mappings) it has.