This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using Pbase bioconductor package to map protein substitutions to chromosome positions

I'm using trying to use this tutorial to map certain protein positions to genomic positions. The whole tutorial is focused on using the Proteins class from the Pbase package, yet they don't show how to construct one. The docs ([1], [2]) have literally no information about the object constructor. They only show how you can import all the information from an mzIdentML (mzid) file:

## Create a Protein object reading all proteins from a fasta file.
    fastaFiles <- list.files(system.file("extdata", package = "Pbase"),
                             pattern = "fasta", full.names = TRUE)
p <- Proteins(fastaFiles)`

## Add indentification data
idfile <- system.file("extdata/Thermo_Hela_PRTC_selected.mzid",
                      package = "Pbase")
p <- addIdentificationData(p, idfile)

This is the most useless piece of documentation I've ever read in my life. And I don't have an mzid file, though I do have all the needed information: UniProt IDs with corresponding Ensembl transcript IDs, peptide coordinates, gene names, etc. Do you have any ideas on how to construct the object with all the data without an mzid file? Thank you in advance.

r mzid pbase bioconductor

Did you check BioConductor's mzID? Does this link help?

0 answers

No answers yet.

Log in to answer this question.