vgteam / vg

tools for working with genome variation graphs
https://biostars.org/tag/vg/
Other
1.09k stars 193 forks source link

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

Open alshai opened 1 month ago

alshai commented 1 month ago

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 commented 2 weeks ago

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:

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