瀏覽代碼

ajout commandes

nicolas-zimmermann 5 年之前
父節點
當前提交
21880df223
共有 1 個文件被更改,包括 8 次插入4 次删除
  1. 8 4
      16S_workflow.sh

+ 8 - 4
16S_workflow.sh 查看文件

@@ -37,9 +37,13 @@ sed -i 's/ //g' fasta/amplicon.fasta
37 37
 # CREATE OTU
38 38
 # DEDUPLICATION
39 39
 mkdir otu
40
-soft/vsearch --derep_fulllength fasta/amplicon.fasta --output otu/fulllength_otu
41
-soft/vsearch --derep_prefix fasta/amplicon.fasta --output otu/prefix_otu
40
+soft/vsearch --derep_fulllength fasta/amplicon.fasta --output otu/fulllength_otu.fasta
41
+soft/vsearch --derep_prefix fasta/amplicon.fasta --output otu/prefix_otu.fasta
42 42
 
43 43
 # REMOVING SINGLETON
44
-soft/vsearch --derep_fulllength fasta/amplicon.fasta --output otu/fulllength_min10_otu --minuniquesize 10
45
-soft/vsearch --derep_prefix fasta/amplicon.fasta --output otu/prefix_min10_otu --minuniquesize 10
44
+soft/vsearch --derep_fulllength fasta/amplicon.fasta --output otu/fulllength_min10_otu.fasta --minuniquesize 10
45
+soft/vsearch --derep_prefix fasta/amplicon.fasta --output otu/prefix_min10_otu.fasta --minuniquesize 10
46
+
47
+# REMOVING CHIMERAS
48
+soft/vsearch --uchime_denovo otu/prefix_min10_otu.fasta --nonchimeras otu/prefix_min10_nochim_otu.fasta
49
+soft/vsearch --uchime_denovo otu/fulllength_min10_otu.fasta --nonchimeras otu/fulllength_min10_nochim_otu.fasta