I'm looking for a tools that detect structural variations (SVs, genome rearrangement) between sequenced genome and its reference.
I have got 45 bp paired-ends reads (Illumina HiSeq), ~40x coverage from 13Mb haploid yeast strain.
I know that there are some tools that can estimate SVs from reads alignments (eg SVDetect, BreakDancer). Another approach is to align de novo assembly onto reference genome and get SVs from there (eg Mummer). I would like you know, which is better approach, and how would you do it ideally? Maybe someone has got ready pipeline for that?
5 answers
Here are a couple of tools for SV detection you may want to look into that I've been reading about lately (just today!), but haven't tried:
Any feedback on how these systems perform would be helpful for the community.
I myself have not tried using any of these tools, but this type of analysis is in my TODO list. I found this a couple of days ago and it looks very promising:
http://bioinf.spbau.ru/en/sibelia
From their website:
Sibelia: A comparative genomics tool: It assists biologists in analysing the genomic variations that correlate with pathogens, or the genomic changes that help microorganisms adapt in different environments. Sibelia will also be helpful for the evolutionary and genome rearrangement studies for multiple strains of microorganisms.
Sibelia is useful in finding: (1) shared regions, (2) regions that present in one group of genomes but not in others, (3) rearrangements that transform one genome to other genomes.
In version 2, Sibelia works with multiple strains of bacteria and partitions their genomes into synteny blocks — blocks of highly conserved regions among all compared genomes. It represents genomes in circos pictures [for publication] or interactive forms [for experts’ analysis].
Sibelia is under active development. If you see that Sibelia has a potential to support your research, please do not hesitate to contact us at vyahhi@bioinf.spbau.ru with a list of features you would like Sibelia to have.
I've tried breakway and breakdancer as well. I was really surprised, that these programmes give completely inconsistent results! Could anyone comment on that?
Even Sibelia has a web interface:
can we use C-Sibelia OR sibelia for more than two genome, if so what would be the command line
Log in to answer this question.