SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Simple Bowtie question nilmot13 Bioinformatics 5 03-07-2012 02:03 AM
question about tophat camelbbs Bioinformatics 7 07-07-2011 06:24 PM
a question about tophat lei RNA Sequencing 12 03-20-2011 12:53 PM
One question about Tophat fgh1124 RNA Sequencing 0 12-15-2010 04:42 PM
A simple question on cDNA and protein sequence cheng Bioinformatics 1 10-06-2010 08:36 PM

Reply
 
Thread Tools
Old 06-02-2011, 06:03 PM   #1
kwatts59
Member
 
Location: nevada

Join Date: Apr 2011
Posts: 46
Default Simple Tophat question

I am using Tophat to perform single-end alignment to a plant genome that has 12 chromosomes. Tophat generated a file called accepted_hits.bam (1.75G). I converted the file to a bed file using the bamToBed utility.
Then I split the bed file based on chromosomes. Below are the number of lines (reads) for each chromosome.
chrom1 - 4.9M
chrom2 - 6.2M
chrom3 - 5.0M
chrom4 - 2.4M
chrom5 - 5.3M
chrom6 - 2.8M
chrom7 - 3.0M
chrom8 - 2.6M
chrom9 - 12.6M
chrom10 - 1.7M
chrom11 - 1.8M
chrom12 - 1.8M

My question: Why are there 12.6 million reads on chromosome 9?
I looked at the genome map and chromosome 9 is the second smallest (chrom 10 is the smallest).
kwatts59 is offline   Reply With Quote
Old 06-02-2011, 06:23 PM   #2
ishmael
Member
 
Location: NY, US

Join Date: Jul 2008
Posts: 17
Default

I think you should check the data first.
Perhaps some high abundant sequences take large proportion of the dataset.(For example, rRNA reducing failed, so some rRNA occupied many reads?)
Check the density distribution of the reads at the chromosomes, it may provide some clues.
ishmael is offline   Reply With Quote
Old 06-02-2011, 07:31 PM   #3
jmartin127
Member
 
Location: San Francisco

Join Date: Mar 2010
Posts: 15
Default

Welcome to RNA-Seq. You may want to first forget everything you learned about genome sequencing.

I would align the reads to the transcriptome FASTA file, assuming you have one, and check which genes are most highly expressed, and if they are on chromosome 9. The size of the chromosome may be irrelevant. It could be that chromosome 9 is much more gene dense than the other chromosomes. Ishmael's suggestion is also good.

Quote:
Originally Posted by kwatts59 View Post
My question: Why are there 12.6 million reads on chromosome 9?
I looked at the genome map and chromosome 9 is the second smallest (chrom 10 is the smallest).
jmartin127 is offline   Reply With Quote
Reply

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 01:48 PM.


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