Mercurial > repos > iuc > qiime_split_libraries
diff generate_test_data.sh @ 5:6f2a0b6a8da2 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/qiime/ commit 3b54163c4f7daff76fcc589c4a9057bb03904380
author | iuc |
---|---|
date | Sat, 05 Aug 2017 07:15:16 -0400 |
parents | a4caa409b9cd |
children |
line wrap: on
line diff
--- a/generate_test_data.sh Mon Jul 10 16:46:48 2017 -0400 +++ b/generate_test_data.sh Sat Aug 05 07:15:16 2017 -0400 @@ -228,6 +228,54 @@ --suppress_group_significance rm -rf core_diversity_analyses_2 +# extract_barcodes +extract_barcodes.py \ + --fastq1 'test-data/extract_barcodes/inseqs.fastq' \ + --input_type 'barcode_single_end' \ + -o extract_barcodes_1 \ + --bc1_len '6' \ + --rev_comp_bc1 +rm -rf extract_barcodes_1 + +extract_barcodes.py \ + --fastq1 'test-data/extract_barcodes/inseqs_R1.fastq' \ + --input_type 'barcode_paired_end' \ + --fastq2 'test-data/extract_barcodes/inseqs_R2.fastq' \ + -o extract_barcodes_2 \ + --bc1_len '6' \ + --bc2_len '6' +rm -rf extract_barcodes_2 + +extract_barcodes.py \ + --fastq1 'test-data/extract_barcodes/inseqs_R1.fastq' \ + --input_type 'barcode_paired_end' \ + --fastq2 'test-data/extract_barcodes/inseqs_R2.fastq' \ + -o extract_barcodes_3 \ + --bc1_len '6' \ + --bc2_len '6' \ + --mapping_fp 'test-data/extract_barcodes/mapping_data.txt' \ + --attempt_read_reorientation \ + --disable_header_match +rm -rf extract_barcodes_3 + +extract_barcodes.py \ + --fastq1 'test-data/extract_barcodes/inseqs_R1.fastq' \ + --input_type 'barcode_paired_stitched' \ + -o extract_barcodes_4 \ + --bc1_len '6' \ + --bc2_len '8' \ + --rev_comp_bc1 \ + --rev_comp_bc2 +rm -rf extract_barcodes_4 + +extract_barcodes.py \ + --fastq1 'test-data/extract_barcodes/inseqs_R1.fastq' \ + --input_type 'barcode_in_label' \ + --char_delineator '#' \ + -o extract_barcodes_5 \ + --bc1_len '6' +rm -rf extract_barcodes_5 + # filter_alignment filter_alignment.py \ --input_fasta_file 'test-data/filter_alignment/alignment.fasta' \ @@ -943,23 +991,79 @@ # split_libraries_fastq split_libraries_fastq.py \ - --sequence_read_fps 'test-data/split_libraries_fastq/forward_reads.fastq' \ - -o split_libraries \ - --mapping_fps 'test-data/map.tsv' \ - --barcode_read_fps 'test-data/split_libraries_fastq/barcodes.fastq' \ + --sequence_read_fps 'test-data/split_libraries_fastq/lane1_read1.fastq.gz' \ + -o split_libraries_1 \ + --mapping_fps 'test-data/split_libraries_fastq/map.txt' \ + --barcode_read_fps 'test-data/split_libraries_fastq/lane1_barcode.fastq.gz' \ + --max_bad_run_length 3 \ + --min_per_read_length_fraction 0.75 \ + --sequence_max_n 0 \ + --start_seq_id 0 \ + --rev_comp_mapping_barcodes \ + --phred_quality_threshold 19 \ + --barcode_type 'golay_12' \ + --max_barcode_errors 1.5 +rm -rf split_libraries_1 + +split_libraries_fastq.py \ + --sequence_read_fps 'test-data/split_libraries_fastq/lane1_read1.fastq.gz' \ + -o split_libraries_2 \ + --mapping_fps 'test-data/split_libraries_fastq/map.txt' \ + --barcode_read_fps 'test-data/split_libraries_fastq/lane1_barcode.fastq.gz' \ --store_qual_scores \ --store_demultiplexed_fastq \ --max_bad_run_length 3 \ --min_per_read_length_fraction 0.75 \ --sequence_max_n 0 \ --start_seq_id 0 \ + --rev_comp_mapping_barcodes \ + --phred_quality_threshold 19 \ + --barcode_type 'golay_12' \ + --max_barcode_errors 1.5 +rm -rf split_libraries_2 + +split_libraries_fastq.py \ + --sequence_read_fps 'test-data/split_libraries_fastq/lane1_read1.fastq.gz,test-data/split_libraries_fastq/lane2_read1.fastq.gz' \ + -o split_libraries_3 \ + --mapping_fps 'test-data/split_libraries_fastq/map.txt,test-data/split_libraries_fastq/map.txt' \ + --barcode_read_fps 'test-data/split_libraries_fastq/lane1_barcode.fastq.gz,test-data/split_libraries_fastq/lane2_barcode.fastq.gz' \ + --max_bad_run_length 3 \ + --min_per_read_length_fraction 0.75 \ + --sequence_max_n 0 \ + --start_seq_id 0 \ + --rev_comp_mapping_barcodes \ + --phred_quality_threshold 19 \ --barcode_type 'golay_12' \ --max_barcode_errors 1.5 -cp split_libraries/histograms.txt 'test-data/split_libraries_fastq/histograms.tabular' -cp split_libraries/seqs.fna 'test-data/split_libraries_fastq/sequences.fasta' -cp split_libraries/seqs.qual 'test-data/split_libraries_fastq/sequence_qualities.qual' -cp split_libraries/seqs.fastq 'test-data/split_libraries_fastq/demultiplexed_sequences.fastq' -rm -rf split_libraries +rm -rf split_libraries_3 + +split_libraries_fastq.py \ + --sequence_read_fps 'test-data/split_libraries_fastq/lane1_read1.fastq.gz' \ + -o split_libraries_4 \ + --sample_ids 'my.sample.1' \ + --max_bad_run_length 3 \ + --min_per_read_length_fraction 0.75 \ + --sequence_max_n 0 \ + --start_seq_id 0 \ + --rev_comp_mapping_barcodes \ + --phred_quality_threshold 19 \ + --barcode_type 'not-barcoded' \ + --max_barcode_errors 1.5 +rm -rf split_libraries_4 + +split_libraries_fastq.py \ + --sequence_read_fps 'test-data/split_libraries_fastq/lane1_read1.fastq.gz,test-data/split_libraries_fastq/lane2_read1.fastq.gz' \ + -o split_libraries_5 \ + --sample_ids 'my.sample.1,my.sample.2' \ + --max_bad_run_length 3 \ + --min_per_read_length_fraction 0.75 \ + --sequence_max_n 0 \ + --start_seq_id 0 \ + --rev_comp_mapping_barcodes \ + --phred_quality_threshold 19 \ + --barcode_type 'not-barcoded' \ + --max_barcode_errors 1.5 +rm -rf split_libraries_5 # summarize_taxa cp 'test-data/core_diversity_analyses/otu_table.biom' 'test-data/summarize_taxa/otu_table.biom'