Mercurial > repos > greg > cfsan_snp_pipeline_call_sites

changeset 3:5bfc9e1ce9e9 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Thu, 23 Nov 2023 19:33:24 +0000
parents 991fbae39ee6
children
files cfsan_snp_pipeline_call_sites.xml
diffstat 1 files changed, 5 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/cfsan_snp_pipeline_call_sites.xml	Wed Nov 22 14:10:59 2023 +0000
+++ b/cfsan_snp_pipeline_call_sites.xml	Thu Nov 23 19:33:24 2023 +0000
@@ -7,11 +7,13 @@
     <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' &&
+    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="bam"/>
+        <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>
@@ -22,7 +24,7 @@
     <tests>
         <test>
             <param name="reference" value="lambda_virus.fasta"/>
-            <param name="input" value="reads.bam"/>
+            <param name="input" value="reads.sam"/>
             <output name="calls" ftype="vcf">
                 <assert_contents>
                     <has_size value="9976" delta="100"/>