comparison sam_dump.xml @ 1:462ee06c9358 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/sra-tools commit 4defaa3ff1c21e2ec39033bfe63ee69471104ede

0:b723c120161a

<tool id="sam_dump" name="Extract reads" version="1.2.5">
    <description>in SAM or BAM format from NCBI SRA.</description>
    <macros>
        <import>sra_macros.xml</import>
    </macros>
    <expand macro="requirements"/>
    <version_command>sam-dump --version</version_command>
    <command>
        <![CDATA[
        ## Need to set the home directory to the current working directory,
        ## else the tool tries to write to home/.ncbi and fails when used 
        ## with a cluster manager. 
        export HOME=\$PWD &&
        vdb-config --restore-defaults &&

        ## Do not use prefetch if region is specified, to avoid downloading
        ## the complete sra file.
        #if ( str( $adv.region ) == "" ):
            ASCP_PATH=`command -v ascp` &&
            ASCP_KEY=`dirname \$ASCP_PATH`/asperaweb_id_dsa.openssh || true &&
            prefetch --ascp-path "\$ASCP_PATH|\$ASCP_KEY" "$input.accession" &&
            ## Duplicate vdb-config, in case settings changed between prefetch and
            ## dump command.
            vdb-config -s "/repository/user/main/public/root=\$PWD" &&
        #end if
        sam-dump --log-level fatal --disable-multithreading

            --primary
        #end if
        #if $input.input_select == "file":
            "$input.file"
        #elif $input.input_select == "accession_number":
            "$input.accession"
        #end if
        #if str( $outputformat ) == "bam":
            | samtools view -Sb - 2> /dev/null
        #end if
        #if $input.input_select == "file":
            > "$output_file"
        #elif $input.input_select == "accession_number":
            > "$output_accession"
        #end if
        ]]>
    </command>
    <inputs>
        <expand macro="input_conditional"/>
        <param name="outputformat" type="select" label="select output format">

            </param>
            <expand macro="minMapq"></expand>
        </section>
    </inputs>
    <outputs>
        <data name="output_accession" format="bam" label="${input.accession}.${outputformat}">
            <filter>input['input_select'] == "accession_number"</filter>
            <change_format>
                <when input="outputformat" value="sam" format="sam"/>
            </change_format>

1:462ee06c9358

<tool id="sam_dump" name="Extract reads" version="1.3.0">
    <description>in SAM or BAM format from NCBI SRA.</description>
    <macros>
        <import>sra_macros.xml</import>
    </macros>
    <expand macro="requirements"/>
    <version_command>sam-dump --version</version_command>
    <command>
        <![CDATA[
    #if $input.input_select=="file_list":
      for acc in `cat $input.file_list` ;
      do
    #elif $input.input_select=="accession_number":
      acc="$input.accession" &&
    #end if

    #if $input.input_select=="file_list" or $input.input_select=="accession_number":
          [ ""\$acc" =~ ^[E|S|D]RR[0-9]{1,}$" ] && (
    #end if



        ## Need to set the home directory to the current working directory,
        ## else the tool tries to write to home/.ncbi and fails when used 
        ## with a cluster manager. 
        export HOME=\$PWD &&
        vdb-config --restore-defaults &&

        ## Do not use prefetch if region is specified, to avoid downloading
        ## the complete sra file.
        #if ( str( $adv.region ) == "" ):
            ASCP_PATH=`command -v ascp` &&
            ASCP_KEY=`dirname \$ASCP_PATH`/asperaweb_id_dsa.openssh || true &&
            prefetch --ascp-path "\$ASCP_PATH|\$ASCP_KEY" "\$acc" &&
            ## Duplicate vdb-config, in case settings changed between prefetch and
            ## dump command.
            vdb-config -s "/repository/user/main/public/root=\$PWD" &&
        #end if
        sam-dump --log-level fatal --disable-multithreading

            --primary
        #end if
        #if $input.input_select == "file":
            "$input.file"
        #elif $input.input_select == "accession_number":
            "\$acc"
        #elif $input.input_select=="file_list":
            "\$acc"
        #end if
        
        #if str( $outputformat ) == "bam":
            | samtools view -Sb - 2> /dev/null
        #end if
        #if $input.input_select == "file":
            > "$output_file"
        #elif $input.input_select == "accession_number":
            > "$output_accession" )
        #end if

        #if $input.input_select=="file_list":
                 #if str( $outputformat ) == "bam":
                      > "\$acc.bam"
                 #elif str( $outputformat ) == "sam":
                      > "\$acc.sam"
                 #end if
        ) ; done
        #end if


        ]]>
    </command>
    <inputs>
        <expand macro="input_conditional"/>
        <param name="outputformat" type="select" label="select output format">

            </param>
            <expand macro="minMapq"></expand>
        </section>
    </inputs>
    <outputs>
        <collection name="output_collection" type='list'>
          <filter>input['input_select'] == "file_list"</filter>
          <discover_datasets pattern="(?P&lt;designation&gt;.+)\.bam" directory="." ext='bam'/>
          <discover_datasets pattern="(?P&lt;designation&gt;.+)\.sam" directory="." ext='sam'/>
        </collection>      
        <data name="output_accession" format="bam" label="${input.accession}.${outputformat}">
            <filter>input['input_select'] == "accession_number"</filter>
            <change_format>
                <when input="outputformat" value="sam" format="sam"/>
            </change_format>