peptide-shaker icon indicating copy to clipboard operation
peptide-shaker copied to clipboard

Spectrum not found when parsing sage results

Open samgregoire opened this issue 3 months ago • 3 comments

Hello. I'm trying to use PeptideShaker to process identification files obtained by SearchGUI using several search engines (Comet, OMSSA, Tide, X! Tandem and Sage). PeptideShaker encounters the following error when trying to process the Sage output file.

Report:

------------------------------------------------------------------

PeptideShaker 3.0.11 Report File

Originally saved by: sgregoire @ Precision-5820-Tower-X-Series

on: 22 Oct 2024, 20:17

------------------------------------------------------------------

Tue Oct 22 20:17:23 CEST 2024 Unzipping searchgui_out.zip. Tue Oct 22 20:17:29 CEST 2024 Import process for scPTM1

Tue Oct 22 20:17:29 CEST 2024 Importing sequences from uniprotkb_proteome_UP000005640_AND_revi_2024_08_26.fasta. Tue Oct 22 20:17:30 CEST 2024 Importing gene mappings. Tue Oct 22 20:17:31 CEST 2024 Establishing local database connection. Tue Oct 22 20:17:31 CEST 2024 Reading identification files. Tue Oct 22 20:17:31 CEST 2024 Parsing 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.comet.pep.xml.gz. Tue Oct 22 20:17:31 CEST 2024 No PSM found in 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.comet.pep.xml.gz. Tue Oct 22 20:17:31 CEST 2024 Parsing 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.omx.gz. Tue Oct 22 20:17:51 CEST 2024 Checking spectra for 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.omx.gz. Tue Oct 22 20:17:51 CEST 2024 Importing PSMs from 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.omx.gz Tue Oct 22 20:17:55 CEST 2024 744 identified spectra (2.9%) did not present a valid peptide. Tue Oct 22 20:17:55 CEST 2024 3745 of the best scoring peptides were excluded by the import filters: Tue Oct 22 20:17:55 CEST 2024 - 67.6% peptide length less than 7 or greater than 30. Tue Oct 22 20:17:55 CEST 2024 - 32.4% peptide presenting high mass or isotopic deviation. Tue Oct 22 20:17:55 CEST 2024 Parsing 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.sage.tsv.gz. Tue Oct 22 20:17:55 CEST 2024 Checking spectra for 240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.sage.tsv.gz. Tue Oct 22 20:17:55 CEST 2024 Spectrum with title 'frame=9411 scan=411' in file named '240726_MClassChip_24_001721_CBio11_LSM_1_1_4350' required to parse '240726_MClassChip_24_001721_CBio11_LSM_1_1_4350.sage.tsv.gz' not found.

Tue Oct 22 20:17:55 CEST 2024 Importing Data Canceled!

Please, ignore the fact that no PSM was found in the Comet file, I encoutered an error during the search.

I tried to run PeptideShaker without the Sage identification file and everything seems to work fine. I tried to process only the Sage identification file alone and ran into the same error. After checking, I can guarantee that the mzML file I'm using contains a 'frame=9411 scan=411' scan. Removing this specific PSM from the identification files only changes the error message to another scan title.

Please, let me know if I can provide any more useful information.

Thanks, Sam

EDIT: perhaps the problem encountered during Comet processing is actually related to this problem? I get the following error message:

Load spectra:free(): invalid pointer

I don't really understand what this means, but it seems to be related to spectra loading, which could be related to the issue I'm describing above.

samgregoire avatar Oct 29 '24 15:10 samgregoire