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

Any option like BWA MEM '-C'? #4370

Open
alshai opened this issue Aug 8, 2024 · 1 comment
Open

Any option like BWA MEM '-C'? #4370

alshai opened this issue Aug 8, 2024 · 1 comment

Comments

@alshai
Copy link

alshai commented Aug 8, 2024

Is there anything similar in VG to the -C option in BWA MEM?

From the BWA mem manual

-C Append append FASTA/Q comment to SAM output. This option can be used to transfer read meta information (e.g. barcode) to the SAM output. Note that the FASTA/Q comment (the string after a space in the header line) must conform the SAM spec (e.g. BC:Z:CGTAC). Malformated comments lead to incorrect SAM output.

Basically we have a comment in our FASTQ that is formatted like a SAM tag that we want to transfer over to the BAM output in VG.

If not already implemented, could I request this as a feature?

@adamnovak
Copy link
Member

We don't have this at the moment. We might be able to turn FASTQ comments into one of our internal annotations and then represent them in SAM output, but we'd have to pull the flag for this down into the general readers and writers somehow.

I'd recommend doing something like:

  • Pull out all the FASTQ header lines and sort by read name
  • Get all the SAM lines and sort by read name
  • Use Unix join to graft the FASTQ comments onto the ends of the corresponding SAM lines.

But that's definitely way more passes than just bringing the information through during mapping.

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

2 participants