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

summarizeFasta error with combination sources #827

Closed
lydiayliu opened this issue Nov 14, 2023 · 1 comment
Closed

summarizeFasta error with combination sources #827

lydiayliu opened this issue Nov 14, 2023 · 1 comment
Assignees
Labels
bug Something isn't working

Comments

@lydiayliu
Copy link
Collaborator

as mentioned in https://github.com/uclahs-cds/project-MissingPeptides-Method/issues/254#issuecomment-1811354683

summarizeFasta doesn't work with combination sources it seems? I'm getting this error

Traceback (most recent call last):
  File "/usr/local/bin/moPepGen", line 8, in <module>
    sys.exit(main())
  File "/usr/local/lib/python3.8/site-packages/moPepGen/cli/__main__.py", line 92, in main
    args.func(args)
  File "/usr/local/lib/python3.8/site-packages/moPepGen/cli/summarize_fasta.py", line 221, in summarize_fasta
    summarizer.write_summary_table(handle)
  File "/usr/local/lib/python3.8/site-packages/moPepGen/aa/PeptidePoolSummarizer.py", line 313, in write_summary_table
    if self.contains_exclusive_sources(comb):
  File "/usr/local/lib/python3.8/site-packages/moPepGen/aa/PeptidePoolSummarizer.py", line 275, in contains_exclusive_sources
    parsers = self.get_parsers_from_source(source)
  File "/usr/local/lib/python3.8/site-packages/moPepGen/aa/PeptidePoolSummarizer.py", line 206, in get_parsers_from_source
    parser = self.source_parser_map[s]
KeyError: frozenset({'Variant', 'Alt'})

The fastas are in /hot/project/method/AlgorithmDevelopment/ALGO-000074-moPepGen/CPCGENE/processed/noncanonical-database/call-nonCanonicalPeptide/2023-11-14_split/pipeline-meta-call-NonCanonicalPeptide-0.0.1/CPCG0100/call-NonCanonicalPeptide-1.0.0/CPCG0100/moPepGen-dev/output/split/

My command is

docker run -u $(id -u):$(id -g) --rm -v /hot/:/hot/ ghcr.io/uclahs-cds/mopepgen:dev \
moPepGen summarizeFasta \
    --variant-peptides ${d} \
    --gvf ${c} \
    --index-dir /hot/project/method/AlgorithmDevelopment/ALGO-000074-moPepGen/ref/GRCh38-EBI-GENCODE34/index/ \
    --group-source Alt:SECT,CodonReassign Variant:gSNP,gIndel,sSNV,sIndel,Fusion,altSplice,RNAEdit \
    --order-source Variant,Alt,Variant-Alt,Noncoding,Noncoding-Variant,Noncoding-Alt,Noncoding-Variant-Alt,circRNA,circRNA-Variant,circRNA-Alt,circRNA-Variant-Alt,Noncoding-circRNA,Noncoding-circRNA-Variant,Noncoding-circRNA-Alt,Noncoding-circRNA-Variant-Alt \
    --output-path ${d}.group.txt \
@lydiayliu lydiayliu added the bug Something isn't working label Nov 14, 2023
@lydiayliu
Copy link
Collaborator Author

Should be address in #828

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants