# HG changeset patch # User jackcurragh # Date 1663941345 0 # Node ID 190fa66855a803e82b50b005f20ef0940ab87b66 # Parent c81aa97acf65b36b1af0b868a1891694e10d22ef Uploaded diff -r c81aa97acf65 -r 190fa66855a8 get_chrom_sizes/get_chrom_sizes.xml --- a/get_chrom_sizes/get_chrom_sizes.xml Tue Aug 30 10:54:24 2022 +0000 +++ b/get_chrom_sizes/get_chrom_sizes.xml Fri Sep 23 13:55:45 2022 +0000 @@ -1,4 +1,4 @@ - + Generates a TAB Delimited chrom.sizes File from an Inputted FASTA #if $refGenomeSource.genomeSource == "builtin": - cat ${refGenomeSource.input1_builtin.fields.path} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf ${input2} substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 + if ${input2} == "chr": + cat ${refGenomeSource.input1_builtin.fields.path} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf "chr" substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 + else: + cat ${refGenomeSource.input1_builtin.fields.path} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 + #else: + if ${input2} == "chr": + cat ${refGenomeSource.input1_file} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf "chr" substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 + else: + cat ${refGenomeSource.input1_file} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 + #end if + + + + cat ${refGenomeSource.input1_builtin.fields.path} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf ${input2} substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 #else: cat ${refGenomeSource.input1_file} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf ${input2} substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1 #end if @@ -35,7 +48,6 @@ -