Comments (8)
Hi!
Can you share the input genomes and the command you used so I can try to reproduce the error?
from superpang.
Species: Mycobacterium tuberculosis
Taxid: 1773
Genome:316 complete genomes in GenBank (https://ftp.ncbi.nlm.nih.gov/genomes/genbank/assembly_summary_genbank.txt)
Command: SuperPang.py --fasta genome/*.fa --output-dir 1773/pangenome --force-overwrite -t 20 --assume-complete -b 0.95 -i 0.95 -k 301
from superpang.
You can also test E.coli with all complete genomes in GenBank, get the same error.
from superpang.
Can reproduce with the Mycobacterium tuberculosis dataset
from superpang.
Full errors:
multiprocessing.pool.RemoteTraceback:
"""
Traceback (most recent call last):
File ".../anaconda3/envs/SuperPang/lib/python3.8/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File ".../anaconda3/envs/SuperPang/lib/python3.8/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File ".../anaconda3/envs/SuperPang/lib/python3.8/site-packages/SuperPang/lib/Assembler.py", line 846, in reconstruct_sequence
if k[:-1] == kmers[-1][1:]:
IndexError: list index out of range
"""
The above exception was the direct cause of the following exception:
Traceback (most recent call last):
File ".../anaconda3/envs/SuperPang/bin/SuperPang.py", line 291, in
main(parse_args())
File ".../anaconda3/envs/SuperPang/bin/SuperPang.py", line 144, in main
contigs = Assembler(input_minimap2, args.ksize, args.threads).run(args.minlen, args.mincov, args.bubble_identity_threshold, args.genome_assignment_threshold, args.threads)
File ".../anaconda3/envs/SuperPang/lib/python3.8/site-packages/SuperPang/lib/Assembler.py", line 259, in run
NBP2seq = dict(zip(NBPs, self.multimap(self.reconstruct_sequence, threads, NBPs, self.vertex2hash, self.compressor)))
File ".../anaconda3/envs/SuperPang/lib/python3.8/site-packages/SuperPang/lib/Assembler.py", line 55, in multimap
res = pool.map(fun, range(len(iterable)))
File ".../anaconda3/envs/SuperPang/lib/python3.8/multiprocessing/pool.py", line 364, in map
return self._map_async(func, iterable, mapstar, chunksize).get()
File ".../anaconda3/envs/SuperPang/lib/python3.8/multiprocessing/pool.py", line 771, in get
raise self._value
IndexError: list index out of range
from superpang.
Thanks! I'm running it right now in my computers, adding some traces to see what could be happening.
Also, wild guess here, but does this also happen if you call it like
SuperPang.py --fasta genome/*.fa --output-dir 1773/pangenome --force-overwrite -t 20 --assume-complete -b 0.95 -i 0.95 -k 301 -n 20
?
from superpang.
Nevermind, I now think that the problem is due to the presence of ambiguous bases in the input fasta files
(e.g. "M" to represent either "A" or "C"). This happens for these genomes and would definitely break my code in a few places. Will think of a fix
from superpang.
Hi again! This issue is hopefully fixed in the new version 0.9.3
from superpang.
Related Issues (15)
- IndexError: list index out of range HOT 2
- Read checkm file at the begining even if not yet needed
- multiple file paths for `--fasta` HOT 1
- Error after identifying connected components HOT 15
- the pangeome size is larger than all size of input genomes HOT 7
- How to improve the mOTUpan core-genome estimation for viruses species? HOT 4
- mOTUpan HOT 4
- errors when running the test-SuperPang.py HOT 4
- There was an error running homogenize.py HOT 12
- There was an error running homogenize.py HOT 13
- ValueError: invalid vertex: 3; condense-edges.py error. HOT 3
- assert len(kmers) == 2 AssertionError HOT 5
- input genome position information HOT 2
- NBP2origins.tsv and graph.NBP2origins.csv contain corrupted SourceContig names HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from superpang.