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

LoadSQM issue #920

Open
tboonf opened this issue Jan 13, 2025 · 6 comments
Open

LoadSQM issue #920

tboonf opened this issue Jan 13, 2025 · 6 comments

Comments

@tboonf
Copy link

tboonf commented Jan 13, 2025

SqueezeMeta completed successfully with all the necessary outputs. However i have issue with loading it to SQMTools (both squeezemeta and sqmtools are the latest versions):

ABC = loadSQM('./Project', engine = 'data.table')
Generating tabular outputs for project in ./Project
Traceback (most recent call last):
File "/home/nova/miniconda3/envs/SqueezeMeta/SqueezeMeta/utils/sqm2tables.py", line 223, in
main(parse_args())
File "/home/nova/miniconda3/envs/SqueezeMeta/SqueezeMeta/utils/sqm2tables.py", line 70, in main
sampleNames, orfs, kegg, cog, pfam, custom, noCDSorfs, noCDScontigs = parse_orf_table(perlVars['$mergedfile'], total_reads, total_bases, nokegg, nocog, nopfam,
File "/home/nova/miniconda3/envs/SqueezeMeta/SqueezeMeta/utils/../lib/utils.py", line 199, in parse_orf_table
abundances = array([int(line[idx[sample]]) for sample in samples])
File "/home/nova/miniconda3/envs/SqueezeMeta/SqueezeMeta/utils/../lib/utils.py", line 199, in
abundances = array([int(line[idx[sample]]) for sample in samples])
IndexError: list index out of range

@fpusan
Copy link
Collaborator

fpusan commented Jan 13, 2025

What version of SqueezeMeta did you use to generate this project?
Any chance you can share the full directory with us using something like dropbox, wetransfer, etc?

@tboonf
Copy link
Author

tboonf commented Jan 14, 2025

Thanks for fast responding..is it possible to send the link to your private email due to data confidentiality?

@fpusan
Copy link
Collaborator

fpusan commented Jan 14, 2025

Of course [email protected]

@fpusan
Copy link
Collaborator

fpusan commented Jan 16, 2025

Did you send it in the end? it doesn't seem so but maybe it went to my spam folder?

@fpusan
Copy link
Collaborator

fpusan commented Jan 17, 2025

I see what's happening now. The last line on your results/13.*.orftable file is truncated. I wonder how the data got corrupted.
Maybe your process was killed in the middle of that step, and you restarted the pipeline?

@tboonf
Copy link
Author

tboonf commented Jan 22, 2025 via email

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

2 participants