Question: Cufflinks Advice
0
gravatar for Michael Gooch
7.5 years ago by
Michael Gooch10 wrote:
I attempted to run cufflinks on some RNA sequencing data. It seemed to complete without any errors, but the output files were empty. I am trying to figure out if I did something wrong or whether my data needs some additional processing before cufflinks will be able to use it. (Or whether the data is unsuitable for cufflinks.) The data is paired end reads. M. Gooch
rna-seq cufflinks • 821 views
ADD COMMENTlink modified 7.5 years ago by David Matthews630 • written 7.5 years ago by Michael Gooch10
0
gravatar for David Matthews
7.5 years ago by
United Kingdom
David Matthews630 wrote:
Hi, As a first guess I would say that your chromosome names do not match somewhere along the line. If you look at your sam file and the fasta of the genome you are working with (and the gtf file as well if you are using it) you may find, for example, one refers to chromosome 1 as "chr1" whilst the other refers to chromosome 1 as "1" or even "Chr1" or some other way of referring to the chromosome - any of these mismatches can cause you to get an empty output. If you are using a built in index it may be you need to change your chromosome names from "1" to "chr1" for example. Amazingly, the names of human chromosomes are apparently not yet standardised across all databases for the human genome (and I presume this may be the case for other genomes as well). Best Wishes, David. __________________________________ Dr David A. Matthews Senior Lecturer in Virology Room E49 Department of Cellular and Molecular Medicine, School of Medical Sciences University Walk, University of Bristol Bristol. BS8 1TD U.K. Tel. +44 117 3312058 Fax. +44 117 3312091 D.A.Matthews@bristol.ac.uk
ADD COMMENTlink written 7.5 years ago by David Matthews630
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 16.09
Traffic: 178 users visited in the last hour