assembly_accuracy icon indicating copy to clipboard operation
assembly_accuracy copied to clipboard

Fastmer.py error

Open Ashma45 opened this issue 5 years ago • 5 comments

Hi, I am using assembly_accuracy first time. I am getting the following error. Can you help me fix the issue?

Thank you

(base) [achakraw@vacc-user1 assembly_accuracy]$ ~/WGS/test/raw_reads/prokka_reads/assembly_accuracy/fastmer.py --reference GCF_000013425.1_ASM1342v1_genomic.fna --assembly 105_r2_prokka/*fsa --min-mapping-quality 10

Traceback (most recent call last): File "/users/a/c/achakraw/WGS/test/raw_reads/prokka_reads/assembly_accuracy/fastmer.py", line 304, in subprocess.check_call(mm2_cmd, stdout=devnull, stderr=devnull, shell=True) File "/users/a/c/achakraw/miniconda3/lib/python3.6/subprocess.py", line 311, in check_call raise CalledProcessError(retcode, cmd) subprocess.CalledProcessError: Command 'minimap2 -Y -a -x asm5 GCF_000013425.1_ASM1342v1_genomic.fna 105_r2_prokka/*fsa | samtools sort -T 105_r2_prokka/*fsa.assembly_analysis.sorted.bam.tmp -o 105_r2_prokka/*fsa.assembly_analysis.sorted.bam -' returned non-zero exit status 127.

Ashma45 avatar Sep 03 '20 00:09 Ashma45

Hi,

This repository is largely unmaintained so I apologize for the poor documentation and lack of error handling.

This error is likely because minimap2 is not found on your PATH.

Jared

On Sep 2, 2020, at 8:28 PM, Ashma Chakrawarti [email protected] wrote:

 Hi, I am using assembly_accuracy first time. I am getting the following error. Can you help me fix the issue?

Thank you

(base) [achakraw@vacc-user1 assembly_accuracy]$ ~/WGS/test/raw_reads/prokka_reads/assembly_accuracy/fastmer.py --reference GCF_000013425.1_ASM1342v1_genomic.fna --assembly 105_r2_prokka/*fsa --min-mapping-quality 10 Traceback (most recent call last): File "/users/a/c/achakraw/WGS/test/raw_reads/prokka_reads/assembly_accuracy/fastmer.py", line 304, in subprocess.check_call(mm2_cmd, stdout=devnull, stderr=devnull, shell=True) File "/users/a/c/achakraw/miniconda3/lib/python3.6/subprocess.py", line 311, in check_call raise CalledProcessError(retcode, cmd) subprocess.CalledProcessError: Command 'minimap2 -Y -a -x asm5 GCF_000013425.1_ASM1342v1_genomic.fna 105_r2_prokka/*fsa | samtools sort -T 105_r2_prokka/*fsa.assembly_analysis.sorted.bam.tmp -o 105_r2_prokka/*fsa.assembly_analysis.sorted.bam -' returned non-zero exit status 127.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or unsubscribe.

jts avatar Sep 03 '20 00:09 jts

Thank you, Jared!

I too doubted that minimap2 was not present on my PATH. But it is present when I check using minimap -h I see the documentation mentions pyvcf and pysam. I installed both using conda but when I call pyvcf -h it returns as command not found.

Ashma45 avatar Sep 03 '20 13:09 Ashma45

Is samtools available on your PATH as well?

jts avatar Sep 05 '20 14:09 jts

Yes, the samtools is available on PATH. This time the error is as following

Traceback (most recent call last): File "./fastmer.py", line 6, in import vcf ModuleNotFoundError: No module named 'vcf'

From: Jared Simpson [email protected] Sent: Saturday, September 5, 2020 10:30 AM To: jts/assembly_accuracy [email protected] Cc: Ashma Chakrawarti [email protected]; Author [email protected] Subject: Re: [jts/assembly_accuracy] Fastmer.py error (#1)

Is samtools available on your PATH as well?

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/jts/assembly_accuracy/issues/1#issuecomment-687619072, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AOGDQ5LTC7IBSWMXP2Y4PFTSEJDQDANCNFSM4QUAO34Q.

Ashma45 avatar Sep 06 '20 14:09 Ashma45

Both of those are python packages - can you try to install them using pip?

Jared

On Thu, Sep 10, 2020 at 1:20 PM Ashma Chakrawarti [email protected] wrote:

Thank you, Jared!

I too doubted that minimap2 was not present on my PATH. But it is present when I check using minimap -h I see the documentation mentions pyvcf and pysam. I installed both using conda but when I call pyvcf -h it returns as command not found.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/jts/assembly_accuracy/issues/1#issuecomment-686476132, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAC7DH4V7SSWBJJBFDLYMWTSFEDFRANCNFSM4QUAO34Q .

jts avatar Sep 10 '20 19:09 jts