search

ksahlin / strobealign

Aligns short reads using dynamic seed size with strobemers
MIT License
141 stars 17 forks source link

Add option to add FASTQ comment to SAM output. #378

Closed ghuls closed 9 months ago

ghuls commented 9 months ago

Add option to add FASTQ comment to SAM output.

BWA supports this with the -C option:

-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.

This is useful when cell barcode information was added to the FASTQ read names in the comment part (read name + space + SAM tags separated by TABs):

@A00305:504:HLMGMDRXY:1:2101:23023:1031 CR:Z:CAACGGCCACACATGT   CY:Z:FFFFFFFFFFFFFFFF   CB:Z:CAACGGCCACACATGT-1
GTCTTGGCTTTCTGTGCGGAAGTGGGGCTGGCTGGCATAGAATTCCTTTG
+
FFFFFFFF,:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFF
@A00305:504:HLMGMDRXY:1:2101:24415:1031 CR:Z:TGTAAGCAGACACAAA   CY:Z:FFFFFFFFFFFFFF,,   CB:Z:TGTAAGCAGACACAAT-1
ATGTCACACTCTGTGTCTTCTGCTAGCCCAGTCCTGTTTGGCAGCTCTAG
+
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF