扩增子分析

1. 质控
2. 拼接 join_paired_ends.py -f 1_1.fq.gz -r 2_2.fq.gz -m fastq-join -o join
3. extract_barcodes.py -f join_paired/fastqjoin.join.fastq -c barcode_single_end --bc1_len 12 -o barcodes
4. split_libraries_fastq.py -i barcodes/reads.fastq -b barcodes/barcodes.fastq -o split -m mapping.txt -q 19

wget https://pypi.python.org/packages/source/s/setuptools/setuptools-19.2.zip
    unzip setuptools-19.2.zip
    cd setuptools-19.2
    python setup.py install

wget https://pypi.python.org/packages/16/e3/06b45eea35359833e7c6fac824b604f1551c2fc7ba0f2bd318d8dd883eb9/cutadapt-1.14.tar.gz 
 cutadapt -g CCGGACTACVSGGGTAT -a WTTCKCGMCTGSTGCRC -e 0.15 -m 300 --discard-untrimmed seqs.fna -o cut_seq.fna

usearch
sed 's/ ./;/g;s/>./&&/g;s/;>/;barcodelabel=/g;s/_[0-9]*;$/;/g' cut_seq.fna > usearch_seq.fa

identify_chimeric_seqs.py -m usearch61 -i split/usearch_seq.fa -o chimeras --suppress_usearch61_ref
filter_fasta.py -f split/usearch_seq.fa -o non_chimeras/non_chimeras_seq.fasta -s chimeras/chimeras.txt -n

pick_otus.py -i non_chimeras/non_chimeras_seq.fasta -o pick_otus_default
pick_otus.py -i non_chimeras/non_chimeras_seq.fasta -m usearch --db /home/llt/database/gold.fa --minsize 2 -o pick_otus_usrarch_ref_gold
pick_rep_set.py -i pick_otus_default/non_chimeras_seq_otus.txt -f non_chimeras/non_chimeras_seq.fasta -o pick_otus_default/rep_set.fna

assign_taxonomy.py -i rep_set.fna -r /home/llt/database/gg_13_8_otus/rep_set/99_otus.fasta -t /home/llt/database/gg_13_8_otus/taxonomy/99_otu_taxonomy.txt -m blast

make_otu_table.py -i non_chimeras_seq_otus.txt -t blast_assigned_taxonomy/rep_set_tax_assignments.txt -o otu_table.biom
biom convert -i otu_table.biom -o table.from_biom_w_consensuslineage.txt --to-tsv --header-key taxonomy --output-metadata-id "ConsensusLineage" --table-type="OTU table"