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

Regarding output of nanosplicer (missing BED file) #11

Open
PRIYANKA-22091995 opened this issue Dec 4, 2023 · 4 comments
Open

Regarding output of nanosplicer (missing BED file) #11

PRIYANKA-22091995 opened this issue Dec 4, 2023 · 4 comments
Labels
bug Something isn't working

Comments

@PRIYANKA-22091995
Copy link

Hello Team,

I have got the nanosplicer output, but it is only the tsv file, there was no bed file generated. What is the reason for that.
The other question i have from the tsv file, based on the tsv file headers, should i consider the JAQ score or the SIQ score for the prediction of the splicing junction.
It would be nice to have your comments on the same.

Thanks & Regards
Priyanka Roy

@youyupei
Copy link
Member

youyupei commented Dec 4, 2023

Hi @PRIYANKA-22091995,

Thanks for being interested in NanoSplicer. For the bed output, would you be about to run the example script (script.sh under example) and check whether you can get the .bed. I can't be certain at the moment without further information but if you could get the bed from the example, I assume the reason you didn't get the .bed was because the program has NOT been finished successfully, you could send me any command line output/error message you got .

For the downstream usage, generally results with higher JAQ , higher SIQ, and higher "best_prob" tend to be more accurate. When you have enough reads, I would recommend filtering some results to get more confident results (usually JAQ > 0.95 SIQ > -0.4 and best_prob>0.95 will mean confident results, but you could adjust the number based on your depth).

Regards,
Yupei

@youyupei youyupei changed the title Regarding output of nanosplicer Regarding output of nanosplicer (missing BED file) Dec 5, 2023
@youyupei youyupei added the bug Something isn't working label Dec 5, 2023
@PRIYANKA-22091995
Copy link
Author

PRIYANKA-22091995 commented Dec 5, 2023 via email

@youyupei
Copy link
Member

youyupei commented Dec 7, 2023

I would highly recommend putting an SIQ threshold. It indicates whether or not the squiggle used to characterize the junctions is informative.
However the threshold doesn't have to be -0.4, I think you could relax it to ~-0.8.

@PRIYANKA-22091995
Copy link
Author

Thank you so much for making it more clear.

Thanks & Regards
Priyanka Roy

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants