Copy number (repeat count) of de nove assembly contigs
I have a set of assembled contigs (genomic sequences). These contigs were assembled using a de novo short read assembler . I map/align the original reads data to these contigs using a short read mapper like bwa. 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?
• 929 views
•
link
written
by
Caronkey •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to remove circular contigs from assembly.fasta
written by kristina.mahan
-
Extracting alignment coordinates from SAM file
written by dasakan •
-
Extract the mapped contig sequences from SAM/BAM file
written by BDK_compbio
-
Does SPades contigs.fasta contain unassembled reads?
written by ciclistadan •
-
Inferring most probable integer copy number (repeat count) of contigs of a draft genome
written by misaghb •
-
Copy number (repeat count) of assembled contigs
written by misaghb •
-
Is there any open source assembly softwares that takes in de novo assembly generated contigs and re…
written by ajagannath.patro •
-
To remove or to keep duplicates in alignment of NGS paired reads to a draft genome (set of assemble…
written by misaghb •
-
De Novo Genome Assembler Preferring Shorter Error-Free Contigs
written by misaghb •
-
Tools For Visualizing Rna-Seq Splicing Graphs From De Novo Assembly Without A Reference Genome?
written by Tianyang Li •
Might be me, but I'm a bit confused why you map new data to the assembly to get an idea on repeated contigs in your assembly? Why don't you align all contigs to themselves (eg. using minimap?) and extract that kind of info form that result?