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?
•
written
by
dzb •
1 answer
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:
Log in to answer this question.
-
Multiple unique reads sharing same UMI barcode
written by ww22runner •
-
VCF generation of UMI tagged data
written by Saurabh •
-
How Cell ranger count UMI from intronic region in the latest version
written by tien •
-
STAR solo paramters for 10X chromium single cell rna seq data alignment with R1 reads 150bp long
written by bp22 •
-
Cellranger output .bam file does Not contain barcodes ( CB and UB )
written by Gabriel
-
How can I align single RNA seq(scRNAseq) data by using STAR or another aligner?
written by dhkwnr95 •
-
Editing BAM file using pysam
written by sam •
-
STARsolo read1 length greater than cell barcode length and UMI length
written by from the mountains
-
split BAM file with same cell-barcode and UMI pair
written by newbinf •
-
Parsing BAM to genome feature
written by dzb •
Please post example lines from the two files.