This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Generating Counts Data From Fastq Sequence Files

I'm new to sequencing and I've currently got several FASTQ files containing data corresponding to sequencing experiments (sequenced using Illumina miseq).

I was hoping to carry out some expression analysis (with edgeR, probably) using this data, but I'll need to generate a counts matrix from this data. Could somebody provide some instruction on how to generate counts data from a FASTQ file?

sequencing fastq rna-seq counts differential-expression
• 11,881 views • link
updated by Xingyu Yang Xingyu Yang • written by josph.sh • josph.sh

2 answers

  1. You will have to first align those fastq files against the reference genome and produce SAM/BAM files.Tophat, STAR and many other splice aware RNA-seq aligners are available for this task. It is always good to preprocess your read data including QC, trimming off the low quality bases etc.

  2. Then you need to use some tool that will generate count data for you. Basically you will have to provide the aligned BAM file and the gene annotation file (gff3, gtf,bed format) for your reference genome. HTSeq, Cufflinks are some tools available for this task. Search "Biostar" and you will get names of other tools.

• 99 views • link
updated by Ram Ram • written by Ashutosh Pandey Ashutosh Pandey

http://www.nature.com/nprot/journal/v7/n3/abs/nprot.2012.016.html

• 434 views • link
written by Xingyu Yang Xingyu Yang

Log in to answer this question.

More posts like this