dedup STAR transcriptome file using umi_tools
1
0
Entering edit mode
5 weeks ago
ATRX ★ 1.1k

Hi,

I am interested in dedup the transcriptome output file (bulk RNA-seq) from STAR using umi_tools. I am using the following dedup function from umi_tools. Here is the command:

umi_tools dedup --paired
--stdin=B1-Cond1_Aligned.toTranscriptome.out.bam --log=B1-Cond1_dedup.txt --umi-separator=":" --output-stats=B1-Cond1_ > .dedup.bam

However, I am getting the following error:

ValueError: fetch called on bamfile without index

I don't think we can index transcriptome file from STAR. I read the umi_tools document and didn't much options for bulk RNA-seq libraries. Do you know what is the best way to dedup using umi_tools? Thanks!

rna-seq umi transcriptome umi_tools • 1.4k views
ADD COMMENT
1
Entering edit mode

How about deduping based on the genomic bam, then filtering the transcriptome bam based on the reads leftover. Maybe some extra work, but might be faster if you can't figure out how to index/dedupe the transcriptome bam directly.

ADD REPLY
0
Entering edit mode

Yes, you are right. I plan to do this. Thanks!

ADD REPLY
0
Entering edit mode

Very helpful, thanks a lot!

ADD REPLY
6
Entering edit mode
5 weeks ago

I don't know if any reason you shouldn't be able to index the bam produced by STAR, although you will need to sort them first.

ADD COMMENT
0
Entering edit mode

Thanks for the reply. STAR generates two different bam files. One is the genome based and the other one is the transcriptome based. I am able to dedup of genome based bam file but not for transcriptome based.

ADD REPLY
3
Entering edit mode

You should definitely be able to use samtools ti sort and then index the transcriptome bam. We've done this many times.

ADD REPLY
0
Entering edit mode

It worked. thanks a lot!

ADD REPLY
0
Entering edit mode

You can go ahead and accept @Ian answer to provide closure to this thread (green check mark).

ADD REPLY
0
Entering edit mode

Thank you, I didn't know about the green check mark. I did it now, thanks!

ADD REPLY

Login before adding your answer.

Traffic: 1847 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