Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Failing trypsin digestion #15

Open
gwajnberg opened this issue Mar 11, 2020 · 2 comments
Open

Failing trypsin digestion #15

gwajnberg opened this issue Mar 11, 2020 · 2 comments

Comments

@gwajnberg
Copy link

This is the output:

NovelJunc-chr14-106175001-106331409|chr=14|type=novel|site1=106175001|site2=106331409|reads=74|frame=1-
WYFDLWGRGTLVTVSSASPTSPKVFPLSLCSTQ

where this sequence should be cleaved in three :
WYFDLWGR = GTLVTVSSASPTSPK=VFPLSLCSTQ

I tried using with and without --no_missed_cleveage

those were teh parameters : quilts.py --output_dir --genome --proteome --somatic -- germline --junction --junction_file_type star [--no_missed_clevage]

@ekawaler
Copy link
Owner

When there is a proline on the C-term side of a K/R trypsin will not cleave. Also, make sure you spell the parameter correctly (--no_missed_cleavage)

@gwajnberg
Copy link
Author

sorry i didn't misspelled in the command line !

I was trying to modify the code the trypsinize function, but I am still getting the same result .

in the condition of line 1048 I removed the " and next_letter != P"

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants