Skip to content

Commit

Permalink
Fixed RP-SR support problem
Browse files Browse the repository at this point in the history
  • Loading branch information
asylvz committed Nov 29, 2019
1 parent bebc056 commit aa9e5d7
Showing 1 changed file with 2 additions and 1 deletion.
3 changes: 2 additions & 1 deletion vh/vh_setcover.c
Original file line number Diff line number Diff line change
Expand Up @@ -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.
Expand Down

0 comments on commit aa9e5d7

Please sign in to comment.