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

CIGAR string #5

Open
fhach opened this issue Jun 25, 2020 · 1 comment
Open

CIGAR string #5

fhach opened this issue Jun 25, 2020 · 1 comment
Assignees
Labels
enhancement New feature or request

Comments

@fhach
Copy link
Collaborator

fhach commented Jun 25, 2020

Any plans on adding CIGAR strings to SAM output so I can load the alignments in a genome browser?

Originally posted by @yhoogstrate in #4 (comment)

@fhach fhach added the enhancement New feature or request label Jun 25, 2020
@HosseinAsghari HosseinAsghari changed the title CIGRA string CIGAR string Jun 25, 2020
@yhoogstrate
Copy link

yhoogstrate commented Jun 26, 2020

Since the CIGAR's N blocks do not support negative numbers (for back splicing), it is in practice impossible to use one SAM segment to describe the alignment of a back spliced read entirely.

I believe that it is common practice to split a backspliced read in two and use two aligned segments of which the counter-part is soft-clipped.
Of those two segments, one (typically last part of the read) is set to 'supplementary/secondary alignment'. Both split alignments/segments then should have an 'SA' tag in place to reference each other. (http://samtools.github.io/hts-specs/SAMtags.pdf ; page 4)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

3 participants