Msconvert Capabilities
From SPCTools
Revision as of 23:29, 19 May 2010 Gsun (Talk | contribs) ← Previous diff |
Revision as of 23:29, 19 May 2010 Gsun (Talk | contribs) Next diff → |
||
Line 22: | Line 22: | ||
<br>X!Tandem search: Tandem Search completed | <br>X!Tandem search: Tandem Search completed | ||
<br>TPP:Results follow below: | <br>TPP:Results follow below: | ||
- | INFO: Reading ProteinProphet file 'interact-prob.prot.xml' | + | INFO: Reading ProteinProphet file 'interact-prob.prot.xml' |
- | PepP 7/ 58 0.121 | ProP 1 ( 6, 6) /proteomics/gsun/trapper-test/msconvert/mzML | + | PepP 7/ 58 0.121 | ProP 1 ( 6, 6) /proteomics/gsun/trapper-test/msconvert/mzML |
- | P threshold for FDR 0.01=0.25 | + | P threshold for FDR 0.01=0.25 |
- | Decoy results: | + | Decoy results: |
PP Incor FDR 0.01=2 | PP Incor FDR 0.01=2 | ||
Decoy IDs at FDR 0.01=0 | Decoy IDs at FDR 0.01=0 | ||
Line 33: | Line 33: | ||
At PepPro FDR 0.010 FDR based on decoys=-0.0000 | At PepPro FDR 0.010 FDR based on decoys=-0.0000 | ||
FDR after discard=-0.0000 | FDR after discard=-0.0000 | ||
- | For FDR 0.01, sensitivity=1.000 at P=0.25 | + | For FDR 0.01, sensitivity=1.000 at P=0.25 |
- | Peptide Count at FDR 0.01: | + | Peptide Count at FDR 0.01: |
Distinct peptides = 6 | Distinct peptides = 6 | ||
- | eptide Count at FDR 0.01: | + | Peptide Count at FDR 0.01: |
Distinct peptides = 6 | Distinct peptides = 6 | ||
Singletons = 5 (0.83) | Singletons = 5 (0.83) | ||
Doubletons = 1 (0.17) | Doubletons = 1 (0.17) | ||
- | Number of PSMs per MSrun: | + | Number of PSMs per MSrun: |
NEP_Plate1_5pmol1.mzML : 7 | NEP_Plate1_5pmol1.mzML : 7 | ||
Charge 3: 7 (FDR 0.01) Queries=58 | Charge 3: 7 (FDR 0.01) Queries=58 | ||
Line 47: | Line 47: | ||
Fraction semi-tryptics: 0.7143 | Fraction semi-tryptics: 0.7143 | ||
N missed cleavages 0: 7 (FDR 0.01) | N missed cleavages 0: 7 (FDR 0.01) | ||
- | Total Modifications at FDR 0.01: | + | Total Modifications at FDR 0.01: |
Modification: None: 7 | Modification: None: 7 | ||
- | Mass diffs for P > 0.25: Mean: 0.393142857142857 StDev: 1.03355058392118 | + | Mass diffs for P > 0.25: Mean: 0.393142857142857 StDev: 1.03355058392118 |
- | Mass diffs for P > 0.25: Min: -0.003 Max: 2.737 | + | Mass diffs for P > 0.25: Min: -0.003 Max: 2.737 |
- | -0.750 -> -0.250 0 | | + | -0.750 -> -0.250 0 | |
- | -0.250 -> 0.250 6 |************************************************** | + | -0.250 -> 0.250 6 |************************************************** |
- | 0.250 -> 0.750 0 | | + | 0.250 -> 0.750 0 | |
- | 0.750 -> 1.250 0 | | + | 0.750 -> 1.250 0 | |
- | 1.250 -> 1.750 0 | | + | 1.250 -> 1.750 0 | |
- | 1.750 -> 2.250 0 | | + | 1.750 -> 2.250 0 | |
- | 2.250 -> 2.750 1 |******** | + | 2.250 -> 2.750 1 |******** |
- | 2.750 -> 3.250 0 | | + | 2.750 -> 3.250 0 | |
- | ProteinProphet info: | + | ProteinProphet info: |
Proteins identified at P>0.9=1 | Proteins identified at P>0.9=1 | ||
Distinct modified peptides at P>0.25=6 | Distinct modified peptides at P>0.25=6 |
Revision as of 23:29, 19 May 2010
Testing of MSconvert is in progress and results will be posted here as they come.
Contents |
Agilent
Raw file type: .d directory
raw-to-mzXML converter: trapper
trapper
mzXML
Conversion: OK
X!Tandem search: Tandem Search completed
TPP:
msconvert
mzXML
Conversion: Conversion completed
X!Tandem search:Tandem Search completed
TPP:
mzML
Conversion: Conversion Completed
X!Tandem search: Tandem Search completed
TPP:Results follow below:
INFO: Reading ProteinProphet file 'interact-prob.prot.xml' PepP 7/ 58 0.121 | ProP 1 ( 6, 6) /proteomics/gsun/trapper-test/msconvert/mzML P threshold for FDR 0.01=0.25 Decoy results: PP Incor FDR 0.01=2 Decoy IDs at FDR 0.01=0 Decoy IDs at P<0.001=-1 Total IDs at P<0.001=47 Decoy fraction P<0.001=-0.0213 At PepPro FDR 0.010 FDR based on decoys=-0.0000 FDR after discard=-0.0000 For FDR 0.01, sensitivity=1.000 at P=0.25 Peptide Count at FDR 0.01: Distinct peptides = 6 Peptide Count at FDR 0.01: Distinct peptides = 6 Singletons = 5 (0.83) Doubletons = 1 (0.17) Number of PSMs per MSrun: NEP_Plate1_5pmol1.mzML : 7 Charge 3: 7 (FDR 0.01) Queries=58 Enzymatic Termini 1: 5 Enzymatic Termini 2: 2 Fraction semi-tryptics: 0.7143 N missed cleavages 0: 7 (FDR 0.01) Total Modifications at FDR 0.01: Modification: None: 7 Mass diffs for P > 0.25: Mean: 0.393142857142857 StDev: 1.03355058392118 Mass diffs for P > 0.25: Min: -0.003 Max: 2.737 -0.750 -> -0.250 0 | -0.250 -> 0.250 6 |************************************************** 0.250 -> 0.750 0 | 0.750 -> 1.250 0 | 1.250 -> 1.750 0 | 1.750 -> 2.250 0 | 2.250 -> 2.750 1 |******** 2.750 -> 3.250 0 | ProteinProphet info: Proteins identified at P>0.9=1 Distinct modified peptides at P>0.25=6 P threshold for protein group FDR 0.01=0.0000 Protein groups identified at FDR 0.01=0
Thermo
Raw file type: .RAW file
raw-to-mzXML converter: ReAdW
ReAdW
mzXML
Conversion: using remoteconvert.pl --convertWith readw, file successfully converted.
X!Tandem search:search ran successfully, slight differences with the msconvert mzXML file
TPP:TPP Search successful--calctppstat for mzXML by ReAdW similar to that of calctppstat for mzXML file by msconvert, slight differences, histograms are relatively different, numbers are relatively close to one another, through most numbers are different--a few number are the same, results are very close--see below(section msconvert: mzXML)
msconvert
mzXML
Conversion: Appears to be successful
X!Tandem search: Search completed, slight differences with the ReAdW converted mzXML file
TPP:Peptide prophet completed with warnings,protein prophet failed. Search completed----calctppstat for mzXML by ReAdW similar to that of calctppstat for mzXML file by msconvert, slight differences, histograms are relatively different, numbers are relatively close to one another, through most numbers are different--a few number are the same, results are very close, comparisons are below:
Calctppstat results for ReAdW file conversion
138 regis /proteomics/gsun/raw-test/ReAdW> /regis/sbeams/bin/calctppstat.pl --write --full --FDRthresh interact-prob.pep.xml INFO: Reading ProteinProphet file 'interact-prob.prot.xml' PepP 25/ 2621 0.010 | ProP 9 ( 16, 25) /proteomics/gsun/raw-test/ReAdW P threshold for FDR 0.01=0.98 Decoy results: PP Incor FDR 0.01=0 Decoy IDs at FDR 0.01=0 Decoy IDs at P<0.001=-1 Total IDs at P<0.001=1493 Decoy fraction P<0.001=-0.0007 At PepPro FDR 0.010 FDR based on decoys=-0.0000 FDR after discard=-0.0000 For FDR 0.01, sensitivity=0.160 at P=0.98 Peptide Count at FDR 0.01: Distinct peptides = 16 Singletons = 12 (0.75) Doubletons = 2 (0.12) Number of PSMs per MSrun: hui052703_15_1 : 25 Charge 1: 8 (FDR 0.01) Queries=146 Charge 2: 11 (FDR 0.01) Queries=2475 Charge 3: 6 (FDR 0.01) Queries=2475 Enzymatic Termini 1: 17 Enzymatic Termini 2: 8 Fraction semi-tryptics: 0.6800 N missed cleavages 0: 22 (FDR 0.01) N missed cleavages 1: 1 (FDR 0.01) N missed cleavages 2: 2 (FDR 0.01) Fraction of missed cleavages: 0.1200 Total Modifications at FDR 0.01: Modification: C[160]: 5 Modification: M[147]: 8 Modification: None: 14 Mass diffs for P > 0.98: Mean: 1.349 StDev: 0.981512735526137 Mass diffs for P > 0.98: Min: -0.154 Max: 3.230 -0.750 -> -0.250 0 | -0.250 -> 0.250 5 |***************************************** 0.250 -> 0.750 1 |******** 0.750 -> 1.250 6 |************************************************** 1.250 -> 1.750 3 |************************* 1.750 -> 2.250 5 |***************************************** 2.250 -> 2.750 3 |************************* 2.750 -> 3.250 2 |**************** 3.250 -> 3.750 0 | ProteinProphet info: Proteins identified at P>0.9=9 Distinct modified peptides at P>0.98=16 P threshold for protein group FDR 0.01=0.9541 Protein groups identified at FDR 0.01=8 Decoy IDs P>0 at protein group FDR 0.01=0 (0.0000) Decoy IDs P>0 at protein group FDR 0.03=0 (0.0000) Decoy IDs P>0 at protein group FDR 0.04=0 (0.0000)
Calctppstat results for msconvert file conversion
143 regis /proteomics/gsun/raw-test/msconvert> /regis/sbeams/bin/calctppstat.pl --write --full --FDRthresh interact-prob.pep.xml INFO: Reading ProteinProphet file 'interact-prob.prot.xml' PepP 18/ 2621 0.007 | ProP 10 ( 13, 23) /proteomics/gsun/raw-test/msconvert P threshold for FDR 0.01=0.98 Decoy results: PP Incor FDR 0.01=0 Decoy IDs at FDR 0.01=0 Decoy IDs at P<0.001=-1 Total IDs at P<0.001=1481 Decoy fraction P<0.001=-0.0007 At PepPro FDR 0.010 FDR based on decoys=-0.0000 FDR after discard=-0.0000 For FDR 0.01, sensitivity=0.110 at P=0.98 Peptide Count at FDR 0.01: Distinct peptides = 10 Singletons = 6 (0.60) Doubletons = 2 (0.20) Number of PSMs per MSrun: hui052703_15_1 : 18 Charge 1: 8 (FDR 0.01) Queries=146 Charge 2: 4 (FDR 0.01) Queries=2475 Charge 3: 6 (FDR 0.01) Queries=2475 Enzymatic Termini 1: 12 Enzymatic Termini 2: 6 Fraction semi-tryptics: 0.6667 N missed cleavages 0: 15 (FDR 0.01) N missed cleavages 1: 1 (FDR 0.01) N missed cleavages 2: 2 (FDR 0.01) Fraction of missed cleavages: 0.1667 Total Modifications at FDR 0.01: Modification: C[160]: 4 Modification: M[147]: 6 Modification: None: 10 Mass diffs for P > 0.98: Mean: 1.373 StDev: 1.09399618237854 Mass diffs for P > 0.98: Min: -0.154 Max: 3.230 -0.750 -> -0.250 0 | -0.250 -> 0.250 5 |************************************************** 0.250 -> 0.750 0 | 0.750 -> 1.250 2 |******************** 1.250 -> 1.750 3 |****************************** 1.750 -> 2.250 4 |**************************************** 2.250 -> 2.750 2 |******************** 2.750 -> 3.250 2 |******************** 3.250 -> 3.750 0 | ProteinProphet info: Proteins identified at P>0.9=10 Distinct modified peptides at P>0.98=13 P threshold for protein group FDR 0.01=0.9484 Protein groups identified at FDR 0.01=9 Decoy IDs P>0 at protein group FDR 0.01=0 (0.0000) Decoy IDs P>0 at protein group FDR 0.02=0 (0.0000) Decoy IDs P>0 at protein group FDR 0.03=0 (0.0000)
mzML
Conversion: Conversion seemed to work. TPP tests did not work with mzML file format
X!Tandem search:
TPP:TPP Search failed
Analyst
Raw file type: .wiff file
raw-to-mzXML converter: mzWiff
mzWiff
mzXML
Conversion: used /regis/sbeams/bin/remoteconvert.pl --default --convertWith mzWiff QS20061230_P0094_C1_08_AQv1.1.wiff, file converted, but not output files created or received. Zip file empty. Tried using ep server, new problem arose. Zhi says that mzWiff has many problems and some files it just cannot convert. This is probably what is happening.
X!Tandem search:
TPP:
msconvert
mzXML
Conversion:Tried /regis/sbeams/bin/remoteconvert.pl --default --convertWith msconvert QS20061230_P0094_C1_08_AQv1.1.wiff, could not connect with server. Also tried /regis/sbeams/bin/remoteconvert.pl --default --server th --convertWith msconvert QS20061230_P0094_C1_08_AQv1.1.wiff, also could not connect to server.
Tried using both actual file and symbolic link for .wiff file; same problem as above.
X!Tandem search:started runtandemsearch but it didn't complete
TPP: Once runtandemsearch completed, zztandempostproccessing.log file was full of warnings--WARNING: Mixture model quality test failed for charge (2+).
WARNING: Mixture model quality test failed for charge (3+).
mzML
Conversion: File appears to have been converted
X!Tandem search: Tandem search completed
TPP: TPP analysis had warnings--WARNING: Mixture model quality test failed for charge (2+).