A Novel transcript without high-confidence reads support

[yuyun-zhang]

Hi, thank you for developing this amazing tool.

I used 9 long-read RNA-seq samples togather to generate a gtf file and I am very interested in a novel transcript (transcript32428.chr19.nic). When I load the bam files and the isoquant gtf file (annov29_new.transcript_models.gtf) in the IGV, I found that no read can support this novel transcript. The start site, end site, and splicing sites of the reads are not exactly the same as this novel transcript, or not mostly the same as this novel transcript.

I checked the file annov29_new.transcript_model_reads.tsv.gz and found 4 reads are recorded to support transcript32428.chr19.nic. However, these 4 reads are actually different from transcript32428.chr19.nic.

Here is the snapshot of reads and novel transcript transcript32428.chr19.nic. I marked these 4 reads and novel transcript . It seems transcript32428.chr19.nic is the result of the fusion of these 4 reads.

Could you please help me understand this result? Do you have any considerations for outputting results like this case?

Thank you!

Best, Yuyun