@@ -91,22 +91,22 @@ These databases were created from Unite v8.0, first downloading two databases fr
9191.. code-block :: none
9292
9393 #Create full length ITS USEARCH Database, convert taxonomy, and create USEARCH database
94- amptk database -i UNITE_public_all_02.02.2019.fasta -f ITS1-F -r ITS4 \
95- --primer_required none -o ITS --create_db usearch --install --source UNITE:8.0
96-
97- #create SINTAX database
94+ amptk database -i UNITE_public_all_02.02.2019.fasta -f ITS1-F -r ITS4 \
95+ --primer_required none -o ITS --create_db usearch --install --source UNITE:8.0
96+
97+ #create SINTAX database
9898 amptk database -i sh_general_release_dynamic_all_02.02.2019_dev.fasta \
9999 -o ITS_SINTAX --create_db utax -f ITS1-F -r ITS4 --derep_fulllength \
100- --install --source UNITE:8.0 --primer_required none
100+ --install --source UNITE:8.0 --primer_required none
101101
102102 #Create UTAX Databases
103103 amptk database -i sh_general_release_dynamic_all_02.02.2019_dev.fasta \
104104 -o ITS_UTAX --create_db utax -f ITS1-F -r ITS4 \
105105 --derep_fulllength --install --source UNITE:8.0 --primer_required none
106106
107- amptk database -i sh_general_release_dynamic_all_02.02.2019_dev.fasta \
108- -o ITS1_UTAX -f ITS1-F -r ITS2 --primer_required rev --derep_fulllength \
109- --create_db utax --install --subsample 65000 --source UNITE:8.0
107+ amptk database -i sh_general_release_dynamic_all_02.02.2019_dev.fasta \
108+ -o ITS1_UTAX -f ITS1-F -r ITS2 --primer_required rev --derep_fulllength \
109+ --create_db utax --install --subsample 65000 --source UNITE:8.0
110110
111111 amptk database -i sh_general_release_dynamic_all_02.02.2019_dev.fasta \
112112 -o ITS2_UTAX --create_db utax -f fITS7 -r ITS4 --derep_fulllength \
@@ -130,63 +130,63 @@ Since it can literally take days to download the arthropod dataset, if you'd lik
130130 cat arthropods.bold-reformated.fa chordates.bold-reformated.fa > arth-chord.bold-reformated.fasta
131131
132132 #generate global alignment database
133- amptk database -i arth-chord.bold.reformated.fasta -f LCO1490 -r mlCOIintR --primer_required none \
134- --derep_fulllength --format off --primer_mismatch 4 -o COI --min_len 200 --create_db usearch \
135- --install --source BOLD:20190219
133+ amptk database -i arth-chord.bold.reformated.fasta -f LCO1490 -r mlCOIintR --primer_required none \
134+ --derep_fulllength --format off --primer_mismatch 4 -o COI --min_len 200 --create_db usearch \
135+ --install --source BOLD:20190219
136136
137137 The second set of output files from `bold2utax.py ` are named with `.BIN-consensus.fa ` which are the result of 99% clustering for each BIN. We will combine those for the two datasets and then use those data to generate the SINTAX and UTAX databases.
138138
139139.. code-block :: none
140140
141- #combine datasets
142- cat arthropods.BIN-consensus.fa chordates.BIN-consensus.fa > arth-chord.bold.BIN-consensus.fasta
141+ #combine datasets
142+ cat arthropods.BIN-consensus.fa chordates.BIN-consensus.fa > arth-chord.bold.BIN-consensus.fasta
143143
144- #generate SINTAX database
145- amptk database -i arth-chord.bold.BIN-consensus.fasta -f LCO1490 -r mlCOIintR --primer_required none \
146- --derep_fulllength --format off --primer_mismatch 4 -o COI_SINTAX --min_len 200 --create_db sintax \
147- --install --source BOLD:20190219
148-
149- #generate UTAX database, need to subsample for memory issues with 32 bit usearch and we require rev primer match here
150- amptk database -i arth-chord.bold.BIN-consensus.fasta -f LCO1490 -r mlCOIintR --primer_required rev \
151- --derep_fulllength --format off --subsample 30000 --primer_mismatch 4 -o COI_UTAX --min_len 200 \
152- --create_db utax --install --source BOLD:20190219
144+ #generate SINTAX database
145+ amptk database -i arth-chord.bold.BIN-consensus.fasta -f LCO1490 -r mlCOIintR --primer_required none \
146+ --derep_fulllength --format off --primer_mismatch 4 -o COI_SINTAX --min_len 200 --create_db sintax \
147+ --install --source BOLD:20190219
148+
149+ #generate UTAX database, need to subsample for memory issues with 32 bit usearch and we require rev primer match here
150+ amptk database -i arth-chord.bold.BIN-consensus.fasta -f LCO1490 -r mlCOIintR --primer_required rev \
151+ --derep_fulllength --format off --subsample 00000 --primer_mismatch 4 -o COI_UTAX --min_len 200 \
152+ --create_db utax --install --source BOLD:20190219
153153
154154 **LSU database **
155155
156156The fungal 28S database (LSU) was downloaded from `RDP <http://rdp.cme.msu.edu/download/current_Fungi_unaligned.fa.gz >`_. The sequences were then converted into AMPtk databases as follows:
157157
158158.. code-block :: none
159159
160- amptk database -i RDP_v8.0_fungi.fa -o LSU --format rdp2utax --primer_required none \
161- --skip_trimming --create_db usearch --derep_fulllength --install --source RDP:8
160+ amptk database -i RDP_v8.0_fungi.fa -o LSU --format rdp2utax --primer_required none \
161+ --skip_trimming --create_db usearch --derep_fulllength --install --source RDP:8
162162
163- amptk database -i RDP_v8.0_fungi.fa -o LSU_SINTAX --format rdp2utax --primer_required none \
164- --skip_trimming --create_db sintax --derep_fulllength --install --source RDP:8
163+ amptk database -i RDP_v8.0_fungi.fa -o LSU_SINTAX --format rdp2utax --primer_required none \
164+ --skip_trimming --create_db sintax --derep_fulllength --install --source RDP:8
165165
166- amptk database -i RDP_v8.0_fungi.fa -o LSU_UTAX --format rdp2utax --primer_required none \
167- --skip_trimming --create_db utax --derep_fulllength --install --source RDP:8 --subsample 4500
168-
169-
166+ amptk database -i RDP_v8.0_fungi.fa -o LSU_UTAX --format rdp2utax --primer_required none \
167+ --skip_trimming --create_db utax --derep_fulllength --install --source RDP:8 --subsample 45000
168+
169+
170170 To generate a training set for UTAX, the sequences were first dereplicated, and clustered at 97% to get representative sequences for training. This training set was then converted to a UTAX database:
171171
172172.. code-block :: none
173173
174- amptk database -i fungi.trimmed.fa -o LSU_UTAX --format off \
175- --skip_trimming --create_db utax --keep_all --install
174+ amptk database -i fungi.trimmed.fa -o LSU_UTAX --format off \
175+ --skip_trimming --create_db utax --keep_all --install
176176
177177 **16S database **
178178This is downloaded from `R. Edgar's website <http://drive5.com/utax/data/rdp_v16.tar.gz >`_ and then formatted for AMPtk. Note there is room for substantial improvement here, I just don't typically work on 16S - so please let me know if you want some suggestions on what to do here. Here I reformatted the "domain" taxonomy level to "kingdom" for simplicity (even though I know it is taxonomically incorrect).
179179
180180.. code-block :: none
181181
182- amptk database -i rdp_16s_v16_sp.kingdom.fa -o 16S --format off --create_db usearch \
183- --skip_trimming --install --primer_required none --derep_fulllength
184-
185- amptk database -i rdp_16s_v16_sp.kingdom.fa -o 16S --format off --create_db sintax \
186- -f 515FB -r 806RB --install --primer_required for --derep_fulllength
187-
188-
189-
182+ amptk database -i rdp_16s_v16_sp.kingdom.fa -o 16S --format off --create_db usearch \
183+ --skip_trimming --install --primer_required none --derep_fulllength
184+
185+ amptk database -i rdp_16s_v16_sp.kingdom.fa -o 16S_SINTAX --format off --create_db sintax \
186+ -f 515FB -r 806RB --install --primer_required for --derep_fulllength
187+
188+ amptk database -i rdp_16s_v16_sp.kingdom.fa -o 16S_UTAX --format off --create_db sintax \
189+ -f 515FB -r 806RB --install --primer_required for --derep_fulllength
190190
191191
192192 Checking Installed Databases
0 commit comments