Add support for reading Sage PSM files; various minor fixes

[ArthurDeclercq]

Added

• Add reader for Sage PSM files.

Changed

• psm: The default values of PSM.provenance_data, PSM.metadata and PSM.rescoring_features are now dict() instead of None.
• io.mzid.MzidReader: Attempt to parse retention time or scan start time cvParams from both SpectrumIdentificationResult as SpectrumIdentificationItem levels. Note that according to the mzIdentML specification document (v1.1.1) neither cvParams are expected to be present at either level.
• io.mzid.MzidReader: Prefer spectrum title cvParam over spectrumID attribute for PSM.spectrum_id as these titles always match to the peak list files. In this case, spectrumID is saved in metadata["mzid_spectrum_id"]. Fall back to spectrumID if spectrum title is absent.
• io.mzid.MzidWriter: PSM.retention_time is now written as cvParam retention time instead of scan start time, and to the SpectrumIdentificationItem level instead of the SpectrumIdentificationResult level, as theoretically in psm_utils, multiple PSMs for the same spectrum can have different values for retention_time.
• io.mzid.MzidWriter: Write PSM score as cvParam search engine specific score instead of userParam score.
• io.percolator.PercolatorTabWriter: For PIN-style files: Use SpecId instead of PSMId and write PSMScore and ChargeN columns by default.

Fixed

• peptidoform: ProForma mass modifications are now correctly parsed within the rename_modifications function.
• io.maxquant.MSMSReader: Correctly parse empty Proteins column to None
• io.mzid.MzidReader: Set PSM.retention_time to None instead of float('nan') if missing from the PSM file.
• io.percolator.PercolatorTabReader: Correctly parse Percolator peptidoform notation if no leading or trailing amino acids are present (e.g. .ACDK. instead of K.ACDK.E).
• io.percolator.PercolatorTabWriter: ScanNr is now correctly written as an integer counting from the first PSM in the file.
• io.percolator.PercolatorTabWriter: If no protein information is present, write the peptidoform preceded by PEP_ to the Proteins column.

[codecov-commenter]

Codecov Report

Patch coverage: 37.68% and project coverage change: -0.05 :warning:

Comparison is base (638d19b) 41.52% compared to head (4dacc86) 41.48%.

:mega: This organization is not using Codecov’s GitHub App Integration. We recommend you install it so Codecov can continue to function properly for your repositories. Learn more

Additional details and impacted files

```diff @@ Coverage Diff @@ ## main #31 +/- ## ========================================== - Coverage 41.52% 41.48% -0.05% ========================================== Files 18 19 +1 Lines 1416 1468 +52 ========================================== + Hits 588 609 +21 - Misses 828 859 +31 ``` | [Impacted Files](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics) | Coverage Δ | | |---|---|---| | [psm\_utils/io/maxquant.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL2lvL21heHF1YW50LnB5) | `84.61% <ø> (ø)` | | | [psm\_utils/io/mzid.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL2lvL216aWQucHk=) | `29.71% <0.00%> (-0.29%)` | :arrow_down: | | [psm\_utils/io/percolator.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL2lvL3BlcmNvbGF0b3IucHk=) | `28.82% <20.00%> (-0.35%)` | :arrow_down: | | [psm\_utils/io/sage.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL2lvL3NhZ2UucHk=) | `41.02% <41.02%> (ø)` | | | [psm\_utils/peptidoform.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL3BlcHRpZG9mb3JtLnB5) | `31.47% <60.00%> (+0.78%)` | :arrow_up: | | [psm\_utils/io/\_\_init\_\_.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL2lvL19faW5pdF9fLnB5) | `37.03% <100.00%> (+0.78%)` | :arrow_up: | | [psm\_utils/psm.py](https://codecov.io/gh/compomics/psm_utils/pull/31?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics#diff-cHNtX3V0aWxzL3BzbS5weQ==) | `85.71% <100.00%> (ø)` | | Help us with your feedback. Take ten seconds to tell us [how you rate us](https://about.codecov.io/nps?utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics). Have a feature suggestion? [Share it here.](https://app.codecov.io/gh/feedback/?utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=compomics)

:umbrella: View full report in Codecov by Sentry.
:loudspeaker: Do you have feedback about the report comment? Let us know in this issue.

[lazear]

I just opened a fork to do this, should've checked here first! Thanks for implementing this

[RalfG]

@lazear, great to hear that you wanted to add Sage support! Feel free to take a look at the implementation and give feedback where needed. One thing I noticed was that the Sage documentation lists the spectrum_q column, while the actual output file seems to contain a spectrum_fdr column instead.

[lazear]

@lazear, great to hear that you wanted to add Sage support! Feel free to take a look at the implementation and give feedback where needed. One thing I noticed was that the Sage documentation lists the spectrum_q column, while the actual output file seems to contain a spectrum_fdr column instead.

Thanks for the catch - this is what I get for using ChatGPT to write docs :)