This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Rna-Seq Analysis Of 454/Roche Data

hi everyone i am having roche/454 transcripts reads in two different condition. i want to check the gene expression. whether it is good to use TopHat and cufflinks for 454/roche data.or suggest me any other software.

rna-seq tophat2 cufflinks
• 4,016 views • link
updated by JC JC • written by raj.gzra • raj.gzra

How many reads do you have per sample?

• 0 views • link
written by Damian Kao Damian Kao

2 answers

I am not sure if TopHat -> Cufflinks is a good choice since 454 does not generate enough throughput as Illumina and SOLiD. You need to estimate the expression using Tag count methods (such as HTSeq followed by DGESeq or DE-Seq or some other simple statistical test (Fisher's exact test??). There have been several papers using EST data to estimate the expression. You can follow their protocols. Try this: http://www.ncbi.nlm.nih.gov/pubmed/20043850

• 0 views • link
written by vinay kumar • vinay kumar

I will use Blat/BWA-Mem to align the reads, TopHat is not designed to map 454 sequences (long and variable reads, errors, ...), then you can use Cufflinks or HTSeq+EdgeR/DESeq for DEG.

• 0 views • link
written by JC JC

i mapped transcripts reads to already assembled transcriptome with BWA and when samtools flagstat command used, it shows number of input reads is higher than that i have. when i did it with bowtie and then samtools flagstat command used, it shows number of input reads equal to that i have. why samtools flagstat command shows higher number of input reads when reads mapped with BWA???

• 0 views • link
written by raj.gzra • raj.gzra

Log in to answer this question.

More posts like this