Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • libiyagirl
    Member
    • Jan 2011
    • 27

    #1

    can I automate running dindel?

    Hello,

    I am writing a script for automatically running dindel commands. But it never works.

    For example, whenever I use:
    "input_file=/path/input_file.bam
    dindel --bamFile $input_file"

    it can go to the path and find the file, but can not open or read it. I'm pretty sure the method is correct, because it works fine on samtools.

    Any idea regarding this?

    Thanks a lot,
  • Michael.James.Clark
    Senior Member
    • Apr 2009
    • 207

    #2
    Might be the permissions of the file. Try:
    Code:
    chmod 755 /path/input_file.bam
    Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
    Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
    Projects: U87MG whole genome sequence [Website] [Paper]

    Comment

    • libiyagirl
      Member
      • Jan 2011
      • 27

      #3
      Thank you Michael, it seems like chmod 755 can not change the permission of input.bam
      Any solution?

      Comment

      • Michael.James.Clark
        Senior Member
        • Apr 2009
        • 207

        #4
        Try putting the input file path in quotes?

        I.e.

        Code:
        input_file="/path/to/input_file.bam";
        dindel --bamFile $input_file
        Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
        Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
        Projects: U87MG whole genome sequence [Website] [Paper]

        Comment

        • mgogol
          Senior Member
          • Mar 2008
          • 197

          #5
          Can you view the bam file?

          samtools view input_file.bam | head

          If you do ls -l on the bam file, what are the permissions that are shown?

          What error message are you getting?

          Comment

          • libiyagirl
            Member
            • Jan 2011
            • 27

            #6
            Thank you all, I used
            input="$PATH/sorted.bam";
            and it's working now!

            Thanks a lot~!

            Comment

            Latest Articles

            Collapse

            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
              ...
              10-14-2026, 09:55 AM
            • SEQadmin2
              New Genomics Technologies Take Aim at Long-Standing Limits
              by SEQadmin2


              Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

              We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
              ...
              09-28-2026, 10:25 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, Yesterday, 11:13 AM
            0 responses
            16 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, Yesterday, 10:48 AM
            0 responses
            15 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-15-2026, 10:52 AM
            0 responses
            35 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-14-2026, 05:39 AM
            0 responses
            42 views
            0 reactions
            Last Post SEQadmin2  
            Working...