Converting STAR Gene-level alignment to TPM expression
0
0
Entering edit mode
4 months ago
CTLong ▴ 90

Hi,

I have recently performed gene-level alignment with STAR on 20 samples with the parameter --quantMode GeneCounts and --outSAMtype BAM SortedByCoordinate. I have the output files ReadsPerGene.out.tab and Aligned.sortedByCoord.out.bam.

From this, how can I generate reliable TPM values with either the sorted BAM files or the raw counts? I saw posts indicating that TPM values can be manually calculated using the gene length from the GTF/GFF3 file. However, I am hesitant in doing this, as there might be additional factors or parameters that have been overlooked.

Currently, my last resort will be to realign using STAR again, but on transcriptomeSAM mode and then use it as input to cufflink or salmon. Just wondering if there is a more convient and reliable approach in doing this with gene-level STAR alignments. Thanks!

RNA-seq • 347 views
ADD COMMENT
2
Entering edit mode

Ask yourself whether your analysis really needs the added precision in terms of accurate gene length from salmon/kallisto/similarTools rather than just using the length from STAR. If so then use this option to get accurate gene length estimates.

ADD REPLY
1
Entering edit mode

Why is that a last resort? That will generate the most reliable TPM values.

ADD REPLY

Login before adding your answer.

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