Mercurial > repos > jjohnson > shear

changeset 8:8fac6c7d32fd

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Copy reference fasta for shear sv to avoid GATK issue with symlinks
author Jim Johnson <jj@umn.edu>
date Tue, 22 Oct 2013 09:44:32 -0500
parents 8180f6b0249b
children 03071ff2e5c9
files shear_wrapper.py
diffstat 1 files changed, 2 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/shear_wrapper.py	Tue Oct 22 08:20:52 2013 -0500
+++ b/shear_wrapper.py	Tue Oct 22 09:44:32 2013 -0500
@@ -47,7 +47,7 @@
     args.append( command )
     # Check required params for command
     buffsize = 1048576 
-    copy_ref = False
+    copy_ref = True
     prefix = options.prefix if options.prefix and len(options.prefix) > 0 else 'ref'
     if command in ['sv']:
         args.append( '-p' )
@@ -57,6 +57,7 @@
         if options.region:
             args.append('-r')
             args.append(options.region)
+        copy_ref = True # GATK has issue with symlinks for .fa and .fai 
     if options.svidx_dir and command in ['sv']:
         if not os.path.isdir(options.svidx_dir):
             os.makedirs(options.svidx_dir)