comparison htseq-count.xml @ 21:19e80c732870 draft

planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/htseq_count commit 8a40cf16ce0b48cdfda88a505869e77e8826cb23

20:6b61ba6dd9e2

            #set $reference_fasta_filename = str( $samout_conditional.reference_source.ref_file.fields.path )
        #end if
    #end if
    #if $force_sort:
        #if $samfile.extension == 'bam':
            samtools sort -n $samfile "name_sorted_alignment" &&
        #else
            samtools view -Su -t ${reference_fasta_filename}.fai $samfile | samtools sort -n - "name_sorted_alignment" &&
        #end if
    #end if
    htseq-count
    --mode=$mode
    --stranded=$stranded
    --minaqual=$minaqual
    --type=$featuretype
    --idattr=$idattr
    #if $samout_conditional.samout:
        --samout=$__new_file_path__/${samoutfile.id}_tmp
    #end if
    #if $force_sort:
        --order=name

    #else
        --order=pos
        --format=$samfile.extension
        $samfile
    #end if
    $gfffile
    | awk '{if ($1 ~ "no_feature|ambiguous|too_low_aQual|not_aligned|alignment_not_unique") print $0 | "cat 1>&2"; else print $0}' > $counts 2>$othercounts
    #if $samout_conditional.samout:
        && samtools view -Su -t ${reference_fasta_filename}.fai $__new_file_path__/${samoutfile.id}_tmp | samtools sort -o - sorted > $samoutfile
    #end if
    ]]>
    </command>

    <inputs>
        <param format="sam,bam" name="samfile" type="data" label="Aligned SAM/BAM File"/>
        <param format="gff" name="gfffile" type="data" label="GFF File"/>
        <param name="mode" type="select" label="Mode">
            <help>Mode to handle reads overlapping more than one feature.</help>
            <option value="union" selected="true">Union</option>
            <option value="intersection-strict">Intersection (strict)</option>
            <option value="intersection-nonempty">Intersection (nonempty)</option>
        </param>
        <param name="stranded" type="select" label="Stranded">
            <help>Specify whether the data is from a strand-specific assay. 'Reverse' means yes with reversed strand interpretation.</help>
            <option value="yes" selected="true">Yes</option>
            <option value="no">No</option>
            <option value="reverse">Reverse</option>
        </param>
        <param name="minaqual" type="integer" value="10" label="Minimum alignment quality">
            <help>Skip all reads with alignment quality lower than the given minimum value</help>
        </param>
        <param name="featuretype" type="text" value="exon" label="Feature type">
            <help>Feature type (3rd column in GFF file) to be used. All features of other types are ignored. The default, suitable for RNA-Seq and Ensembl GTF files, is exon.</help>
        </param>
        <param name="idattr" type="text" value="gene_id" label="ID Attribute">
            <help>GFF attribute to be used as feature ID. Several GFF lines with the same feature ID will be considered as parts of the same feature. The feature ID is used to identity the counts in the output table. All features of the specified type MUST have a value for this attribute. The default, suitable for RNA-Seq and Ensembl GTF files, is gene_id.</help>
        </param>
        <conditional name="samout_conditional">

21:19e80c732870

            #set $reference_fasta_filename = str( $samout_conditional.reference_source.ref_file.fields.path )
        #end if
    #end if
    #if $force_sort:
        #if $samfile.extension == 'bam':
            samtools sort -n "$samfile" "name_sorted_alignment" &&
        #else
            samtools view -Su -t ${reference_fasta_filename}.fai "$samfile" | samtools sort -n - "name_sorted_alignment" &&
        #end if
    #end if
    htseq-count
    --mode=$mode
    --stranded=$stranded
    --minaqual=$minaqual
    --type="$featuretype"
    --idattr="$idattr"
    #if $samout_conditional.samout:
        --samout=$__new_file_path__/${samoutfile.id}_tmp
    #end if
    #if $force_sort:
        --order=name

    #else
        --order=pos
        --format=$samfile.extension
        $samfile
    #end if
    "$gfffile"
    | awk '{if ($1 ~ "no_feature|ambiguous|too_low_aQual|not_aligned|alignment_not_unique") print $0 | "cat 1>&2"; else print $0}' > $counts 2>$othercounts
    #if $samout_conditional.samout:
        && samtools view -Su -t ${reference_fasta_filename}.fai $__new_file_path__/${samoutfile.id}_tmp | samtools sort -o - sorted > $samoutfile
    #end if
    ]]>
    </command>

    <inputs>
        <param format="sam,bam" name="samfile" type="data" label="Aligned SAM/BAM File"/>
        <param format="gff" name="gfffile" type="data" label="GFF File"/>
        <param name="mode" type="select" label="Mode" help="(--mode)">
            <help>Mode to handle reads overlapping more than one feature.</help>
            <option value="union" selected="true">Union</option>
            <option value="intersection-strict">Intersection (strict)</option>
            <option value="intersection-nonempty">Intersection (nonempty)</option>
        </param>
        <param name="stranded" type="select" label="Stranded" help="(--stranded)">
            <help>Specify whether the data is from a strand-specific assay. 'Reverse' means yes with reversed strand interpretation.</help>
            <option value="yes" selected="true">Yes</option>
            <option value="no">No</option>
            <option value="reverse">Reverse</option>
        </param>
        <param name="minaqual" type="integer" value="10" label="Minimum alignment quality">
            <help>Skip all reads with alignment quality lower than the given minimum value. (-minaqual)</help>
        </param>
        <param name="featuretype" type="text" value="exon" label="Feature type">
            <help>Feature type (3rd column in GFF file) to be used. All features of other types are ignored. The default, suitable for RNA-Seq and Ensembl GTF files, is exon. (--type)</help>
        </param>
        <param name="idattr" type="text" value="gene_id" label="ID Attribute">
            <help>GFF attribute to be used as feature ID. Several GFF lines with the same feature ID will be considered as parts of the same feature. The feature ID is used to identity the counts in the output table. All features of the specified type MUST have a value for this attribute. The default, suitable for RNA-Seq and Ensembl GTF files, is gene_id.</help>
        </param>
        <conditional name="samout_conditional">