ZhongPei

 the PS data goes like this

Also, my study is a low biomass study, does it need to be changed to isNotContaminant

setwd("/") library(phyloseq) library(decontam) mapp

M10 is one of my samples, I used contigs.fa assembled by M10 as index for decontamination test. Then, I used bwa-mem2 to mapping the clean_data of 64 samples (including four...

**Thank you very much for the help you can give me, it will be very helpful!**


I think, every feature specific to a Sample(M10). Because in this test, the TPM is counted by the contig.fa assembled by M10_clean_data.fq

Or you mean the present times of each contig_name in the CSV table? 

**Your answer made me think.** 这是我输入的csv文件，可以看到“M10_1001_510”在各个样品中的TPM都不为零。  Also, for each column, my data is barely zero. So "table(colSums(otu>0))" results in mostly 1's, and I don't know why, **is there something...

Thank you very much for the help you can give me, it will be very helpful!