# HG changeset patch # User davidvanzessen # Date 1448445774 18000 # Node ID 5a886d753c0a3171d59483585e3b38647e0e22fb # Parent 7973dd7a39ab1eee443b80c73032b04cb9ebb30e Uploaded diff -r 7973dd7a39ab -r 5a886d753c0a demultiplex.xml --- a/demultiplex.xml Tue Nov 04 09:16:17 2014 -0500 +++ b/demultiplex.xml Wed Nov 25 05:02:54 2015 -0500 @@ -4,7 +4,7 @@ fastx_toolkit - wrapper.sh $input $out_file $out_file.files_path $where $mismatches $partial $input.name + wrapper.sh "$input" "$out_file" "$out_file.files_path" "$where" "$mismatches" "$partial" "$input.name" #for $i, $b in enumerate($barcodes) "$b.id" "$b.mid" @@ -340,7 +340,7 @@ -- Splitting sff files into FASTQ, FASTA and (optional) trimmed FASTA files with a FASTQC report on the FASTQ file, this tool uses: +- Splitting sff or fastq files into FASTQ, FASTA and (optional) trimmed FASTA files with a FASTQC report on the FASTQ file, this tool uses: - sff2fastq (https://github.com/indraniel/sff2fastq) to extract a fastq file. - fastx_barcode_splitter.pl (http://hannonlab.cshl.edu/fastx_toolkit/commandline.html) to demultiplex. - fastqc (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) to provide analysis of the fastq files. diff -r 7973dd7a39ab -r 5a886d753c0a wrapper.sh --- a/wrapper.sh Tue Nov 04 09:16:17 2014 -0500 +++ b/wrapper.sh Wed Nov 25 05:02:54 2015 -0500 @@ -1,11 +1,11 @@ #!/bin/bash -input=$1 -output=$2 -outDir=$3 -mkdir $outDir -EOL=$4 -mismatches=$5 -partial=$6 +input="$1" +output="$2" +outDir="$3" +mkdir "$outDir" +EOL="$4" +mismatches="$5" +partial="$6" name=$(basename "$7") ext="${name##*.}" name="${name%.*}" @@ -23,7 +23,7 @@ start_int=$((i+2)) end_int=$((i+3)) id="${!i}" - echo "$id" + echo "$id, ${start_int}, ${end_int}" trim_start[$id]=${!start_int} trim_end[$id]=${!end_int} echo -e "$id\t${!j}" >> $outDir/barcodes.txt