Dear BioStars Community:
I performed a de-novo assembly of PacBio reads using Canu for a viral genome. I have my resulting unitig from the canu pipeline. I am now interested in characterizing repeats in my resulting assembly. I'm new to assembly and repeat identification. One simple approach was taking the reads used to form the assembly, align them against the assembly with MUMmer, and take a look at the resulting dot plot? Are there any more sophisticated approaches that would yield better results?
1 answer
Yes, you should probably try RepeatModeler, which can detect repeat families (de novo) and classify them. It has worked well for both model/non-model species and is very easy to run (it does have few dependencies to install though: TRF, RECON, RepeatScount, NSEG).
EDIT: you can find my sample run script here!
Log in to answer this question.
-
How to evaluate and check the correctness of the assembly?
written by JXYang •
-
Canu diploid assembly with nanopore reads
written by Irsan
-
Viral de-novo Variant Call
written by Ribas •
-
Correct alignmnent region from Circlator output mapped to orignal assembly
written by ian.will •
-
Canu: failed with exit code 134
written by biomagician
-
Error correcting PacBio reads prior to alignment with bwa mem
written by aindap
-
Aligning de-novo assembly to reference strain genome
written by aindap
-
Do we have any other de novo assembly approaches than Graph based (OLC and De bruijn), Seed Extensi…
written by saranpons3
-
Principled structural variation detection with assembled genomes or pacbio reads
written by hbw
-
Tool: NGS Data Analysis Pipeline for Virus Identification
written by Yang Li
why you do not use Tandem Repeat Finder or RepeatMasker to do this?
Assuming the assembly is correct, it seems to make more sense to align the assembly to itself rather than aligning reads to the assembly.