Question: Transcripts ID matching in tximport for edgeR
gravatar for Sharon
8 months ago by
Sharon340 wrote:

I am using tximport to convert transcripts counts from Salmon to gene counts for edgeR. Here is how my Salmon counts look like:

Name    Length  EffectiveLength TPM     NumReads
NM_000014.5     4945    4714.81 16.5914 4870.8
NM_000015.2     1317    962.778 0.0333619       2
NM_000016.5     2623    1865.02 5.28102 613.276
NM_000017.3     1964    1925.07 1.67686 201

But my tx2gene is like that, with no version (.x) after the transcript id:


tximport throws the following error :

None of the transcripts in the quantification files are present in the first column of tx2gene. Check to see that you are using the same annotation for both

My question is: Is it safe to remove the version .x numbers from my Salmon counts to match tx2gene file? The transcripts in Salmon will match tx2gene if I removed the version .x. Like for example convert NM_000014.5 in Salmon counts to NM_000014?


edger rna-seq salmon tximport • 553 views
ADD COMMENTlink modified 8 months ago by WouterDeCoster30k • written 8 months ago by Sharon340
gravatar for WouterDeCoster
8 months ago by
WouterDeCoster30k wrote:

Yes, safe.

You are looking for the ignoreTxVersion option of tximport:

logical, whether to split the tx id on the ’.’ character to remove version information, for easier matching with the tx id in gene2tx (default FALSE)

ADD COMMENTlink written 8 months ago by WouterDeCoster30k

Thanks a lot DeCoster. Works well now.

ADD REPLYlink written 8 months ago by Sharon340

You're welcome ron. Glad to help.

ADD REPLYlink written 8 months ago by WouterDeCoster30k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1278 users visited in the last hour