diff --git a/tools/halfdeep/halfdeep.xml b/tools/halfdeep/halfdeep.xml index 0654347134d..db29477d589 100644 --- a/tools/halfdeep/halfdeep.xml +++ b/tools/halfdeep/halfdeep.xml @@ -1,11 +1,11 @@ - identifies genomic regions with half-depth coverage based on sqeuencing read mappings. These regions may reveal insights into heterogametic sex chromosomes, haplotype-specific variation, or potential assembly errors such as heterotypic duplications. + identifies genomic regions with half-depth coverage based on sequencing read mappings. macros.xml > input.fofn && + #set $read_base = re.sub('[^\w\-\s]', '_', str($read.element_identifier)) + ln -s '$read' '$reads_dir/${read_base}.$read.ext' && + echo '$reads_dir/${read_base}.$read.ext' >> input.fofn && ## ## mapped reads ## #for $mapped_read in $mapped_reads - ln -s '$mapped_read' "$mapped_reads_dir/${seq_base}.bam" && - ln -s "${seq_base}.bam" "$mapped_reads_dir/${seq_base}.sort.bam" && - ln -s '$mapped_read.metadata.bam_index' "$mapped_reads_dir/${seq_base}.sort.bam.bai" && + ln -s '$mapped_read' "$mapped_reads_dir/${read_base}.bam" && + ln -s "${read_base}.bam" "$mapped_reads_dir/${read_base}.sort.bam" && + ln -s '$mapped_read.metadata.bam_index' "$mapped_reads_dir/${read_base}.sort.bam.bai" && #end for #end for ## @@ -43,7 +43,7 @@ ]]> - + @@ -64,35 +64,28 @@