Add tags to BAM/SAM file
1
0
Entering edit mode
4.5 years ago
dzb • 0

I'm working with some RNA-seq data. I have alignments done in STAR and the resultant BAM file. I'd like to annotate this BAM alignment data with custom tags using data that are stored in a separate file. The data in the second file contain are a read ID, a barcode, and a UMI. I want to add the barcode and UMI to all reads in the BAM file that match the read ID in the second file.

To summarise:

First file: BAM output from STAR

Second file: Read ID (matching those in STAR BAM file), UMI, barcode.

How do I get the UMI and barcode in file 2 tagged onto the reads in the BAM file?

Intensive Google and forum searching have yielded little info about this but I have a feeling there's a simple answer. Can anyone help?

alignment RNA-Seq sequencing BAM SAM • 5.5k views
ADD COMMENT
0
Entering edit mode

Please post example lines from the two files.

ADD REPLY
3
Entering edit mode
4.5 years ago

I don't think there's a simple answer for this question, but if you want to write a script you might find simplesam (https://github.com/mdshw5/simplesam) useful:

ADD COMMENT
0
Entering edit mode

This looks like an excellent starting point, thanks for sending it my way.

ADD REPLY

Login before adding your answer.

Traffic: 2255 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6