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

Change to support CRAM #6

Closed
4 tasks done
keiranmraine opened this issue Feb 25, 2016 · 1 comment
Closed
4 tasks done

Change to support CRAM #6

keiranmraine opened this issue Feb 25, 2016 · 1 comment
Assignees

Comments

@keiranmraine
Copy link
Contributor

keiranmraine commented Feb 25, 2016

If VerifyBamId adds CRAM support remove dependancy on Bio::DB::Sam and switch to Bio::DB::HTS.

This will then be fully BAM/CRAM agnostic.

Modify the perl wrapper to sample the relevant locations from an input CRAM to a temporary BAM:

samtools view -F 3844 -L loci.bed -q 10 -b orig.cram | tee new.bam | samtools index - > new.bam.bai

Make the applied flags obvious in code:

4 - read unmapped (0x4)
256 - not primary alignment (0x100)
512 - read fails platform/vendor quality checks (0x200)
1024 - read is PCR or optical duplicate (0x400)
2048 - supplementary alignment (0x800)
@keiranmraine
Copy link
Contributor Author

Investigate suggestion here:

statgen/verifyBamID#2

@keiranmraine keiranmraine self-assigned this Feb 14, 2018
keiranmraine added a commit that referenced this issue Feb 14, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant