This is a test version of Biostars. For the public version, visit https://www.biostars.org.
biomaRt correct annotation but...

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
• 2,520 views • link
written by Mozart Mozart

1 answer

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

• 0 views • link
written by haroldpimentel • haroldpimentel

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

• 0 views • link
written by Mozart Mozart

Log in to answer this question.

More posts like this