Skip to content

a big bug? the tmpClyJgf exists, but it is empty, so can not read. how should I reviese it? #30

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

Open
kobejamescurry opened this issue Feb 25, 2019 · 18 comments

Comments

@kobejamescurry
Copy link

kobejamescurry commented Feb 25, 2019

thanks a lot, urgent for help

@ambj
Copy link
Owner

ambj commented Feb 26, 2019

This seems as an error occurring from NetMHCpan.
try running MuPeXI with the -t (temp file option) and look at the NetMHCpan output.

@kobejamescurry
Copy link
Author

I ran with that -f -t -n
this use the allele HLA-A01:01,HLA-B08:01 you supplied in the manual, though there is error in pepmatch * file , it can give the result TCRBOA1.mupexi
image

but when I use other alleles, the pepmatch * file has the errors and there can not give the result TCRBOA1.mupexi

if I use the example.vcf you supplied, the pepmatch * file is ok, also it can give the result TCRBOA1.mupexi
image

@ambj
Copy link
Owner

ambj commented Feb 26, 2019

what about the NetMHCpan output file, any error reported here.
Just to clarify - the '*' is not allowed in the input HLA alleles.

@kobejamescurry
Copy link
Author

after replace , give the result TCRBOA1.mupexi
but all the pepmatch_
is still wrong, HOW to solve this, thanks a lot

image

this is the netMHCpan_EL.txt
image

this is the netMHCpan_BA.txt
image

@ambj
Copy link
Owner

ambj commented Feb 26, 2019

What is the MuPeXI output error?
Do you have permission to the /tmp dir?

@kobejamescurry
Copy link
Author

when I run the command, just through python MuPeXI.py, the stdout is like following , gives no error
image

@kobejamescurry
Copy link
Author

permissoin is 600
image

@ambj
Copy link
Owner

ambj commented Feb 26, 2019

If MuPeXI runs without error and you get the output file and the corresponding normal peptide sequences for frame-shifts and indels - then it is not clear to me what your problem is?
The error regarding HLA-B40:02 is solved?

@kobejamescurry
Copy link
Author

in your MuPeXI.py code, you change the file permission

@kobejamescurry
Copy link
Author

it seems that The error regarding HLA-B40:02 is solved after replace *, I need to test more

@kobejamescurry
Copy link
Author

as the ower, I have R and W, it should be totally ok for the file,

@kobejamescurry
Copy link
Author

kobejamescurry commented Feb 26, 2019

1 in the final result, sometimes can has FS/X in the column of Amino_Acid_Change , it means frame-shifts and indels or other things , thanks a lot.

the following two columns all are FS/X, X means random amino, is FS short for frame-shifts, if so, if it is
indel, what should be in the column of Amino_Acid_Change

VPPAGTWAL VKPAPTWAL
PESTPSTL NTSTPSTL

2 without getting the right pepmatch * file, but give the .mupexi file, is the .mupexi file can be used?

thanks a lot

@kobejamescurry
Copy link
Author

kobejamescurry commented Feb 26, 2019

here you change the permission, in your MuPeXI.py
image

I think the permission maybe can only changed by you, because all the file generated by your procedure, thanks a lot

@kobejamescurry
Copy link
Author

kobejamescurry commented Feb 26, 2019

I tried to revise the origin code, I make the tmp_dir in the my own defined out_dir, and comment the # clean_up(tmp_dir), so I can keep al the tmp file.
but still get the same error, I really can not undestand.

I just ran python MuPeXI.py, the stdout is ok, but the inner is wrong,
I guess you leak some extreme conditions

image

image

image

the tmpClyJgf exists, but it is empty, so can not read. how should I reviese it ?
image

@kobejamescurry kobejamescurry changed the title Allele "HLA-B*40:02" not stated in NetMHCpan output a big bug? the tmpClyJgf exists, but it is empty, so can not read. how should I reviese it? Feb 26, 2019
@kobejamescurry
Copy link
Author

here is another extreme condition you may not considered

when there is no missense, fs , deletion, insertion, you lose this

@kobejamescurry
Copy link
Author

kobejamescurry commented Feb 27, 2019

2 without getting the right pepmatch * file, but give the .mupexi file, is the .mupexi file can be used?

like following
image

thanks a lot

@kobejamescurry
Copy link
Author

I do not know how the tmp file is used, hope you can give me help , thanks a lot
image

@ssudhaman
Copy link

ssudhaman commented Oct 6, 2020

Hi,
I was wondering if you had a fix for this yet, I seem to be getting the same error on my "pepmatch.txt" files.
"Error. Cannot read peplist from /tmp/MuPeXI_G_Iql4/tmpibkyxs"

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

3 participants