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

qbuf problem #59

Open
wuyiwen-genomics-cn opened this issue Sep 23, 2024 · 0 comments
Open

qbuf problem #59

wuyiwen-genomics-cn opened this issue Sep 23, 2024 · 0 comments

Comments

@wuyiwen-genomics-cn
Copy link

wuyiwen-genomics-cn commented Sep 23, 2024

Hi, I just noticed the qbuf in features.rs -> extract_features -> reads[qid as usize].seq.get_subseq does not work as it supposed. If I am not wrong, qbuf should be assigned with seq value in range [qstart, qend). In decode function, it is assigned with correct seq value. However, when I print qbuf outside "match ow.overlap.strand ", qbuf has some different values at some positions.

For example, in decode function, qbuf is assigned with 84, 84, 65, 71, 65, but printed qbuf is 84, 84, 65, 67, 65.
At position 4, the value is 67 instead of 71.

May I know why it happened? I did not find any extra function called.
Thank you!

@wuyiwen-genomics-cn wuyiwen-genomics-cn changed the title qbuf qbuf problem Sep 23, 2024
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

1 participant