Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Difference between the genome browser view of BAM produced by STAR and the result of Salmon #918

Open
NJU-Bio-Info opened this issue Mar 16, 2024 · 0 comments

Comments

@NJU-Bio-Info
Copy link

Hi, such a fantastic RNA-seq analyze tool Salmon! Recently, I used Salmon to analyze transcript level expression and the information of my analyze tool:

  • Salmon: v1.6.0 (install from bioconda)
  • STAR: 2.7.9a (install from bioconda)

And the detailed commands:

  • Salmon:
    image
    and the transcript index is constructed with hg38 transcripts fasta download from GENCODE:
    image

  • STAR:
    image

All commands were evaluated with Centos system.

After get the result file of Salmon, I tried to check the result in genome browser using BAM only with unique mapping reads produced by STAR. But I am confused with the following result:

According to salmon transript level TPM value, PPIE gene is tend to use longer isoform (the labelled value is TPM):
image
But when I visualized this gene in genome browser, I found that their are few reads aligned to specific region of the longest transcript:
image

So, I'm confused about this result and search form some possible reasons. There are some people saying that because Salmon uses transcript index rather than genome index like STAR, so it can rescue more reads. But I would like to ask you for more advices, great thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant