Circos plot using vcf : Mutect2 MultiSample VCF out file
1
0
Entering edit mode
11 months ago
m. bioinfo • 0

Hi

I have multisample vcf out file from Mutect2 and would like to make circos plot showing variants in multiple sample.

Any suggestions or directions please ?

Mutect2 vcf • 1.0k views
ADD COMMENT
0
Entering edit mode

what kind of graphics do you want and what did you try so far ?

ADD REPLY
0
Entering edit mode

Thanks for reply! My expectations are to create such a figure using my multisample vcf file which is a MuTect2 output. First choice of image is published here link Other images are from published articles and mentioned online at circos webpage.

I have tried circos examples and able to create for fusions only on single sample.

Would you have something to suggest or direct me ?

enter image description here

ADD REPLY
0
Entering edit mode
11 months ago

If you want to create one track per sample, the idea is to loop over the samples and to create a file for each samples. Something like (not tested)

bcftools query -l in.vcf.gz | while read S
do
  bcftools view --samples "${S}" in.vcf.gz | bcftools query -f '[%CHROM\t%POS0\t%END\t%SAMPLE\t%GT\n] | grep -v '0/0$'  > "${S}.bed"

done
ADD COMMENT
0
Entering edit mode

thanks Pierre ! and I suspect this .bed output should be used as input for Circos. Do I need any corrections here ? 😊

ADD REPLY
0
Entering edit mode

May I know if you solved this? Thanks a lot!

ADD REPLY

Login before adding your answer.

Traffic: 1462 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6