DiaNN icon indicating copy to clipboard operation
DiaNN copied to clipboard

--min-pep-len 6 but no peptides of lenght 6 are found

Open wolski opened this issue 1 year ago • 1 comments

We are interested in a protein that produces mostly short peptides. Therefore, we did try to make a DIANN search which includes peptides of length 6 by setting the parameter --min-pep-len 6 However, no peptides of length 6 are reported by DIANN.

( I am attaching the log) and the R code we used to determine the length of the peptides reported by DIANN.

report.log.txt

> xx <- readr::read_tsv("o32924_DIANN_GUI/report.tsv")
Rows: 1074785 Columns: 58                                                                                                                                                        
── Column specification ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: "\t"
chr (13): File.Name, Run, Protein.Group, Protein.Ids, Protein.Names, Genes, Modified.Sequence, Stripped.Sequence, Precursor.Id, First.Prot...
dbl (44): PG.Quantity, PG.Normalised, PG.MaxLFQ, Genes.Quantity, Genes.Normalised, Genes.MaxLFQ, Genes.MaxLFQ.Unique, Precursor.Charge, Q....
lgl  (1): Translated.Quality

ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
> xl <- xx$Stripped.Sequence |> unique() |> stringi::stri_length()
> xl |> table()
xl
   7    8    9   10   11   12   13   14   15   16   17   18   19   20   21   22   23   24   25 
2379 4883 6823 8145 8560 7933 7015 6219 5541 5097 4223 3571 3035 2454 1896 1376 1088  857  720 

wolski avatar Oct 26 '23 09:10 wolski