How to draw circular plot to show genomewide snps and indels distribution?
1
0
Entering edit mode
7.6 years ago
mirza ▴ 180

Hi,

I have 2 related fungal isolates and am working on variant calling. I have called variants using samtools-bcftools and have annotated them using snpEff. What tools can I use to draw a circular plot to represent snps and indels distribution?

I want to generate a plot like this Genome analysis of rice-blast fungus Magnaporthe oryzae field isolates, Figure 3

Here is the link to the paper https://doi.org/10.1016/j.gdata.2015.06.018

circular plot snp distribution • 3.0k views
ADD COMMENT
1
Entering edit mode

The paper you linked includes a citation for circos which was used to generate that figure.

ADD REPLY
0
Entering edit mode

Another option is BRIG (though the resolution of the image on large genomes can make snps etc hard to see)

ADD REPLY
1
Entering edit mode

If you want to do it in R, I found circlize quite useful and easy to understand. circlize Cran

ADD REPLY
0
Entering edit mode

thank you all for your reply.

ADD REPLY
2
Entering edit mode
7.6 years ago
Titus ▴ 910

Hi , Looks like you are looking for some circos plot , there is this discussion on biostars which can interest you : Circos-How To Generate Circos Image Of Snp And Indel Across Chromosomes

ADD COMMENT

Login before adding your answer.

Traffic: 1825 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