command_line = ("{7} {0} -Xmx{1} -jar {2} -T Pileup -R {3} -I {4} -L {5} -o {6} " +
"-verbose -rf DuplicateRead --filter_reads_with_N_cigar " +
"--filter_mismatching_base_and_quals").format(JAVA_TEMP, opts.xmx_java, GATK, REFERENCE, opts.bam, MARKER_FILE, opts.outfile, opts.java)

os.system(command_line)

相关文章:

  • 2022-12-23
  • 2021-11-20
  • 2021-07-22
  • 2022-02-11
  • 2021-11-18
  • 2021-05-21
  • 2022-01-04
  • 2022-12-23
猜你喜欢
  • 2021-09-20
  • 2022-12-23
  • 2021-12-17
  • 2022-12-23
  • 2022-12-23
相关资源
相似解决方案