nzimmermann
/
16S_workflow
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Compare commits

base: nzimmermann:21880df2239a140c9f6d87898fd64225249a608b
Branches Tags
nzimmermann:master
...
compare: nzimmermann:8ff28a2a579f92fc54d8833017344375b1aa33b3
Branches Tags
nzimmermann:master

2 Commits
21880df223 ... 8ff28a2a57

Author SHA1 Message Date
nicolas-zimmermann 8ff28a2a57 added otu table creation 2019-10-24 11:59:23 +02:00
nicolas-zimmermann 8165ca41a8 ajout etape clustering 2019-10-24 11:38:01 +02:00
1 changed files with 11 additions and 6 deletions
Show Stats Expand all files Collapse all files

17
16S_workflow.sh
View File

@ -37,13 +37,18 @@ sed -i 's/ //g' fasta/amplicon.fasta
# CREATE OTU # CREATE OTU
# DEDUPLICATION # DEDUPLICATION
mkdir otu mkdir otu
soft/vsearch --derep_fulllength fasta/amplicon.fasta --output otu/fulllength_otu.fasta soft/vsearch --derep_fulllength fasta/amplicon.fasta --output fasta/fulllength.fasta
soft/vsearch --derep_prefix fasta/amplicon.fasta --output otu/prefix_otu.fasta soft/vsearch --derep_prefix fasta/amplicon.fasta --output fasta/prefix.fasta
# REMOVING SINGLETON # REMOVING SINGLETON
soft/vsearch --derep_fulllength fasta/amplicon.fasta --output otu/fulllength_min10_otu.fasta --minuniquesize 10 soft/vsearch --derep_prefix fasta/amplicon.fasta --output fasta/prefix_min10.fasta --minuniquesize 10
soft/vsearch --derep_prefix fasta/amplicon.fasta --output otu/prefix_min10_otu.fasta --minuniquesize 10
# REMOVING CHIMERAS # REMOVING CHIMERAS
soft/vsearch --uchime_denovo otu/prefix_min10_otu.fasta --nonchimeras otu/prefix_min10_nochim_otu.fasta soft/vsearch --uchime_denovo fasta/prefix.fasta --nonchimeras fasta/prefix_nochim.fasta
soft/vsearch --uchime_denovo otu/fulllength_min10_otu.fasta --nonchimeras otu/fulllength_min10_nochim_otu.fasta
# CLUSTERING
otu=">OTU_"
vsearch --cluster_fast fasta/prefix_nochim.fasta --id 0.97 --centroids otu/centroids.fasta --relabel $otu
# ABONDANCE
vsearch --usearch_global fasta/amplicon.fasta --db otu/clusters.fasta --id 0.97 --otutabout otu/abondance_table.txt