INNER CODE UNIT · C

gap

aquaskyline/SOAPdenovo2 · fusion/attachPEinfo.c:338

  gap = insert_size - overlaplen + pre_pos + pos - contig_array[e1].length - contig_array[e2].length;

  //fprintf(stderr,"[%s]\tgap\t%d\t%d\t%f\t%f\n",__FUNCTION__,gap,insert_size,close_threshold,insert_size*close_threshold);
  if(gap < -(insert_size * close_threshold))
    {
      ignorePE2++;
      return 0;
    }

  if(gap > insert_size)
    {
      ignorePE3++;
      return 0;
    }

  add1AccuConnect(e1, e2, gap, 1);
  add1AccuConnect(bal_e2, bal_e1, gap, 1);

View source record →