I used tophat-cufflinks pipeline for RNA-seq analysis.
In the output file, there is only temporarily test ids such as XLOC_000001.
Since I use GRCh38 annotation file from NCBI and there are ensemble ids from gtf file.
How can I keep the transcript ids such as ensemble ids during cufflinks?
Or how can I retrieve ensemble ids from the output files of cuffdiff?
I did search and found this useful post from seqanswers:
http://seqanswers.com/forums/showthread.php?t=51071
I save this script as "m.py" run it. then it shows that
File "m.py", line 12 for line in fh: ^ IndentationError: expected an indented block
I am a newbie for python, anyone can help me ?
Thanks so much~~
1 answer
The cuffdiff output should have "nearest gene" annotation, as an XLOC code can encompass more than one gene, it's a sliding window of sorts that arbitrarily* chooses sections of the genome to quantify against.
*It's probably not arbitrarily, it's just not documented and I haven't looked at the code.
Log in to answer this question.