search

gpertea / stringtie

Transcript assembly and quantification for RNA-Seq
MIT License
385 stars 78 forks source link

StringTie gives incorrect 5 end. #380

Open Pentayouth opened 2 years ago

Pentayouth commented 2 years ago

Dear gpertea and all Stringtie maintainer,

I performed single sample assembly with StringTie however the reconstructed 5 end is incorrect even the read coverage support in bigwig is obvious.

The parameters I used are --rf -p 1 -j 2 -c 2 -s 5 -f 0.05 (without -G because I perfer de novo assembly)

Here is the example of gene KRT18: 图片 The grey track is gencode v40 reference gtf. The red arrow denotes the transcript which I think should be reconstructed by StringTie

图片 However, in all single sample assemblies, StringTie gave a shorter, more downstream 5 end.

I couldn't figure out why could it happen. Is there anything to do with my parameter?

Best regards, Wang