Aquarium data – Round 2

Guillaume re-demultiplexed the Aquarium data for us – we found one barcode that was used twice (so had to throw out those samples, 42.SedimentCoralpond.2 and 0.WaterFreshwater.2). Newest demultiplexed data allows 2 mismatches per merged dual index long barcode (so 1 per mate barcode). Concatenated data on Merlot and downloaded to my iMac:

cat /share/eisen-z2/gjospin/Slims_Eisenlab_backup/AQUEXP1/demul_9_reverse/qiime_ready/AQUEXP1.MPF.*.faa > AQUEXP1.MPF.all_reads_merged_20Jun.faa
cat /share/eisen-z2/gjospin/Slims_Eisenlab_backup/AQUEXP1/demul_9_reverse/qiime_ready/AQUEXP1.MPR.*.faa > AQUEXP1.MPR.all_reads_merged_20Jun.faa
cat /share/eisen-z2/gjospin/Slims_Eisenlab_backup/AQUEXP1/demul_9_reverse/qiime_ready/AQUEXP1.M.*.faa > AQUEXP1.M.all_reads_merged_20Jun.faa

And then kicked off the QIIME analysis (uclust open reference OTU picking, with parameters file enabling reverse strand match):

pick_open_reference_otus.py -i /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/AQUEXP1.MPF.all_reads_merged_20Jun.faa -r /macqiime/greengenes/gg_12_10_otus/rep_set/99_otus.fasta -o /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun --parallel -O 2 -p /Users/hollybik/Desktop/Data/aquarium_project/mapping_files/aquarium_parametersfile.txt

Continuation 6/28/13 – Next filtered out the Koi pond samples, and any samples that had less than 1000 observations (sequences):

filter_samples_from_otu_table.py -m /Users/hollybik/Desktop/Data/aquarium_project/mapping_files/AquariumSampleMap8.txt -n 1000 -o /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/filtered_otu_tables/otu_table_mc2_w_tax_no_pynast_failures_filtered1000seqs_noKoiPond.biom -i /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/otu_table_mc2_w_tax_no_pynast_failures.biom --sample_id_fp /Users/hollybik/Desktop/Data/aquarium_project/mapping_files/samples_to_keep_noKoi.txt

Removed genus Brachybacterium from samples:

filter_taxa_from_otu_table.py -i /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/filtered_otu_tables/otu_table_mc2_w_tax_no_pynast_failures_filtered1000seqs_noKoiPond.biom -o /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/filtered_otu_tables/otu_table_mc2_w_tax_no_pynast_failures_filtered1000seqs_noKoiPond_noBrachybacterium.biom -n g__Brachybacterium

Filtered out “dead OTUs” (no sequences left after filtering):

filter_otus_from_otu_table.py -i /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/filtered_otu_tables/otu_table_mc2_w_tax_no_pynast_failures_filtered1000seqs_noKoiPond_noBrachybacterium.biom -o /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/filtered_otu_tables/otu_table_mc2_w_tax_no_pynast_failures_filtered1000seqs_noKoiPond_noBrachybacterium_nodeadOTUs.biom -n 1

Then running the core_diversity_analyses.py script:

core_diversity_analyses.py -i /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/filtered_otu_tables/otu_table_mc2_w_tax_no_pynast_failures_filtered1000seqs_noKoiPond_noBrachybacterium_nodeadOTUs.biom -o /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/core_diversity_analyses -m /Users/hollybik/Desktop/Data/aquarium_project/mapping_files/AquariumSampleMap8.txt -e 90 -t /Users/hollybik/Desktop/Data/aquarium_project/demul_trim_prep_2mismatcheslongBC/uclust_openref_20Jun/rep_set.tre -c SampleType,Location,SampleTypeLocation,Days
Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: