STAR output processing
0
0
Entering edit mode
4.4 years ago
kintany ▴ 10

Hi all,

I'm running STAR to align RNA-seq reads separately to paternal and maternal genomes for subsequent analysis of allele-specific expression with GATK tools. So as an output I have two BAM files.

I want to keep reads aligned uniquely to only one genome and in cases where reads mapped uniquely to both genomes, use the alignment with the higher alignment quality.

Are there any existing tools for it? I can write a script doing this but I would prefer to use existing tool if there is one.

Thank you!

RNA-Seq ASE STAR • 2.6k views
ADD COMMENT
0
Entering edit mode

BBMap/BBSplit from the BBTools suite might be able to do this.

ADD REPLY
0
Entering edit mode

So you have two separate reference genomes? While BBSplit would certainly be an option (as indicated by @Sinji) multi-mapping is going to make this tricky. What exactly are you trying to do?

Will go ahead and tag : Brian Bushnell

ADD REPLY
0
Entering edit mode

Yes, I'm mapping to two separate genomes. Then I want to use ASEReadCounter from GATK for ASE statistics. It requires BAM file as input so I need to find a way to create this file from two BAM files that I have now

ADD REPLY
0
Entering edit mode

I don't think that strategy will work. ASEReadCounter uses a single BAM. The strand is based on the alleles of the reads ("tool calculates allele counts at a set of positions"). Even if you could use two BAMs or merge them as one, the paternal and maternal reads would not align to the same position.

ADD REPLY
0
Entering edit mode

one set of fastqs mapped to two REF genome = two BAM, isn't it ?

ADD REPLY
0
Entering edit mode

yes, that's what I'm saying

ADD REPLY

Login before adding your answer.

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