comparison cfsan_snp_pipeline_call_sites.xml @ 3:5bfc9e1ce9e9 draft default tip

Uploaded

2:991fbae39ee6

    </macros>
    <expand macro="requirements"/>
    <command detect_errors="exit_code"><![CDATA[
    cp '$reference' 'reference.fasta' &&
    samtools faidx 'reference.fasta' &&  
    picard MarkDuplicates INPUT='$input' OUTPUT=/dev/stdout METRICS_FILE='$metrics' | samtools mpileup -f ./reference.fasta - -o '$pileup' &&
    varscan mpileup2snp '$pileup' --min_var-freq 0.90 --output-vcf 1 > '$calls' 
    ]]></command>
    <inputs>
        <param name="input" label="Read alignment to reference" type="data" format="bam"/>
        <param name="reference" label="FASTA Reference from your history" type="data" format="fasta"/>
    </inputs>
    <outputs>
        <data name="calls" label="Call Sites (unfiltered SNPs)" format="vcf"/>
        <data name="pileup" label="Call Sites (pileup)" format="pileup"/>
        <data name="metrics" label="Call sites (Picard metrics)" format="txt"/>
    </outputs>
    <tests>
        <test>
            <param name="reference" value="lambda_virus.fasta"/>
            <param name="input" value="reads.bam"/>
            <output name="calls" ftype="vcf">
                <assert_contents>
                    <has_size value="9976" delta="100"/>
                </assert_contents>
            </output>

3:5bfc9e1ce9e9

    </macros>
    <expand macro="requirements"/>
    <command detect_errors="exit_code"><![CDATA[
    cp '$reference' 'reference.fasta' &&
    samtools faidx 'reference.fasta' &&  
    samtools view -b -F 4 '$input' -o ./unsorted.bam &&
    samtools sort -o ./sorted.bam -O BAM ./unsorted.bam &&
    picard MarkDuplicates INPUT=./sorted.bam OUTPUT=/dev/stdout METRICS_FILE='$metrics' | samtools mpileup -f ./reference.fasta - -o '$pileup' &&
    varscan mpileup2snp '$pileup' --min_var-freq 0.90 --output-vcf 1 > '$calls' 
    ]]></command>
    <inputs>
        <param name="input" label="Read alignment to reference" type="data" format="sam"/>
        <param name="reference" label="FASTA Reference from your history" type="data" format="fasta"/>
    </inputs>
    <outputs>
        <data name="calls" label="Call Sites (unfiltered SNPs)" format="vcf"/>
        <data name="pileup" label="Call Sites (pileup)" format="pileup"/>
        <data name="metrics" label="Call sites (Picard metrics)" format="txt"/>
    </outputs>
    <tests>
        <test>
            <param name="reference" value="lambda_virus.fasta"/>
            <param name="input" value="reads.sam"/>
            <output name="calls" ftype="vcf">
                <assert_contents>
                    <has_size value="9976" delta="100"/>
                </assert_contents>
            </output>