ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/gclib/tophat_cpp/fusions.cpp
(Generate patch)
# Line 488 | Line 488
488      }  
489   }
490  
491 < void pair_support(const HitsForRead& left_hits, const HitsForRead& right_hits, FusionSet& fusions, FusionSet& fusions_ref)
491 > void pair_support(const vector<pair<BowtieHit, BowtieHit> >& best_hits, FusionSet& fusions, FusionSet& fusions_ref)
492   {
493 <  const vector<BowtieHit>& left = left_hits.hits;
494 <  const vector<BowtieHit>& right = right_hits.hits;
495 <
496 <  if (left.size() > fusion_multipairs || left.size() != right.size())
493 >  if (best_hits.size() > fusion_multipairs)
494      return;
495  
496 <  for (size_t i = 0; i < left.size(); ++i)
496 >  for (size_t i = 0; i < best_hits.size(); ++i)
497      {
498 <      const BowtieHit& lh = left[i];
499 <      const BowtieHit& rh = right[i];
498 >      const BowtieHit& lh = best_hits[i].first;
499 >      const BowtieHit& rh = best_hits[i].second;
500  
501        bool left_fusionSpanned = lh.fusion_opcode() != FUSION_NOTHING;
502        bool right_fusionSpanned = rh.fusion_opcode() != FUSION_NOTHING;

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines