Skip to content
Snippets Groups Projects
Commit 0ad05745 authored by Langella Olivier's avatar Langella Olivier
Browse files

fix adding modification at position -1 in pwiz reader

parent 4b2e43a6
No related branches found
No related tags found
No related merge requests found
...@@ -247,7 +247,8 @@ getQVariantDoubleParam(pwiz::data::ParamContainer *item, pwiz::cv::CVID param) ...@@ -247,7 +247,8 @@ getQVariantDoubleParam(pwiz::data::ParamContainer *item, pwiz::cv::CVID param)
void void
IdentificationPwizReader::read( IdentificationPwizReader::read(
IdentificationDataSource *p_identification_data_source, Project *p_project, IdentificationDataSource *p_identification_data_source,
Project *p_project,
IdentificationGroup *p_identification_group) IdentificationGroup *p_identification_group)
{ {
qDebug() << "IdentificationPwizReader::read begin"; qDebug() << "IdentificationPwizReader::read begin";
...@@ -314,9 +315,29 @@ IdentificationPwizReader::read( ...@@ -314,9 +315,29 @@ IdentificationPwizReader::read(
pappso::pappso_double mass_delta = pappso::pappso_double mass_delta =
mod_ptr.get()->monoisotopicMassDelta; mod_ptr.get()->monoisotopicMassDelta;
peptide.addAaModification( if(location > 0)
Utils::guessAaModificationPbyMonoisotopicMassDelta(mass_delta), {
location - 1); location = location - 1;
}
try
{
peptide.addAaModification(
Utils::guessAaModificationPbyMonoisotopicMassDelta(
mass_delta),
location);
}
catch(pappso::PappsoException &error)
{
throw pappso::PappsoException(
QObject::tr("ERROR adding modification mass delta %1 to "
"peptide %2 at "
"position %3 :\n %4")
.arg(mass_delta)
.arg(peptide.toAbsoluteString())
.arg(location)
.arg(error.qwhat()));
}
// std::vector<CVParam> cvParams; // std::vector<CVParam> cvParams;
for(const pwiz::data::CVParam &param : mod_ptr.get()->cvParams) for(const pwiz::data::CVParam &param : mod_ptr.get()->cvParams)
{ {
......
...@@ -911,16 +911,12 @@ MascotDatParser::setVariableModifications(PeptideXtpSp &peptide_sp, ...@@ -911,16 +911,12 @@ MascotDatParser::setVariableModifications(PeptideXtpSp &peptide_sp,
[variable_modifications_string.mid(i, 1).toUInt() - 1]; [variable_modifications_string.mid(i, 1).toUInt() - 1];
try try
{ {
if(i == 0) if(i > 0)
{ {
peptide_sp.get()->addAaModification( i = i - 1;
mascot_modification.modification, 0);
}
else
{
peptide_sp.get()->addAaModification(
mascot_modification.modification, i - 1);
} }
peptide_sp.get()->addAaModification(
mascot_modification.modification, i);
} }
catch(pappso::PappsoException &error) catch(pappso::PappsoException &error)
{ {
...@@ -929,7 +925,7 @@ MascotDatParser::setVariableModifications(PeptideXtpSp &peptide_sp, ...@@ -929,7 +925,7 @@ MascotDatParser::setVariableModifications(PeptideXtpSp &peptide_sp,
"position %3 :\n %4") "position %3 :\n %4")
.arg(mascot_modification.modification->getAccession()) .arg(mascot_modification.modification->getAccession())
.arg(peptide_sp.get()->toAbsoluteString()) .arg(peptide_sp.get()->toAbsoluteString())
.arg(i - 1) .arg(i)
.arg(error.qwhat())); .arg(error.qwhat()));
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment