TransWikia.com

BAM file filteing to remain best isoform

Bioinformatics Asked by user977828 on January 10, 2021

I ran HiSat2, MarkDuplicate, removed reads with the lower quality score than 40 and finally only kept properly paired reads.

After the BAM filtering steps, I used the Scallop results with TransDecoder. Next, I used bamCoverage to create RNA-Seq profile. I noticed that not always the first isoform is the best one compared to the RNA-Seq profile as could be seen below a few examples:

1.
Screen Shot 2020-08-22 at 4 35 26 PM

2.
Screen Shot 2020-08-22 at 4 38 45 PM

3.
Screen Shot 2020-08-22 at 4 51 37 PM

4.
Screen Shot 2020-08-22 at 4 58 42 PM

5.
Screen Shot 2020-08-22 at 5 04 34 PM

By any chance, is there a way to keep only the best isoform which is in red box?

Thank you in advance,

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP