script to calculate posistional coverage of next gen contigs from mapping data of reads that constructed contigs