I want to know if there is some kind of pipeline i can use in order to annotate a given list of structural variants detected by NGS.
specifically, i just ran Breakdancer on a paired end sequencing data and recieved the output specifying the different types of predicted variations (insertions / deletions / inversions / translocations) example :
1 1413910 9+0- 1 1452528 0+10- DEL 38636 99
1 1588571 0+17- 1 1653965 17+0- ITX 64960 99
the columns are : chr1, position1, supportingreads(forward/reverse), chr2, position2, supportingreads(forward/reverse), variationtype, variationsize, quality
i basically want to annotate these variants in regards to gene name, exon/intron, transcription binding sites, ncRNA (miRNA) and any other known feature regarding the genetic region affected by them.
is there a known tool for this? if not, i would appreciate any suggestions on how to perform such a task (for example some sort of outline for a workflow i can build)
thanks.
•
written
by
Doctoroots
Bedtools may help you
There happens to be a complete pipeline that results in IGV visualization as well as flat files for tertiary analysis.
Look at this publication:
http://bioinformatics.oxfordjournals.org/content/early/2011/11/14/bioinformatics.btr612.full.pdf
Can be downloaded Here:
http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm
--> This is a front to back full pipeline. However you can see that there are option capabilities to run from any step. Either all, just alignment, just variants, just annotation.
** It's important to note that this pipeline is based on using .BAMS and .vcf as standard file types for modular employment. You'll need to refer to pac-bio or some other source for converters.
•
written
by
Raymond301 •
*** Note that there is also a Cloud deployment in case you don't have the computational resources to run this pipeline.
hi raymond, although the tool looks useful for other purposes i couldnt see how it handles structural variants such as translocations and inversions
You are correct. Upon closer inspection translocations & inversions are not readily annotated through this software. I apologize.
I have basically the same question. I have to create a SV pipeline as a project on my master. We also have to include some sort of annotation. For now we annotate the output via a GFF file. The GFF contains all the positions of genes, exons, introns etc. We compare these positions with our SV list.
Furthermore we compare SV calls from different software tools to see which SV are supported by multiple SV callers.
I am really wondering what other people's thoughts are about annotating SV's.
Log in to answer this question.
How did you end up annotating your SV data. I find myself in the same kind of situation.
SV annotation (with OMIM, DGV, 1000g, haploinsufficiency, TAD, ... and also with your own in-house information) can be easily automated !
You can look at this post describing the annotSV tool: Annotation for SV and CNV