This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VEP annotation of Manta BND

I am using VEP to annotate SVs from Manta vcf files. For the translocations (MantaBND) only the first reported chr region gets a gene annoation. Is there a way to get gene annotations for both chromosomes involved in the translocation?

annnotation vep

duplicate the variant to create another one with the 2nd position.

What format are your translocations inputted in?

The structure of the VCF is:

chr1 84799235 MantaBND:6124:1:2:0:0:0:0 A A]chr14:85914438] . PASS SVTYPE=BND;MATEID=MantaBND:6124:1:2:0:0:0:1;IMPRECISE;CIPOS=-92,92;SOMATIC;SOMATICSCORE=34;BND_DEPTH=27;MATE_BND_DEPTH=31 PR

@Emily_Ensembl the breakpoint is defined in the ALT column. A]chr14:85914438]

Yes, I can see that.

Sorry, I see now that Manta actually reports breakend translocations on double rows so 'chr14 85914438' comes further down a second row, which means that it gets a gene annotatation by VEP. Sorry, my bad...

0 answers

No answers yet.

Log in to answer this question.