Skip to content

purge failed to get fasta ids #1

@suvrathaprasad

Description

@suvrathaprasad

I'm in the last step of purge haplotigs. I ran the below command.

purge_haplotigs purge -g m54151_190317_115243_merged_assembly.fasta -c coverage_stats.csv -b 243.sorted.bam

this is the error i get -

30-03-2019 15:18:30 CHECKING bedtools... OK
30-03-2019 15:18:30 CHECKING blastn... OK
30-03-2019 15:18:30 CHECKING lastz... OK
30-03-2019 15:18:30 CHECKING makeblastdb... OK
30-03-2019 15:18:30 CHECKING samtools... OK
30-03-2019 15:18:30 INFO: All required dependencies OK
30-03-2019 15:18:30 CHECKING parallel... OK
30-03-2019 15:18:30 INFO:
Genome fasta: m54151_190317_115243_merged_assembly.fasta
Coverage csv: coverage_stats.csv
Bam file: 243.sorted.bam
Threads: 4
Cutoff, alignment: 70 %
Cutoff, repeat: 250 %
Cutoff, suspect: 40 %
Out prefix: curated
Coverage window len: 9000 bp
Window step dist: 3000 bp
Blastn parameters: -evalue 0.00001 -max_target_seqs 50 -max_hsps 1000 -word_size 24 -culling_limit 10
Lastz paramters: --gfextend --chain --gapped --seed=match14 --allocate:traceback=200.0M

Running using command:
purge_haplotigs purge -g m54151_190317_115243_merged_assembly.fasta -c coverage_stats.csv -b 243.sorted.bam

30-03-2019 15:18:30 INFO:

BEGINNING PIPELINE!

30-03-2019 15:18:30 INFO: reading in genome index file: m54151_190317_115243_merged_assembly.fasta.fai
30-03-2019 15:18:30 INFO: reading in contig coverage stats file: coverage_stats.csv
30-03-2019 15:18:30 INFO: "mincing" genome
30-03-2019 15:18:30 ERROR: failed to get id from fasta file

PIPELINE FAILURE

any idea how to go about this?
thanks in advance.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions