Entering edit mode
Hi, I obtained BAM files from minimap2 aligner and I wish to conduct the SGSeq pipeline with them.
However, the getBamInfo()
requires custom tag from the BAM files. The BAM files produced with minimap2 uses "ts" tag instead.
Custom tag 'XS' not found in BAM file:
''
Spliced alignments must include a custom tag 'XS'.
Compatible BAM files can be obtained with an alignment
program for RNA-seq data (e.g. GSNAP, HISAT, STAR, TopHat).
The recommended aligner may not suitable for my case. How could I fix this? I appreciate your help. Thank you very much.