Question: biomaRt correct annotation but...
0
gravatar for Mozart
2.9 years ago by
Mozart200
Mozart200 wrote:

Hi everyone, I am a bit confused about a painful step of my pipeline using Sleuth. I am trying to move from the transcript level to the gene level and that's why I am using biomaRt. PS: please pardon my ignorance in term of bioinformatic jargon

rna-seq • 1.3k views
ADD COMMENTlink modified 3 months ago • written 2.9 years ago by Mozart200
3
gravatar for haroldpimentel
2.9 years ago by
haroldpimentel110 wrote:

Hey there,

Do you get results after running a test?

Basically, there is some weirdness between some of the ensembl annotations. Your transcriptome probably says something like 'ENSMUST00000148678.2' for that specific transcript. When we try to merge that with the target_mapping, it ends up being empty because the target_id is 'ENSMUST00000148678' (note the missing .2). It's basically just telling you that we removed .N and it seems to have fixed it. When we do that and there is no intersection mapping still, then the program will completely error.

If you still get results after running a test and calling sleuth_results(), then everything should be fine.

Best,

Harold

ADD COMMENTlink written 2.9 years ago by haroldpimentel110

Thanks a lot Harold, actually everything should be fine even if I am getting the error.

ADD REPLYlink written 2.9 years ago by Mozart200
Please log in to add an answer.

Help
Access

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