-
Notifications
You must be signed in to change notification settings - Fork 12
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
Prepend chromosomes with chr to match the SAM/BAM files. #132
base: main
Are you sure you want to change the base?
Conversation
Sorry, I'm terrible and added a separate commit on to this cosmetic one. If you like one and not the other feel free to ignore one of them. |
Y 9196964 9197108 Y-9196963-9197108 0 + Y-CCDS59533.1 Y-TSPY8 0.0 0.0 | ||
Y 9197216 9197296 Y-9197215-9197296 0 + Y-CCDS59533.1 Y-TSPY8 0.0 0.0 | ||
Y 9197992 9198013 Y-9197991-9198013 0 + Y-CCDS59533.1 Y-TSPY8 0.0 0.0 | ||
chr#{"group_id": null, "sample_id": "tiriwiro", "coverage_source": "/Users/robinandeer/Projects/clinical-genomics/qxt/bam/Sample_E2-QXT-prep2/sorted.bam", "cutoff": 10, "extension": 0} |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This leading "chr" must be a mistake, right? 😃
Thanks a lot @roryk - good catch! I'll be happy to test the PR locally after you've made the few suggested changes If you have time you could also update the |
apply_bed_stream takes a function and a BED stream applies the function to every interval in the BED stream Also fixes the test BED file not having chr prepended to the chromosome.
Superb Rory! Very nice touch with the rebase - I should keep learning more about how to do that And I wonder what it say of me if I'm blushing a little reading you last comment? 😅 Thanks! |
Hi Robin,
The demo was failing because the chromosome names didn't match up with the alignment files. This fixes that. Hope you are well.