1
1
Entering edit mode
2.0 years ago

Hi all,

I extracted some reads from a region (~200kb) of interest (including the unmapped ones) and all the pairs of unmapped reads from my WGS BAM: I'd like to know if it is possible to automatically (using vg?) detect any structural variant using the spades assembler (http://cab.spbu.ru/files/release3.12.0/manual.html ) for this region of interest .

The output directory looks like:

assembly_graph_after_simplification.gfa
assembly_graph.fastg
assembly_graph_with_scaffolds.gfa
before_rr.fasta
contigs.fasta
contigs.paths
corrected
dataset.info
input_dataset.yaml
K21
K33
K55
K77
misc
params.txt
pipeline_state
scaffolds.fasta
scaffolds.paths
tmp


is there any way to automatically detect any large inversion/insertion/deletion/etc.. from the assembly graph ?

spades vg sv • 837 views
0
Entering edit mode

Do you want to identify SVs by comparing your assembly with the reference genome?

0
Entering edit mode

yes. I wonder if the unmapped reads can give some new insights compared to the usual algorithms (manta, etc...) .

1
Entering edit mode

I don't know what you already tried, but Assemblytics was developed for that purpose, should work but is not very recent. Might be okay for the length of contigs you have, but more problematic with contigs from long read sequencing. Other alternatives are SVIM-asm (which I haven't tried, but SVIM is good) and dipcall (which I know does only identify SVs span by an alignment). Please let me know what works for you :)

0
Entering edit mode

very interesting, I'll try it tomorrow

1
Entering edit mode
22 months ago

Was just browsing and saw this question, maybe it is too late but oh well, I wrote a tool called MUM&Co to compare de novo assemblies and I've tried it on many illumina assemblies and it seem to work well. Also there is Syri that does similar work. As they both explore the full range of SV types unlike Assemblytics.

Hope that helps