Copy number (repeat count) of assembled contigs

Question

[misaghb]

I have a set of assembled contigs (genomic sequences). These contigs were assembled using a de novo short read assembler (e.g. Velvet, ABySS, etc). I map/align a new short read library (e.g. Illumina paired end or mate pair library) against these contigs using a short read mapper like BOWTIE2. So now I have a SAM/BAM alignment file at hand.

Is there any software (or a pipeline using different tools) to take this alignment SAM/BAM file as input and outputs which contigs are repeated in the genome and estimate those integer repeat counts (copy numbers) using coverage information in the alignment input file?

Notes:

1. No reference genome available.
2. Magnolya is a software that can be used for this purpose although it's not specifically developed for this problem. But it's limited to ace format and some depreciated scripts from tools like ABySS and Newbler.

[Caronkey]

Hi, this is what I want too. I wonder if you get a solution!!

Answer 1

[apelin20]

You can use https://github.com/lindenb/jvarkit/wiki/BamStats04 to caluculate coverage for regions specified in BED file, these regions can be your entire contigs. Go go Pierre.