Perform a theoretical digest of a protein database.trans_proteomic_pipeline
digestdb -m$missed_cleavages -l$min_mass -h$max_mass -r${str($enzyme).split(":")[0]} -n${str($enzyme).split(":")[1]} $sequences > $output
**What it does**
Produces a tabular description of peptides resulting from theoretical
digest of a given protein database by the selected enzyme agent using
the digestdb program distributed with the Transproteomic Pipeline.
** Output **
Results are spit out to stdout in a tab-delimited format.
Asterisks (*) in sequence are treated as proper break points
Output columns:
- peptide length
- protein reference
- peptide mass
- previous amino acid before peptide
- peptide sequence
- next amino acid after peptide
- peptide start location
- peptide end location
- pI
------
**Citation**
For the underlying tool, please cite `Proteomics. 2010 Mar;10(6):1150-9. A guided tour of the Trans-Proteomic Pipeline. Deutsch EW, Mendoza L, Shteynberg D, Farrah T, Lam H, Tasman N, Sun Z, Nilsson E, Pratt B, Prazen B, Eng JK, Martin DB, Nesvizhskii AI, Aebersold R. PMID 20101611`
If you use this tool in Galaxy, please cite Chilton J, et al. https://bitbucket.org/galaxyp/galaxyp-toolshed-digestdb