RPKM with bowtie output
1
0
Entering edit mode
24 months ago
flogin ▴ 270

Hello guys,

There is a way to evaluate RPKM direct with bowtie 2 output (e.g. bam or sorted bam)? like an R package or some scripts to do this?

bowtie rpkm rna count • 919 views
ADD COMMENT
2
Entering edit mode

Is this RNA-seq? Are you aware bowtie2 is not splice-aware? What organism? Do you have a GTF file or list of reference peaks (if not RNA-seq)?

ADD REPLY
0
Entering edit mode

It's a mapping of RNA-seq against a set of contigs, I just have the contigs, rna-seq files and bam files.

ADD REPLY
4
Entering edit mode
24 months ago

RSEM will take a bam file as input and return raw counts, RPKM and TPM

ADD COMMENT
0
Entering edit mode

thanks swbarnes2, I'll check this.

ADD REPLY
0
Entering edit mode

Hello, Can you please provide the steps used for calculating expression using RSEM, after building the reference. The script is throwing errors.

ADD REPLY
0
Entering edit mode

Please read the manual and then come back with more specific questions while opening a new question, providing the necessary details.

ADD REPLY
0
Entering edit mode

The script is throwing errors.

How do you expect us to help if you don't tell us what's going wrong?

ADD REPLY
0
Entering edit mode

I used the same commands given in the manual, with the BAM files I have,

software/RSEM-1.2.25/rsem-calculate-expression -p 8 --paired-end \ --bam \ --estimate-rspd \ --append-names \ --output-genome-bam \ exp/LPS_6h.bam \ ref/mouse_ref exp/LPS_6h

Error:

Warning: The SAM/BAM file declares less reference sequences (955) than RSEM knows (41078)! Please make sure that you aligned your reads against transcript sequences instead of genome. RSEM can not recognize reference sequence name 1!

"rsem-parse-alignments /home/user/data/RSEM/Rat /home/user/data/RSEM/Rat.temp/Rat /home/user/data/RSEM/Rat.stat/Rat /home/ec2-user/data/Rat/*.bam 3 -tag XM" failed! Plase check if you provide correct parameters/options for the pipeline!

ADD REPLY

Login before adding your answer.

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