Browse Source

added annotation part

nicolas-zimmermann 4 years ago
parent
commit
721e404673
1 changed files with 5 additions and 2 deletions
  1. 5 2
      16S_workflow.sh

+ 5 - 2
16S_workflow.sh View File

@@ -48,7 +48,10 @@ soft/vsearch --uchime_denovo fasta/prefix.fasta --nonchimeras fasta/prefix_nochi
48 48
 
49 49
 # CLUSTERING
50 50
 otu=">OTU_"
51
-vsearch --cluster_fast fasta/prefix_nochim.fasta --id 0.97 --centroids otu/centroids.fasta --relabel $otu
51
+vsearch --cluster_fast fasta/prefix_nochim.fasta --id 0.97 --centroids otu/otu.fasta --relabel $otu
52 52
 
53 53
 # ABONDANCE
54
-vsearch --usearch_global fasta/amplicon.fasta --db otu/clusters.fasta --id 0.97 --otutabout otu/abondance_table.txt
54
+vsearch --usearch_global fasta/amplicon.fasta --db otu/otu.fasta --id 0.97 --otutabout otu/abondance_table.txt
55
+
56
+# ANNOTATION
57
+soft/vsearch --usearch_global otu/otu.fasta --db databases/mock_16S_18S.fasta --id 0.97 --top_hits_only --userfields query+target --userout otu/annotation.txt