# HG changeset patch # User brigidar # Date 1446510663 18000 # Node ID a928631d829741284d3956c602e1140172ae94bc # Parent ea2f686dfd4a3861cd0d683b06a5350febe4444c new wrapper diff -r ea2f686dfd4a -r a928631d8297 vcf_to_snp.xml --- a/vcf_to_snp.xml Mon Nov 02 19:30:43 2015 -0500 +++ b/vcf_to_snp.xml Mon Nov 02 19:31:03 2015 -0500 @@ -1,4 +1,4 @@ - + vcf to snp for fasta @@ -12,11 +12,7 @@ vcf_snp.py -s ${vcf} - -o output.txt; - cut -f1 output.txt > header.txt | cut -f2- output.txt > seq.txt; - cat seq.txt | sed 's/[[:blank:]]//g' >seq2.txt; - paste header.txt seq2.txt >tab.txt; - + -o output.txt @@ -40,7 +36,7 @@ **What it does** -Transforms a vcf into a snp tab file +Transforms a vcf into a snp tab file. Has to take the nucleotide from the reference column when identical. -------------