Question: Qiime file for Krona pie chart?
1
gravatar for manekineko
2.8 years ago by
manekineko130
Bulgaria
manekineko130 wrote:

Hi,

I have analysis with Qiime of 16s reads, and I want to make a pie chart with Krona but do not know which file from Qiime to use for Krona?

qiime krona • 2.5k views
ADD COMMENTlink modified 19 months ago by ferstar40 • written 2.8 years ago by manekineko130

You'll be more successful asking at the qiime forum: https://groups.google.com/forum/#!forum/qiime-forum

ADD REPLYlink written 2.8 years ago by Daniel3.7k
5
gravatar for ferstar
19 months ago by
ferstar40
China/Xi'an/NWU
ferstar40 wrote:

Hi, hope this answer is not too late. :)

convert biom to tsv format

biom convert -i single_sample.biom -o single_sample.tsv --to-tsv --table-type "OTU table" --header-key taxonomy

remove the first two rows with comment

# take a look at single_sample.tsv
$ head single_sample.tsv 
# Constructed from biom file
#OTU ID H2O taxonomy
346085  140.0   Bacteria; Proteobacteria; Alphaproteobacteria; Caulobacterales; Caulobacteraceae; Brevundimonas; Brevundimonas_bullata
10298   2.0 Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacteriales; Enterobacteriaceae; Photorhabdus; Photorhabdus_temperata
122823  3.0 Bacteria; Proteobacteria; Betaproteobacteria; Burkholderiales; Oxalobacteraceae; Massilia; EF516371_s
130468  2.0 Bacteria; Proteobacteria; Alphaproteobacteria; Sphingomonadales; Sphingomonadaceae; Sphingopyxis; Sphingopyxis_witflariensis
139977  38.0    Bacteria; Proteobacteria; Alphaproteobacteria; Sphingomonadales; Erythrobacteraceae; Erythrobacter; Erythrobacter_flavus
121751  2.0 Bacteria; Firmicutes; Clostridia; Clostridiales; Lachnospiraceae; Catonella; JX096343_s
96934   10.0    Bacteria; Proteobacteria; Gammaproteobacteria; Pseudomonadales; Pseudomonadaceae; Pseudomonas; Pseudomonas_guguanensis
95181   4.0 Bacteria; Proteobacteria; Gammaproteobacteria; Pseudomonadales; Moraxellaceae; Acinetobacter; Acinetobacter_radioresistens

# remove the comments
$ egrep -v "^#" single_sample.tsv > sample_no_comment.tsv

remove the first column

$ cut sample_no_comment.tsv -f 2- > sample.tsv
$ head -n2 sample.tsv 
140.0   Bacteria; Proteobacteria; Alphaproteobacteria; Caulobacterales; Caulobacteraceae; Brevundimonas; Brevundimonas_bullata
2.0 Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacteriales; Enterobacteriaceae; Photorhabdus; Photorhabdus_temperata

replace the "; " with "\t"

$ sed -i 's/;\s*/\t/g' sample.tsv
$ head -n2 sample.tsv 
140.0   Bacteria    Proteobacteria  Alphaproteobacteria Caulobacterales Caulobacteraceae    Brevundimonas   Brevundimonas_bullata
2.0 Bacteria    Proteobacteria  Gammaproteobacteria Enterobacteriales      Enterobacteriaceae   Photorhabdus    Photorhabdus_temperata

now we can make a pie chart with ktImportText(One of the Krona tools)

$ ktImportText -n enjoy_your_life -o sample.krona.html sample.tsv

open sample.krona.html with your web browser, you'll see this

Krona_snapshot_enjoy_your_life

all in one command(optional)

$ awk -F '\t|;' 'BEGIN{OFS="\t"} FNR > 2 {$1=""; print $0}' single_sample.tsv | cut -f 2- > sample.tsv
ADD COMMENTlink modified 19 months ago • written 19 months ago by ferstar40

krona is awesome! I used it with kraken output glad to know qiime output can be used as well

ADD REPLYlink written 19 months ago by badribio240

Yeah, its really awesome and I'm glad my solution is useful to you:)

ADD REPLYlink modified 19 months ago • written 19 months ago by ferstar40

Hello badribio, Could you please explain how to create kraken output,what is the command to be used??

ADD REPLYlink written 13 months ago by arunsaikirans0

What have you done so far, and which step are you stuck?

ADD REPLYlink written 13 months ago by badribio240

I have finished qiime analysis,but I want to know what output file generated by qiime to provide as an input in kraken.

ADD REPLYlink written 13 months ago by arunsaikirans0
1
gravatar for biokoala
2.2 years ago by
biokoala3.6k
biokoala3.6k wrote:

otu_table_L6.txt file generated by the script summarize_taxa.py will be the best file to be used for generating Krona pie chart for obvious reason (having information for all the taxonomic ranks right from phylum to genus).

ADD COMMENTlink written 2.2 years ago by biokoala3.6k

Hi Vijay, I have installed Krona and Kraken what is the input file for kraken and krona??

ADD REPLYlink written 13 months ago by arunsaikirans0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2049 users visited in the last hour