Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Handling of singletons in Cufflinks jhb1980 Bioinformatics 0 10-28-2011 05:36 AM
Handling alignments spanning chromosomes boundaries in BWA nirav99 Bioinformatics 13 08-14-2011 12:48 AM
illumina single-end reads run cufflink louis7781x Bioinformatics 3 04-23-2011 06:05 AM
R error handling Stiff0810 Bioinformatics 0 03-24-2011 07:49 AM
Any good programs for handling solexa data? Burkholderia Introductions 1 08-13-2010 01:50 PM

Thread Tools
Old 11-17-2010, 04:02 AM   #1
Location: Denmark

Join Date: Mar 2010
Posts: 10
Default Cufflink's handling of reads

I am using Cufflinks to analyze RNAseq data and comparing different samples.
I run the programs with the following options:
cufflinks --num-threads 5 --GTF $ref_transcriptome --verbose --output-dir lane5/ $gene_name.lane5.bam >cufflinks.lane5.log 2>&1
cuffcompare -o summary.stats -r $ref_transcriptome -V lane5/transcripts.gtf lane6/transcripts.gtf
1) How does cufflink deal with reads only partially overlapping an exon (from a GTF file)?
2) How are reads administered between overlapping genes with overlapping exons (the RNA data are unstranded)?
3) If a read is only partially overlapping with the exon, will that part contribute to the FPKM value?
4) How is the minimum number of reads needed in a transcript to produce an FPKM value controlled? I can often visually se reads overlapping a transcript in my GTF file, but still giving scores of zero in the transcripts.expr
vang42 is offline   Reply With Quote

cufflink, rna seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:06 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO