Problem of running samtools in python using subprocess.run()

Question

I'm using subprocesss.run() to run samtools command in python. The code is as below:

result = subprocess.run(['samtools', 'faidx', 'hg38.fa.gz', 'chr1:169699712-169699719'], check=True, stdout = subprocess.PIPE)

and I got the following problem:

CalledProcessError: Command '['samtools', 'faidx', 'hg38.fa.gz', 'chr1:169699712-169699719']' died with <Signals.SIGABRT: 6>.

The samtools command ran successfully in terminal, but failed in subprocess.run.

Does anyone know the cause of the bug? Thank you so much.

Try printing the `results.stderr` and `results.stdout` attributes of your result to see what error samtools threw exactly. — Pallie, Feb 09 '21 at 13:16

score 0 · Answer 1 · answered Mar 25 '21 at 15:19

0

Maybe try with this code, I use a similar one to run commands with bwa, samtools and others without issues so far:

import subprocess

def run(cmd) :
   proc = subprocess.Popen(cmd, shell=True, stdout=subprocess.PIPE)
   proc.communicate()

cmd = "samtools faidx /path/to/file.fa.gz" # change for desired command line here
run(cmd)

answered Mar 25 '21 at 15:19

AnHo

123
6

Thank you so much! The problem is solved. – Aragaki Mar 26 '21 at 03:11

score 0 · Answer 2 · answered Apr 16 '22 at 20:09

Try the command below:

P1 = subprocess.run('samtools faidx hg38.fa.gz chr1:169699712-169699719', capture_output=True, text=True, shell=True)
print(P1.returncode) # will show 0 if execution is successful
print(P1.stdout) # will contain the output in str format

Problem of running samtools in python using subprocess.run()

2 Answers2