Question: What is "Error in value[[3L]](cond) : 'asBam' failed to open SAM/BAM file"
0
gravatar for rj.innet
6 months ago by
rj.innet0
rj.innet0 wrote:

I am encountering the below error message while creating bampaths in systemPipeR workflow running in the RStudio.

   > bampaths <- runCommandline(args=args)
  Error in value[[3L]](cond) : 'asBam' failed to open SAM/BAM file

     file: '/home/rna/sample/results/SRR446028_1.fastq.gz.hisat.sam'
      SAM file: '/home/rna/sample/results/SRR446028_1.fastq.gz.hisat.sam'
    In addition: Warning message:
    In system2(command, args = commandargs, stdout = TRUE, stderr = TRUE) :
      running command ''/home/rna/tool/hisat2-2.1.0' -p 4 -k 1 --min-intronlen 30 --max-intronlen 3000 -S /home/rna/sample/results/SRR446028_1.fastq.gz.hisat.sam /home/rna/sample/data/tair10.fasta -U /home/rna/sample//home/rna/sample/data/SRR446028_1.fastq.gz /home/rna/sample/ 2>&1' had status 126

What is this [[3L]] (cond).... error?? Any suggestions to solve this error will be very helpful. Thanks!!

sam/bam rna-seq systempiper • 438 views
ADD COMMENTlink written 6 months ago by rj.innet0
1

As people already suggested in your previous post about bowtie2, do yourself a favor and run bash commands via the command line, rather than from inside R. It only adds unnecessary complication.

ADD REPLYlink written 6 months ago by ATpoint15k

exit status: 126 for “command not executable". See if you can execute in shell:

/home/rna/tool/hisat2-2.1.0 -p 4 -k 1 --min-intronlen 30 --max-intronlen 3000 -S /home/rna/sample/results/SRR446028_1.fastq.gz.hisat.sam /home/rna/sample/data/tair10.fasta -U /home/rna/sample//home/rna/sample/data/SRR446028_1.fastq.gz /home/rna/sample/ 2>&1

Check if hisat is executable.

ADD REPLYlink modified 6 months ago • written 6 months ago by cpad011211k

I have ran the following command and it gave me the result shown below.

hisat2 -p 4 -k 1 --min-intronlen 30 --max-intronlen 3000 -S /home/rna/sample/results/SRR446027_1.fastq.gz.hisat.sam /home/rna/sample/data/tair10.fasta -U /home/rna/sample/data/SRR446027_1.fastq.gz

96459 reads; of these: 96459 (100.00%) were unpaired; of these: 3968 (4.11%) aligned 0 times 92491 (95.89%) aligned exactly 1 time 0 (0.00%) aligned >1 times 95.89% overall alignment rate

I guess hisat2 is working fine but 3L error is comming when I run

bampaths <- runCommandline(args=args)
ADD REPLYlink written 6 months ago by rj.innet0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1234 users visited in the last hour