Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error with Upper Bound Assembly #232

Open
MartaBenegas opened this issue Dec 23, 2022 · 0 comments
Open

Error with Upper Bound Assembly #232

MartaBenegas opened this issue Dec 23, 2022 · 0 comments

Comments

@MartaBenegas
Copy link

Dear QUAST team,

I'm running QUAST v 5.2.0 inside a docker and testing it with the test data. I haven't managed to generate the Upper Bound Assembly:

root@246f1e3d414e:/data# python3 /quast/quast.py contigs_1.fasta contigs_2.fasta -r reference.fasta.gz --upper-bound-assembly --mp1 reads1.fastq.gz --mp2 reads2.fastq.gz -o output/
WARNING: Python locale settings can't be changed
/quast/quast.py contigs_1.fasta contigs_2.fasta -r reference.fasta.gz --upper-bound-assembly --mp1 reads1.fastq.gz --mp2 reads2.fastq.gz -o output/

Version: 5.2.0, 3d87c606

System information:
  OS: Linux-5.15.0-56-generic-x86_64-with-glibc2.29 (linux_64)
  Python version: 3.8.10
  CPUs number: 12

Started: 2022-12-23 11:01:11

Logging to /data/output/quast.log
NOTICE: Maximum number of threads is set to 3 (use --threads option to set it manually)

CWD: /data
Main parameters: 
  MODE: default, threads: 3, min contig length: 500, min alignment length: 65, min alignment IDY: 95.0, \
  ambiguity: one, min local misassembly length: 200, min extensive misassembly length: 1000

Reference:
  /data/reference.fasta.gz ==> reference

2022-12-23 11:01:12
Generating Upper Bound Assembly...
Compiling BWA (details are in /quast/quast_libs/bwa/make.log and make.err)
Compiling BEDtools (details are in /quast/quast_libs/bedtools/make.log and make.err)
Compiling Minimap2 (details are in /quast/quast_libs/minimap2/make.log and make.err)
Downloading red (file: Red)...
 0.0% of 4468453 bytes
 1.0% of 4468453 bytes
 2.0% of 4468453 bytes
[....]
98.0% of 4468453 bytes
 98.0% of 4468453 bytes
 99.0% of 4468453 bytes
red successfully downloaded!
  Insert size, which will be used for Upper Bound Assembly: auto
  Pre-processing reads...
  Running BWA for reference...
  Done.
  Sorting SAM-file for reference...
  Calculating reads coverage...
  Running repeat masking tool...

ERROR!   Failed to create Upper Bound Assembly, see log for details: /data/output/upper_bound_assembly/upper_bound_assembly.log
[the rest ends correctly]

upper_bound_assembly.log file:

Could not remove: /data/output/upper_bound_assembly/tmp/tmp_red
Invalid State Exception.
The number of states must be positive.
terminate called after throwing an instance of 'exception::InvalidStateException'
  what():  The number of states must be positive.
the default background order: 5.
Using the default threshold: 2.
Calculating GC content ...
Using the default half width: 20 based on the GC content of 52.07

List of final parameters: 
-frm: 2
-gau: 20
-gnm: /data/output/upper_bound_assembly/tmp/tmp_red
-len: 12
-min: 5
-ord: 5
-rpt: /data/output/upper_bound_assembly/tmp
-thr: 2

Deleting pre-existing files under /data/output/upper_bound_assembly/tmp


Stage 1: Building the table ...
Counting k-mers in /data/output/upper_bound_assembly/tmp/tmp_red/reference.fa ...
Processing 0 keys out of 16777216
Processing 10000000 keys out of 16777216
Enrichment view is ready.


Stage 2: Calculating the percentage ...
Calculating the percentage in: /data/output/upper_bound_assembly/tmp/tmp_red/reference.fa ...
The percentage is 100


Stage 3: Training ...
The number of states is: -2147483648

Any help would be appreciated!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant