# HG changeset patch # User greg # Date 1628018690 0 # Node ID be5875f29ea48b2a62561b978add67c3992d194f # Parent 0fe292b20b9dab01a77df67a6847e90e5f16b2b4 Uploaded diff -r 0fe292b20b9d -r be5875f29ea4 .shed.yml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/.shed.yml Tue Aug 03 19:24:50 2021 +0000 @@ -0,0 +1,11 @@ +name: vsnp_get_snps +owner: greg +description: | + Contains a tool that collects quality parsimonious SNPs from vcf files and outputs alignment files in fasta format. +homepage_url: https://github.com/USDA-VS/vSNP +long_description: | + Contains a tool that collects quality parsimonious SNPs from vcf files and outputs alignment files in fasta format. +remote_repository_url: https://github.com/gregvonkuster/galaxy_tools/tree/master/tools/sequence_analysis/vsnp/vsnp_get_snps +type: unrestricted +categories: + - Sequence Analysis diff -r 0fe292b20b9d -r be5875f29ea4 macros.xml --- a/macros.xml Thu Jul 29 13:16:03 2021 +0000 +++ b/macros.xml Tue Aug 03 19:24:50 2021 +0000 @@ -1,7 +1,34 @@ 1.0 - 19.09 + 20.09 + + biopython + + + numpy + + + openpyxl + + + pandas + + + pysam + + + pyvcf + + + pyyaml + + + xlrd + + + xlsxwriter + diff -r 0fe292b20b9d -r be5875f29ea4 vsnp_get_snps.py --- a/vsnp_get_snps.py Thu Jul 29 13:16:03 2021 +0000 +++ b/vsnp_get_snps.py Tue Aug 03 19:24:50 2021 +0000 @@ -440,7 +440,7 @@ parser = argparse.ArgumentParser() parser.add_argument('--ac', action='store', dest='ac', type=int, help='Allele count value'), - parser.add_argument('--all_isolates', action='store_true', dest='all_isolates', required=False, default=False, help='Create table with all isolates'), + parser.add_argument('--all_isolates', action='store_true', dest='all_isolates', help='Create table with all isolates'), parser.add_argument('--input_excel', action='store', dest='input_excel', required=False, default=None, help='Optional Excel filter file'), parser.add_argument('--input_vcf_dir', action='store', dest='input_vcf_dir', help='Input vcf directory'), parser.add_argument('--min_mq', action='store', dest='min_mq', type=int, help='Minimum map quality value'), diff -r 0fe292b20b9d -r be5875f29ea4 vsnp_get_snps.xml --- a/vsnp_get_snps.xml Thu Jul 29 13:16:03 2021 +0000 +++ b/vsnp_get_snps.xml Tue Aug 03 19:24:50 2021 +0000 @@ -1,13 +1,13 @@ - + macros.xml - openpyxl - pandas - pyvcf - xlrd + + + + &2 echo "The dbkeys associated with the zero coverage VCF files with SNPs found in closely related isolate groups are not unique" && -exit 1 + >&2 echo "The dbkeys associated with the zero coverage VCF files with SNPs found in closely related isolate groups are not unique" && exit 1 #end if #set vcf_identifier = re.sub('[^\s\w\-]', '_', str($i.element_identifier)) ln -s '${i}' '$input_vcf_dir/${vcf_identifier}' && @@ -49,7 +48,7 @@ #if str($input_excel_cond.excel_source_cond.excel_source) == 'cached': #set excel_file = 'No genome specified for input VCF (database) file(s)' #set excel_fields = $__app__.tool_data_tables['vsnp_excel'].get_fields() - ## The value of excel_fields is a nested list that looks like this. + ## The value of excel_fields is a nested list that looks like this. ## [['AF2122', 'Mbovis_define_filter.xlsx', '~/tool-data/vsnp/AF2122/excel/Mbovis_define_filter.xlsx', 'Excel file for AF2122'],...] #for $i in $excel_fields: #if str($i[0]) == $dbkey: @@ -110,7 +109,6 @@ - @@ -122,7 +120,7 @@ - + @@ -154,6 +152,7 @@ + @@ -182,6 +181,7 @@ + @@ -211,6 +211,7 @@ +