Assembly Aligned Paired-End Reads

Hi all,

I have a set of mapped paired-end reads and I would like to assemble the ones that overlap respecting the pairing information.

This means assemblying only pairs when the 2 first mates overlap and the 2 second mates overlap too. The reads are already mapped on a genome, there is nothing more to do with the sequences, only with the positions.

The goal is to get the extended positions with the count information.

Pairs example:

chr5:1456-1498,+     chr5:1654-1702,+
chr5:958-1012,+      chr5:1318-1388,+
chr5:1423-1478,+     chr5:1612-1667,+

I would like to get:

2     chr5:1423-1498,+     chr5:1612-1702,+
1     chr5:958-1012,+       chr5:1318-1388,+

I can't find any software working on the positions, all I can find is FLASH, PEAR, etc. which are working on the fastq files.

Cheers