How to get gene abundance visualization and venn diagram from PROKKA output file?
0
0
Entering edit mode
4.7 years ago
nisrinalulu ▴ 10

Hi Everyone,

I hope you all doing well.

I work with metagenomics data from soil samples and have annotated it using PROKKA. I need visualize PROKKA output to get gene abundance visualization and venn diagram. Do you have a suggestion for me what software or what R packages that I can use to get gene abundance visualization and venn diagram from PROKKA output?

Thank for your help and time.

Regards, Lulu

metagenomics R • 2.3k views
ADD COMMENT
0
Entering edit mode

what exactly do you want to see visualised in a venn diagram?

if you're not afraid of some file reformatting, this webtool might be of use?: http://bioinformatics.psb.ugent.be/webtools/Venn/

ADD REPLY
0
Entering edit mode

Dear lieven.sterck,

Thank you for your reply and your web tool recommendation. The web tool that you mentioned looks easy to use but I have to convert the file format.

I want to see gene in venn diagram among 4 my samples.

Do you have another software recommendation?

Thank you

ADD REPLY
0
Entering edit mode

well, in any case if you want to represent it in a venn diagram, you will need some way of linking the genes to each other (based on naming for instance, but that will likely not work in this case).

ADD REPLY
0
Entering edit mode

I see. Thank you for your reply. I think I can use R packages to produce gene abundance visualization and venn diagram from my data set, but I still looking for the best packages that I can use to produced it.

ADD REPLY
0
Entering edit mode

Dear lieven.streck,

Do you have a tutorial or paper to reformatting PROKKA output file to suitable format file for http://bioinformatics.psb.ugent.be/webtools/Venn/?

Thank you

ADD REPLY
0
Entering edit mode

no sorry, don't have such a document.

however, input is pretty simple, it's just a list of IDs (== elements you want to compare, can be anything actually) in plain text format.

key issue remains that you will need to find some way to link the genes (?0 to each other. perhaps do a AllVsAll blast search and take the best hits? I assume that simply taking the predicted gene IDs will not do it as you can't be sure they point to the same gene.

ADD REPLY

Login before adding your answer.

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