Guest User

Untitled

a guest
Nov 15th, 2018
112
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.48 KB | None | 0 0
  1. EXITING: FATAL INPUT ERROR: duplicate parameter "outSAMtype" in input "Command-Line"
  2. SOLUTION: keep only one definition of input parameters in each input source
  3.  
  4. Nov 15 09:46:07 ...... FATAL ERROR, exiting
  5. logs/star/se/UY2_S7.log (END)
  6.  
  7. __author__ = "Johannes Köster"
  8. __copyright__ = "Copyright 2016, Johannes Köster"
  9. __email__ = "koester@jimmy.harvard.edu"
  10. __license__ = "MIT"
  11.  
  12. import os
  13. from snakemake.shell import shell
  14.  
  15. extra = snakemake.params.get("extra", "")
  16. log = snakemake.log_fmt_shell(stdout=True, stderr=True)
  17.  
  18. fq1 = snakemake.input.get("fq1")
  19. assert fq1 is not None, "input-> fq1 is a required input parameter"
  20. fq1 = [snakemake.input.fq1] if isinstance(snakemake.input.fq1, str) else snakemake.input.fq1
  21. fq2 = snakemake.input.get("fq2")
  22. if fq2:
  23. fq2 = [snakemake.input.fq2] if isinstance(snakemake.input.fq2, str) else snakemake.input.fq2
  24. assert len(fq1) == len(fq2), "input-> equal number of files required for fq1 and fq2"
  25. input_str_fq1 = ",".join(fq1)
  26. input_str_fq2 = ",".join(fq2) if fq2 is not None else ""
  27. input_str = " ".join([input_str_fq1, input_str_fq2])
  28.  
  29. if fq1[0].endswith(".gz"):
  30. readcmd = "--readFilesCommand zcat"
  31. else:
  32. readcmd = ""
  33.  
  34. outprefix = os.path.dirname(snakemake.output[0]) + "/"
  35.  
  36. shell(
  37. "STAR "
  38. "{extra} "
  39. "--runThreadN {snakemake.threads} "
  40. "--genomeDir {snakemake.params.index} "
  41. "--readFilesIn {input_str} "
  42. "{readcmd} "
  43. "--outSAMtype BAM Unsorted "
  44. "--outFileNamePrefix {outprefix} "
  45. "--outStd Log "
  46. "{log}")
Add Comment
Please, Sign In to add comment