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

            • kinpict.com
              Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by kinpict.com
              The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
              Today, 06:37 PM
            • SEQadmin2
              16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by SEQadmin2





              NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
              10-30-2026, 06:46 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 11-12-2026, 07:59 AM
            0 responses
            28 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-05-2026, 09:25 AM
            0 responses
            35 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-03-2026, 11:54 AM
            2 responses
            56 views
            0 reactions
            Last Post kinpict.com  
            Started by SEQadmin2, 10-28-2026, 09:57 AM
            0 responses
            175 views
            0 reactions
            Last Post SEQadmin2  
            Working...