...式如下:列之间用tab分隔,第一行为SNP在基因组上的位置信息;从第二行开始,第一列为样品名称,后面为对应SNP。
Strain 325 968 2351 4768 4980 6971 9323 1319
ARK0000038 G C C G G A G T
ARK0000056 G C C G G A G T
ARK0000062 G C C G G A G T
ARK0000066 G C C G ...
#去掉序列ID中注释信息,避免错误:no ASV ids overlap between input FASTA and sequence abundance table
get_fa_by_id.pl $workdir/5.pick_otu_qiime/pick_de_novo_otus/otu_table_clean.txt $workdir/5.pick_otu_qiime/pick_de_novo_otus/rep_set/qiime_rep_set.fasta rep_set.farm -rf picrust2_out...