# HG changeset patch # User Daniel Blankenberg # Date 1368710879 14400 # Node ID 4a67169be0eedbaed7b9309783cac8370eb861b2 # Parent 1c23da12b7681c0d6c12c8a6f7fea62a594c8cf5 Update to allow out of bounds positions to be not be fatal errors. diff -r 1c23da12b768 -r 4a67169be0ee dependency_configs/tool_dependencies.xml --- a/dependency_configs/tool_dependencies.xml Wed May 15 09:37:01 2013 -0400 +++ b/dependency_configs/tool_dependencies.xml Thu May 16 09:27:59 2013 -0400 @@ -3,9 +3,9 @@ - + - + diff -r 1c23da12b768 -r 4a67169be0ee tools/naive_variant_detector.py --- a/tools/naive_variant_detector.py Wed May 15 09:37:01 2013 -0400 +++ b/tools/naive_variant_detector.py Thu May 16 09:27:59 2013 -0400 @@ -19,6 +19,7 @@ parser.add_option( '-q', '--min_base_quality', dest='min_base_quality', action='store', type="int", default=None, help='Minimum base quality.' ) parser.add_option( '-m', '--min_mapping_quality', dest='min_mapping_quality', action='store', type="int", default=None, help='Minimum mapping.' ) parser.add_option( '-t', '--coverage_dtype', dest='coverage_dtype', action='store', type="string", default='uint8', help='dtype to use for coverage array' ) + parser.add_option( '--allow_out_of_bounds_positions', dest='allow_out_of_bounds_positions', action='store_true', default = False, help='Allows out of bounds positions to not throw fatal errors' ) parser.add_option( '--region', dest='region', action='append', type="string", default=[], help='region' ) (options, args) = parser.parse_args() @@ -55,7 +56,7 @@ coverage = VCFReadGroupGenotyper( map( lambda x: Reader( *x ), bam_files ), options.reference_genome_filename, dtype=options.coverage_dtype, min_support_depth=options.min_support_depth, min_base_quality=options.min_base_quality, min_mapping_quality=options.min_mapping_quality, - restrict_regions=regions, use_strand=options.use_strand ) + restrict_regions=regions, use_strand=options.use_strand, allow_out_of_bounds_positions=options.allow_out_of_bounds_positions ) for line in coverage.iter_vcf( ploidy=options.ploidy, variants_only=options.variants_only ): out.write( "%s\n" % line ) out.close() diff -r 1c23da12b768 -r 4a67169be0ee tools/naive_variant_detector.xml --- a/tools/naive_variant_detector.xml Wed May 15 09:37:01 2013 -0400 +++ b/tools/naive_variant_detector.xml Thu May 16 09:27:59 2013 -0400 @@ -45,6 +45,8 @@ --coverage_dtype "${coverage_dtype}" + --allow_out_of_bounds_positions +