Entering edit mode
6.4 years ago
mks002
▴
220
Dear all, I have performed SNV calling on plasmodium samples. However I am facing some difficulty to interpreting certain fields. I wish to annotate this .vcf file further.
Here is an example from one of the .vcf file
Pf3D7_02_v3 903199 BND00000329 A A]Pf3D7_01_v3:84434] . LowQual IMPRECISE;SVTYPE=BND;SVMETHOD=EMBL.DELLYv0.7.7;CHR2=Pf
3D7_01_v3;END=84434;PE=3;MAPQ=1;CT=3to3;CIPOS=-115,115;CIEND=-115,115 GT:GL:GQ:FT:RCL:RC:RCR:CN:DR:DV:RR:RV 0/0:0,-9.17065,-38.1212:92:P
ASS:12520:0:21281:0:48:0:0:0
Pf3D7_03_v3 90256 INV00000039 A <INV> . LowQual IMPRECISE;SVTYPE=INV;SVMETHOD=EMBL.DELLYv0.7.7;CHR2=Pf3D7_03_v3;END=90
384;PE=2;MAPQ=40;CT=5to5;CIPOS=-55,55;CIEND=-55,55 GT:GL:GQ:FT:RCL:RC:RCR:CN:DR:DV:RR:RV 0/0:0,-0.927125,-97.7996:10:LowQual:2513:642
2:3688:2:25:2:0:0
Pf3D7_03_v3 282969 DEL00000154 A <DEL> . LowQual IMPRECISE;SVTYPE=DEL;SVMETHOD=EMBL.DELLYv0.7.7;CHR2=Pf3D7_03_v3;END=28
3430;PE=2;MAPQ=42;CT=3to5;CIPOS=-153,153;CIEND=-153,153 GT:GL:GQ:FT:RCL:RC:RCR:CN:DR:DV:RR:RV 0/0:0,-6.04816,-184.599:60:PASS:13141:28020:
14221:2:46:2:0:0
Pf3D7_03_v3 598312 INV00000040 T <INV> . LowQual IMPRECISE;SVTYPE=INV;SVMETHOD=EMBL.DELLYv0.7.7;CHR2=Pf3D7_03_v3;END=598470;PE=2;MAPQ=36;CT=5to5;CIPOS=-236,236;CIEND=-236,236 GT:GL:GQ:FT:RCL:RC:RCR:CN:DR:DV:RR:RV 0/0:0,-3.72035,-65.2937:37:PASS:5282:4380:3205:1:25:1:0:0
Here are my specific doubts
- CT=5to5,3to5,5to5,3to3 what do these correspond to?
- CIEND=-55,55 why are these values in negative value and positive value of the same number.
- I performed svprops to convert .vcf to .bed and then to perform annotation did a bed intersect with the gtf file. The difficulty arises in case of Interchromosomal translocations (represented by tag BND). Because the start field corresponds to the original chromosome and end corresponds to the end co-ordinates of the chromosome of translocation( correct if this point is not right). So how to annotate these translocation events?
4.How well does latest SNPeff annotates SVs?
Thanks for help!