changeset 1:82f2b6f20fa2 draft

planemo upload commit c295a06009ccafd0c1dc8aba5dd87239795d2b61
author iuc
date Sat, 26 Jan 2019 14:36:48 -0500
parents c9a8ef2aa380
children 776ebd1239da
files snippy.xml test-data/fna_ref/snps.gff test-data/fna_ref/snps.txt test-data/gbk_ref/snps.gff test-data/gbk_ref/snps.txt test-data/map_qual/snps.gff test-data/map_qual/snps.txt test-data/min_cov/snps.gff test-data/min_cov/snps.txt test-data/min_frac/snps.gff test-data/min_frac/snps.txt
diffstat 11 files changed, 107 insertions(+), 83 deletions(-) [+]
line wrap: on
line diff
--- a/snippy.xml	Fri Feb 16 13:40:16 2018 -0500
+++ b/snippy.xml	Sat Jan 26 14:36:48 2019 -0500
@@ -1,4 +1,4 @@
-<tool id="snippy" name="snippy" version="@VERSION@">
+<tool id="snippy" name="snippy" version="@VERSION@+galaxy1">
   <description>
       Snippy finds SNPs between a haploid reference genome and your NGS sequence reads.
   </description>
@@ -40,14 +40,14 @@
                 --pe2 '$fastq_input.fastq_input2'
             #end if
             #if str( $fastq_input.fastq_input_selector ) == "paired_collection"
-                --pe1 '$fastq_input.fastq_input1.forward'
-                --pe2 '$fastq_input.fastq_input1.reverse'
+                --pe1 '$fastq_input.fastq_input.forward'
+                --pe2 '$fastq_input.fastq_input.reverse'
             #end if
             #if str( $fastq_input.fastq_input_selector ) == "single"
-                --se '$fastq_input.fastq_input1'
+                --se '$fastq_input.fastq_input'
             #end if
             #if str( $fastq_input.fastq_input_selector ) == "paired_iv"
-                --peil '$fastq_input.fastq_input1'
+                --peil '$fastq_input.fastq_input'
             #end if
 
         &&
@@ -57,9 +57,13 @@
         &&
 
         #import re
-        #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input1.element_identifier)
-        mkdir -p ${dir_name}/reference && cp out/snps.tab out/snps.aligned.fa ${dir_name}/ && cp out/reference/ref.fa ${dir_name}/reference/ &&
+        #if str( $fastq_input.fastq_input_selector ) == "paired"
+            #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input1.element_identifier)
+        #else
+            #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input.element_identifier)
+        #end if
 
+        mkdir -p ${dir_name} && cp -r out/reference out/snps.tab out/snps.aligned.fa ${dir_name}/ &&
         tar -czf out.tgz ${dir_name}
 
 
@@ -81,13 +85,13 @@
                 <param name="fastq_input2" type="data" format="fastqsanger,fasta" label="Select second set of reads" help="Specify dataset with reverse reads"/>
             </when>
             <when value="single">
-                <param name="fastq_input1" type="data" format="fastqsanger,fasta" label="Select fastq dataset" help="Specify dataset with single reads"/>
+                <param name="fastq_input" type="data" format="fastqsanger,fasta" label="Select fastq dataset" help="Specify dataset with single reads"/>
             </when>
             <when value="paired_collection">
-                <param name="fastq_input1" format="fastqsanger,fasta" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
+                <param name="fastq_input" format="fastqsanger,fasta" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
             </when>
             <when value="paired_iv">
-                <param name="fastq_input1" type="data" format="fastqsanger" label="Select fastq dataset" help="Specify dataset with interleaved reads"/>
+                <param name="fastq_input" type="data" format="fastqsanger" label="Select fastq dataset" help="Specify dataset with interleaved reads"/>
             </when>
         </conditional>
 
@@ -204,6 +208,20 @@
             <output name="snpgff" ftype="gff3" file="min_frac/snps.gff" />
         </test>
 
+        <test> <!-- test 6 - fasta ref default paired_collection -->
+	    <param name="ref" value="wildtype.fna" ftype="fasta" />
+            <param name="fastq_input_selector" value="paired_collection" />
+            <param name="fastq_input">
+                <collection type="paired">
+                    <element name="forward" ftype="fastqsanger" value="mutant_R1.fastq" />
+                    <element name="reverse" ftype="fastqsanger" value="mutant_R2.fastq" />
+                </collection>
+            </param>
+            <param name="outputs" value="outgff,outsum" />
+            <output name="snpsum" ftype="tabular" file="fna_ref/snps.txt" lines_diff="6" />
+            <output name="snpgff" ftype="gff3" file="fna_ref/snps.gff" />
+        </test>
+
     </tests>
 
 
--- a/test-data/fna_ref/snps.gff	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/fna_ref/snps.gff	Sat Jan 26 14:36:48 2019 -0500
@@ -1,16 +1,17 @@
 ##gff-version 3
 Wildtype	snippy:3.2-dev	variation	24388	24388	.	.	0	note=snp A=>G G:22 A:0
 Wildtype	snippy:3.2-dev	variation	29479	29479	.	.	0	note=snp T=>G G:21 T:0
-Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:24 T:0
+Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:25 T:0
 Wildtype	snippy:3.2-dev	variation	102969	102969	.	.	0	note=snp G=>C C:16 G:0
-Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:20 T:0
-Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAA=>GA GA:11 GAA:0
-Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:21 T:0
-Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:16 T:0
-Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATT=>AT AT:25 ATT:0
-Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:18 GT:0
+Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:21 T:0
+Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAAAATAG=>GAAATAG GAAATAG:11 GAAAATAG:0
+Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:22 T:0
+Wildtype	snippy:3.2-dev	variation	106722	106722	.	.	0	note=ins GTTTTTAAC=>GTTTTTTAAC GTTTTTTAAC:11 GTTTTTAAC:0
+Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:17 T:0
+Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATTGTGA=>ATGTGA ATGTGA:26 ATTGTGA:0
+Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:19 GT:0
 Wildtype	snippy:3.2-dev	variation	138877	138877	.	.	0	note=snp G=>C C:14 G:0
 Wildtype	snippy:3.2-dev	variation	138920	138920	.	.	0	note=snp A=>G G:10 A:0
+Wildtype	snippy:3.2-dev	variation	138954	138954	.	.	0	note=snp A=>G G:10 A:0
 Wildtype	snippy:3.2-dev	variation	160547	160547	.	.	0	note=del GTC=>GC GC:18 GTC:0
-Wildtype	snippy:3.2-dev	variation	160552	160552	.	.	0	note=del CTA=>CA CA:20 CTA:0
-Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTT=>GT GT:18 GTT:0
+Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTTTTCTTTAA=>GTTTCTTTAA GTTTCTTTAA:17 GTTTTCTTTAA:0
--- a/test-data/fna_ref/snps.txt	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/fna_ref/snps.txt	Sat Jan 26 14:36:48 2019 -0500
@@ -1,9 +1,10 @@
-DateTime	2017-07-12T04:51:24
-ReadFiles	/home/ubuntu/snippy/test-data/mutant_R1.fastq /home/ubuntu/snippy/test-data/mutant_R2.fastq
-Reference	/home/ubuntu/snippy/test-data/wildtype.fna
+DateTime	2019-01-24T13:08:34
+ReadFiles	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R1.fastq /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R2.fastq
+Reference	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/wildtype.fna
 ReferenceSize	200694
 Software	snippy 3.2-dev
-Variant-DEL	5
+Variant-DEL	4
+Variant-INS	1
 Variant-MNP	1
-Variant-SNP	9
-VariantTotal	15
+Variant-SNP	10
+VariantTotal	16
--- a/test-data/gbk_ref/snps.gff	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/gbk_ref/snps.gff	Sat Jan 26 14:36:48 2019 -0500
@@ -1,16 +1,17 @@
 ##gff-version 3
 Wildtype	snippy:3.2-dev	variation	24388	24388	.	.	0	note=snp A=>G G:22 A:0 CDS + 1/702 1/233 initiator_codon_variant c.1A>G p.Met1? WILD_00022 walR Transcriptional regulatory protein WalR
 Wildtype	snippy:3.2-dev	variation	29479	29479	.	.	0	note=snp T=>G G:21 T:0 CDS + 39/792 13/263 synonymous_variant c.39T>G p.Gly13Gly WILD_00026 yycJ Putative metallo-hydrolase YycJ
-Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:24 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
+Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:25 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
 Wildtype	snippy:3.2-dev	variation	102969	102969	.	.	0	note=snp G=>C C:16 G:0 CDS - 87/1281 29/426 synonymous_variant c.87C>G p.Gly29Gly WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:20 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAA=>GA GA:11 GAA:0
-Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:21 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
-Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:16 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
-Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATT=>AT AT:25 ATT:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
-Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:18 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
+Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:21 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
+Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAAAATAG=>GAAATAG GAAATAG:11 GAAAATAG:0
+Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:22 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	106722	106722	.	.	0	note=ins GTTTTTAAC=>GTTTTTTAAC GTTTTTTAAC:11 GTTTTTAAC:0 CDS - 867/993 289/330 frameshift_variant c.867_868insA p.Asn291fs WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:17 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
+Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATTGTGA=>ATGTGA ATGTGA:26 ATTGTGA:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
+Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:19 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
 Wildtype	snippy:3.2-dev	variation	138877	138877	.	.	0	note=snp G=>C C:14 G:0 CDS + 1119/1545 373/514 missense_variant c.1119G>C p.Trp373Cys WILD_00125  hypothetical protein
 Wildtype	snippy:3.2-dev	variation	138920	138920	.	.	0	note=snp A=>G G:10 A:0 CDS + 1162/1545 388/514 missense_variant c.1162A>G p.Lys388Glu WILD_00125  hypothetical protein
+Wildtype	snippy:3.2-dev	variation	138954	138954	.	.	0	note=snp A=>G G:10 A:0 CDS + 1196/1545 399/514 missense_variant c.1196A>G p.Asp399Gly WILD_00125  hypothetical protein
 Wildtype	snippy:3.2-dev	variation	160547	160547	.	.	0	note=del GTC=>GC GC:18 GTC:0
-Wildtype	snippy:3.2-dev	variation	160552	160552	.	.	0	note=del CTA=>CA CA:20 CTA:0
-Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTT=>GT GT:18 GTT:0 CDS - 28/1356 10/451 frameshift_variant c.28delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
+Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTTTTCTTTAA=>GTTTCTTTAA GTTTCTTTAA:17 GTTTTCTTTAA:0 CDS - 26/1356 9/451 frameshift_variant c.26delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
--- a/test-data/gbk_ref/snps.txt	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/gbk_ref/snps.txt	Sat Jan 26 14:36:48 2019 -0500
@@ -1,9 +1,10 @@
-DateTime	2017-07-12T04:52:49
-ReadFiles	/home/ubuntu/snippy/test-data/mutant_R1.fastq /home/ubuntu/snippy/test-data/mutant_R2.fastq
-Reference	/home/ubuntu/snippy/test-data/wildtype.gbk
+DateTime	2019-01-24T13:29:29
+ReadFiles	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R1.fastq /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R2.fastq
+Reference	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/wildtype.gbk
 ReferenceSize	200727
 Software	snippy 3.2-dev
-Variant-DEL	5
+Variant-DEL	4
+Variant-INS	1
 Variant-MNP	1
-Variant-SNP	9
-VariantTotal	15
+Variant-SNP	10
+VariantTotal	16
--- a/test-data/map_qual/snps.gff	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/map_qual/snps.gff	Sat Jan 26 14:36:48 2019 -0500
@@ -1,16 +1,17 @@
 ##gff-version 3
 Wildtype	snippy:3.2-dev	variation	24388	24388	.	.	0	note=snp A=>G G:22 A:0 CDS + 1/702 1/233 initiator_codon_variant c.1A>G p.Met1? WILD_00022 walR Transcriptional regulatory protein WalR
 Wildtype	snippy:3.2-dev	variation	29479	29479	.	.	0	note=snp T=>G G:21 T:0 CDS + 39/792 13/263 synonymous_variant c.39T>G p.Gly13Gly WILD_00026 yycJ Putative metallo-hydrolase YycJ
-Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:24 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
+Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:25 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
 Wildtype	snippy:3.2-dev	variation	102969	102969	.	.	0	note=snp G=>C C:16 G:0 CDS - 87/1281 29/426 synonymous_variant c.87C>G p.Gly29Gly WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:20 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAA=>GA GA:11 GAA:0
-Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:21 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
-Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:16 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
-Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATT=>AT AT:25 ATT:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
-Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:18 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
+Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:21 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
+Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAAAATAG=>GAAATAG GAAATAG:11 GAAAATAG:0
+Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:22 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	106722	106722	.	.	0	note=ins GTTTTTAAC=>GTTTTTTAAC GTTTTTTAAC:11 GTTTTTAAC:0 CDS - 867/993 289/330 frameshift_variant c.867_868insA p.Asn291fs WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:17 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
+Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATTGTGA=>ATGTGA ATGTGA:26 ATTGTGA:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
+Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:19 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
 Wildtype	snippy:3.2-dev	variation	138877	138877	.	.	0	note=snp G=>C C:14 G:0 CDS + 1119/1545 373/514 missense_variant c.1119G>C p.Trp373Cys WILD_00125  hypothetical protein
 Wildtype	snippy:3.2-dev	variation	138920	138920	.	.	0	note=snp A=>G G:10 A:0 CDS + 1162/1545 388/514 missense_variant c.1162A>G p.Lys388Glu WILD_00125  hypothetical protein
+Wildtype	snippy:3.2-dev	variation	138954	138954	.	.	0	note=snp A=>G G:10 A:0 CDS + 1196/1545 399/514 missense_variant c.1196A>G p.Asp399Gly WILD_00125  hypothetical protein
 Wildtype	snippy:3.2-dev	variation	160547	160547	.	.	0	note=del GTC=>GC GC:18 GTC:0
-Wildtype	snippy:3.2-dev	variation	160552	160552	.	.	0	note=del CTA=>CA CA:20 CTA:0
-Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTT=>GT GT:18 GTT:0 CDS - 28/1356 10/451 frameshift_variant c.28delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
+Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTTTTCTTTAA=>GTTTCTTTAA GTTTCTTTAA:17 GTTTTCTTTAA:0 CDS - 26/1356 9/451 frameshift_variant c.26delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
--- a/test-data/map_qual/snps.txt	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/map_qual/snps.txt	Sat Jan 26 14:36:48 2019 -0500
@@ -1,9 +1,10 @@
-DateTime	2017-07-13T05:37:56
-ReadFiles	/home/ubuntu/snippy/test-data/mutant_R1.fastq /home/ubuntu/snippy/test-data/mutant_R2.fastq
-Reference	/home/ubuntu/snippy/test-data/wildtype.gbk
+DateTime	2019-01-24T14:57:31
+ReadFiles	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R1.fastq /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R2.fastq
+Reference	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/wildtype.gbk
 ReferenceSize	200727
 Software	snippy 3.2-dev
-Variant-DEL	5
+Variant-DEL	4
+Variant-INS	1
 Variant-MNP	1
-Variant-SNP	9
-VariantTotal	15
+Variant-SNP	10
+VariantTotal	16
--- a/test-data/min_cov/snps.gff	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/min_cov/snps.gff	Sat Jan 26 14:36:48 2019 -0500
@@ -1,13 +1,12 @@
 ##gff-version 3
 Wildtype	snippy:3.2-dev	variation	24388	24388	.	.	0	note=snp A=>G G:22 A:0 CDS + 1/702 1/233 initiator_codon_variant c.1A>G p.Met1? WILD_00022 walR Transcriptional regulatory protein WalR
 Wildtype	snippy:3.2-dev	variation	29479	29479	.	.	0	note=snp T=>G G:21 T:0 CDS + 39/792 13/263 synonymous_variant c.39T>G p.Gly13Gly WILD_00026 yycJ Putative metallo-hydrolase YycJ
-Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:24 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
+Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:25 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
 Wildtype	snippy:3.2-dev	variation	102969	102969	.	.	0	note=snp G=>C C:16 G:0 CDS - 87/1281 29/426 synonymous_variant c.87C>G p.Gly29Gly WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:20 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:21 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
-Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:16 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
-Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATT=>AT AT:25 ATT:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
-Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:18 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
+Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:21 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
+Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:22 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:17 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
+Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATTGTGA=>ATGTGA ATGTGA:26 ATTGTGA:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
+Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:19 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
 Wildtype	snippy:3.2-dev	variation	160547	160547	.	.	0	note=del GTC=>GC GC:18 GTC:0
-Wildtype	snippy:3.2-dev	variation	160552	160552	.	.	0	note=del CTA=>CA CA:20 CTA:0
-Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTT=>GT GT:18 GTT:0 CDS - 28/1356 10/451 frameshift_variant c.28delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
+Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTTTTCTTTAA=>GTTTCTTTAA GTTTCTTTAA:17 GTTTTCTTTAA:0 CDS - 26/1356 9/451 frameshift_variant c.26delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
--- a/test-data/min_cov/snps.txt	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/min_cov/snps.txt	Sat Jan 26 14:36:48 2019 -0500
@@ -1,9 +1,9 @@
-DateTime	2017-07-12T04:54:55
-ReadFiles	/home/ubuntu/snippy/test-data/mutant_R1.fastq /home/ubuntu/snippy/test-data/mutant_R2.fastq
-Reference	/home/ubuntu/snippy/test-data/wildtype.gbk
+DateTime	2019-01-24T15:21:35
+ReadFiles	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R1.fastq /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R2.fastq
+Reference	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/wildtype.gbk
 ReferenceSize	200727
 Software	snippy 3.2-dev
-Variant-DEL	4
+Variant-DEL	3
 Variant-MNP	1
 Variant-SNP	7
-VariantTotal	12
+VariantTotal	11
--- a/test-data/min_frac/snps.gff	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/min_frac/snps.gff	Sat Jan 26 14:36:48 2019 -0500
@@ -1,17 +1,18 @@
 ##gff-version 3
 Wildtype	snippy:3.2-dev	variation	24388	24388	.	.	0	note=snp A=>G G:22 A:0 CDS + 1/702 1/233 initiator_codon_variant c.1A>G p.Met1? WILD_00022 walR Transcriptional regulatory protein WalR
 Wildtype	snippy:3.2-dev	variation	29479	29479	.	.	0	note=snp T=>G G:21 T:0 CDS + 39/792 13/263 synonymous_variant c.39T>G p.Gly13Gly WILD_00026 yycJ Putative metallo-hydrolase YycJ
-Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:24 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
+Wildtype	snippy:3.2-dev	variation	47299	47299	.	.	0	note=snp T=>A A:25 T:0 CDS + 54/1758 18/585 stop_gained c.54T>A p.Cys18* WILD_00043 mecR1 Methicillin resistance mecR1 protein
 Wildtype	snippy:3.2-dev	variation	102969	102969	.	.	0	note=snp G=>C C:16 G:0 CDS - 87/1281 29/426 synonymous_variant c.87C>G p.Gly29Gly WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:20 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
-Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAA=>GA GA:11 GAA:0
-Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:21 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
-Wildtype	snippy:3.2-dev	variation	106722	106722	.	.	0	note=ins GTTTTTA=>GTTTTTTA GTTTTTTA:11 GTTTTTA:0 CDS - 867/993 289/330 frameshift_variant c.867_868insA p.Asn291fs WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
-Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:16 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
-Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATT=>AT AT:25 ATT:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
-Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:18 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
+Wildtype	snippy:3.2-dev	variation	103048	103048	.	.	0	note=snp T=>A A:21 T:0 CDS - 8/1281 3/426 missense_variant c.8A>T p.Lys3Met WILD_00093 spa Immunoglobulin G-binding protein A
+Wildtype	snippy:3.2-dev	variation	103379	103379	.	.	0	note=del GAAAATAG=>GAAATAG GAAATAG:11 GAAAATAG:0
+Wildtype	snippy:3.2-dev	variation	106602	106602	.	.	0	note=snp T=>G G:22 T:0 CDS - 993/993 331/330 stop_lost&splice_region_variant c.993A>C p.Ter331Tyrext*? WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	106722	106722	.	.	0	note=ins GTTTTTAAC=>GTTTTTTAAC GTTTTTTAAC:11 GTTTTTAAC:0 CDS - 867/993 289/330 frameshift_variant c.867_868insA p.Asn291fs WILD_00097 yfiY putative siderophore-binding lipoprotein YfiY
+Wildtype	snippy:3.2-dev	variation	109833	109833	.	.	0	note=snp T=>A A:17 T:0 CDS + 1/1755 1/584 initiator_codon_variant c.1T>A p.Leu1? WILD_00100 iucC_1 Aerobactin synthase
+Wildtype	snippy:3.2-dev	variation	114540	114540	.	.	0	note=del ATTGTGA=>ATGTGA ATGTGA:26 ATTGTGA:0 CDS + 1717/1737 573/578 frameshift_variant c.1717delT p.Cys573fs WILD_00102 iucA N(2)-citryl-N(6)-acetyl-N(6)-hydroxylysine synthase
+Wildtype	snippy:3.2-dev	variation	129881	129881	.	.	0	note=mnp GT=>AA AA:19 GT:0 CDS + 55/708 19/235 missense_variant c.55_56delGTinsAA p.Val19Asn WILD_00117 deoD Purine nucleoside phosphorylase DeoD-type
 Wildtype	snippy:3.2-dev	variation	138877	138877	.	.	0	note=snp G=>C C:14 G:0 CDS + 1119/1545 373/514 missense_variant c.1119G>C p.Trp373Cys WILD_00125  hypothetical protein
 Wildtype	snippy:3.2-dev	variation	138920	138920	.	.	0	note=snp A=>G G:10 A:0 CDS + 1162/1545 388/514 missense_variant c.1162A>G p.Lys388Glu WILD_00125  hypothetical protein
+Wildtype	snippy:3.2-dev	variation	138954	138954	.	.	0	note=snp A=>G G:10 A:0 CDS + 1196/1545 399/514 missense_variant c.1196A>G p.Asp399Gly WILD_00125  hypothetical protein
 Wildtype	snippy:3.2-dev	variation	160547	160547	.	.	0	note=del GTC=>GC GC:18 GTC:0
-Wildtype	snippy:3.2-dev	variation	160552	160552	.	.	0	note=del CTA=>CA CA:20 CTA:0
-Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTT=>GT GT:18 GTT:0 CDS - 28/1356 10/451 frameshift_variant c.28delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
+Wildtype	snippy:3.2-dev	variation	160552	160552	.	.	0	note=del CTATTA=>CATTA CATTA:18 CTATTA:0
+Wildtype	snippy:3.2-dev	variation	190866	190866	.	.	0	note=del GTTTTCTTTAA=>GTTTCTTTAA GTTTCTTTAA:17 GTTTTCTTTAA:0 CDS - 26/1356 9/451 frameshift_variant c.26delA p.Asn10fs WILD_00166 brnQ Branched-chain amino acid transport system 2 carrier protein
--- a/test-data/min_frac/snps.txt	Fri Feb 16 13:40:16 2018 -0500
+++ b/test-data/min_frac/snps.txt	Sat Jan 26 14:36:48 2019 -0500
@@ -1,10 +1,10 @@
-DateTime	2017-07-12T04:56:36
-ReadFiles	/home/ubuntu/snippy/test-data/mutant_R1.fastq /home/ubuntu/snippy/test-data/mutant_R2.fastq
-Reference	/home/ubuntu/snippy/test-data/wildtype.gbk
+DateTime	2019-01-24T15:35:47
+ReadFiles	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R1.fastq /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R2.fastq
+Reference	/home/dfornika/Code/tools-iuc/tools/snippy/test-data/wildtype.gbk
 ReferenceSize	200727
 Software	snippy 3.2-dev
 Variant-DEL	5
 Variant-INS	1
 Variant-MNP	1
-Variant-SNP	9
-VariantTotal	16
+Variant-SNP	10
+VariantTotal	17