Question: Multiple mapping to one gene ID using Cufflinks in Galaxy
0
gravatar for mgtrax10
3.5 years ago by
mgtrax10 • 0
United States
mgtrax10 • 0 wrote:

I am currently running RNA-seq analysis on rat mRNA in Galaxy using a flow that uses TopHat2 then cufflinks to determine FPKMs.  When I run cufflinks, however, some of the geneIDs get mapped to multiple loci on the genome with varying hits as shown below.  

NM_001008876   -   -   NM_001008876    -    -    chr10:4284438-4326284    -    -    24.3459    20.5645    25.0729    OK
NM_001008876   -   -   NM_001008876    -    -    chr13:40006284-40007070  -    -    7.2634    4.27328    7.56582    OK
NM_001008876   -   -   NM_001008876    -    -    chr17:51085969-51086776  -   -   17.9851    11.5232    16.8103    OK

As far as I can tell the gene only exists on chr 10 so I am unsure what it is mapping to.  Is this an issue in the file I am inputting or some issue with the parameters? I am aligning to the UCSC, group: Gene and Gene prediction, track: RefSeq Genes, outpuf: GTF.

Thank you,

Matt

cufflinks galaxy • 824 views
ADD COMMENT • link • written 3.5 years ago by mgtrax10 • 0
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: 172 users visited in the last hour