Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

failed with own error model

See original GitHub issue

–Hi, i have tried to create own error model with:

bowtie2-build --threads 16 -f chr8.fasta bowtie2_index/chr8.genome
bowtie2 -p 16 -x bowtie2_index/chr8.genome -1 chr8_paired1.fq.gz -2 chr8_paired2.fq.gz | samtools view -@16 -bS | samtools sort -@16 -o genomes.bam
samtools index genomes.bam
bin/iss model -b genomes.bam -o genomes

(chr8_paired1.fq.gz and chr8_paired2.fq.gz are paired-end reads 50bp length)

then program failed with: bin/iss generate --cpus 8 -n 180M --genomes chr8.fasta --model genomes.npz -z --output reads

INFO:iss.app:Starting iss generate
INFO:iss.app:Using kde ErrorModel
INFO:iss.util:Stitching input files together
INFO:iss.app:Using lognormal abundance distribution
INFO:iss.app:Using 8 cpus for read generation
INFO:iss.app:Generating 180000000 reads
INFO:iss.app:Generating reads for record: chr8
joblib.externals.loky.process_executor._RemoteTraceback:
"""
Traceback (most recent call last):
  File ".local/lib/python3.7/site-packages/joblib/externals/loky/process_executor.py", line 418, in _process_worker
    r = call_item()
  File ".local/lib/python3.7/site-packages/joblib/externals/loky/process_executor.py", line 272, in __call__
    return self.fn(*self.args, **self.kwargs)
  File ".local/lib/python3.7/site-packages/joblib/_parallel_backends.py", line 608, in __call__
    return self.func(*args, **kwargs)
  File ".local/lib/python3.7/site-packages/joblib/parallel.py", line 256, in __call__
    for func, args, kwargs in self.items]
  File ".local/lib/python3.7/site-packages/joblib/parallel.py", line 256, in <listcomp>
    for func, args, kwargs in self.items]
  File ".local/lib/python3.7/site-packages/iss/generator.py", line 65, in reads
    forward, reverse = simulate_read(record, ErrorModel, i, cpu_number)
  File ".local/lib/python3.7/site-packages/iss/generator.py", line 140, in simulate_read
    forward = ErrorModel.introduce_error_scores(forward, 'forward')
  File ".local/lib/python3.7/site-packages/iss/error_models/__init__.py", line 66, in introduce_error_scores
    self.quality_forward, 'forward')
  File ".local/lib/python3.7/site-packages/Bio/SeqRecord.py", line 88, in __setitem__
    "strings) of length {0}.".format(self._length)
TypeError: We only allow python sequences (lists, tuples or strings) of length 50.
"""

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
  File "bin/iss", line 11, in <module>
    load_entry_point('InSilicoSeq==1.4.5', 'console_scripts', 'iss')()
  File ".local/lib/python3.7/site-packages/iss/app.py", line 547, in main
    args.func(args)
  File ".local/lib/python3.7/site-packages/iss/app.py", line 260, in generate_reads
    args.gc_bias, mode) for i in range(cpus))
  File ".local/lib/python3.7/site-packages/joblib/parallel.py", line 1017, in __call__
    self.retrieve()
  File ".local/lib/python3.7/site-packages/joblib/parallel.py", line 909, in retrieve
    self._output.extend(job.get(timeout=self.timeout))
  File ".local/lib/python3.7/site-packages/joblib/_parallel_backends.py", line 562, in wrap_future_result
    return future.result(timeout=timeout)
  File "/tools/python/3.7.4/lib/python3.7/concurrent/futures/_base.py", line 435, in result
    return self.__get_result()
  File "/tools/python/3.7.4/lib/python3.7/concurrent/futures/_base.py", line 384, in __get_result
    raise self._exception
TypeError: We only allow python sequences (lists, tuples or strings) of length 50.

what’s wrong ?