SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   RNA Sequencing (http://seqanswers.com/forums/forumdisplay.php?f=26)
-   -   RNA seq analysis - small number of genes (http://seqanswers.com/forums/showthread.php?t=63169)

sparky 10-02-2015 12:05 AM

RNA seq analysis - small number of genes
 
Hi all,

I am new to RNA seq analysis and want to look at the expression of a small number of genes (8) in some publicly available RNA seq datasets.

I came up with a simple method that avoids me having to align the RNA seq reads to the genome and do a full tophat/cufflinks analysis (or similar). Briefly, what I did was: Download and QC filter SRA dataset > map the reads to a multi-fasta file containing the exonic sequences for my genes of interest with bowtie2 > filter out hits with MAPQ<30 > obtain FPKM values using eXpress.

Does anyone have any thoughts on whether my method is valid or not? My main concern is that mapping reads to such a low complexity reference may artificially inflate the number of reads that map to my genes of interest and bias the FPKM values?

Thanks in advance :)


All times are GMT -8. The time now is 08:56 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.