Changeset 5773
- Timestamp:
- Dec 3, 2019, 1:40:58 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
extensions/net.sf.basedb.reggie/trunk/src/net/sf/basedb/reggie/grid/VariantStatisticsJobCreator.java
r5772 r5773 122 122 script.cmd("cd ${TMPDIR}"); 123 123 script.cmd("cp ${ScriptDir}/mut_stats.py ."); 124 script.cmd("mkdir results"); 125 script.cmd("mkdir tmp"); 124 126 script.newLine(); 125 127 … … 128 130 statCmd += " ${WD}/progress"; 129 131 statCmd += " ${WD}/mut_stats.log"; 130 statCmd += " > mut_stats.vcf";132 statCmd += " > tmp/mut_stats.vcf"; 131 133 script.cmd(statCmd); 132 134 script.progress(95, "Sorting and indexing..."); 133 script.cmd(bedtools_path+" sort -header -i mut_stats.vcf | bgzip -c > mut_stats.vcf.gz"); 134 script.cmd("tabix mut_stats.vcf.gz"); 135 script.cmd(bedtools_path+" sort -header -i tmp/mut_stats.vcf | bgzip -c > results/mut_stats.vcf.gz"); 136 script.cmd("tabix results/mut_stats.vcf.gz"); 137 script.cmd("cp results/* ${WD}"); 135 138 136 139 script.progress(99, "Cleaning up temporary folders"); … … 215 218 216 219 return "Parsed " + stat.numVcf + " VCF files for " + stat.numPatients + " patients." + 217 " Found "+Reggie.formatCount(stat.numVariants) + " variants."; 220 " Found "+Reggie.formatCount(stat.numVariants) + " variants." + 221 " Result files can be found at: " + session.getHost().getWorkFolder(jobName); 218 222 } 219 223 }
Note: See TracChangeset
for help on using the changeset viewer.