diff shear_wrapper.py @ 8:8fac6c7d32fd

Copy reference fasta for shear sv to avoid GATK issue with symlinks

--- 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)