StringTie inconsistent expression value

0

Entering edit mode

6.8 years ago

nand ▴ 10

I'm trying to find novel isoforms from RNA-Seq. I have six samples, 3 control and 3 treatment. I made a merged gtf file after running StringTie with gencode annotation. Then the merged annotation was used to estimate FPKM in all samples. Then novel transcripts expression was tabulated as shown below. The left three columns are control and the right three are treatment. Why is the fpkm absolute zero for in some cases, when it is very high for the same transcript in other samples.

enter image description here

RNA-Seq • 1.8k views

ADD COMMENT • link updated 6.8 years ago by GenoMax 148k • written 6.8 years ago by nand ▴ 10

1

Entering edit mode

Have you loaded the BAM into IGV and gone into that locus to check? Most of the times, inspecting in IGV resolves doubt about differences in expression estimate (from assembler).

ADD REPLY • link 6.8 years ago by Amitm ★ 2.3k

0

Entering edit mode

I loaded BAM and merged GTF to IGV. I am seeing many fragments aligned to that transcript region. But FPKM is absolutely zero.

ADD REPLY • link 6.8 years ago by nand ▴ 10

Login before adding your answer.