Announcement

Collapse

Welcome to the New Seqanswers!

Welcome to the new Seqanswers! We'd love your feedback, please post any you have to this topic: New Seqanswers Feedback.
See more
See less

question about tophat

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • question about tophat

    I am running tophat index/mm9 1.fastq

    to get accepted_hist.bam

    The running process has no error. But 1.fastq has 1Gb, accepted_hist.bam only has 470K.

    So What does it mean? It means only a few reads has mapped to reference genome?

    I downloaded mm9 from bowtie website and contructed a .fa file using bowtie-inspect mm9 > mm9.fa

    Anyone know this problem?

    thanks!

  • #2
    I am curious why no people reply this? anyone can help?

    Comment


    • #3
      Did you check look at your alignment file to see if anything was aligned or check the tophat log files?

      Comment


      • #4
        thanks, I have checked my result , I really got some result but the result number is so small.

        Comment


        • #5
          Originally posted by chadn737 View Post
          Did you check look at your alignment file to see if anything was aligned or check the tophat log files?
          Hi chadn737,

          May you recommend some good method to extract gene expression from RNA seq data.

          Now I am using tophat and cufflinks but the result seems not I want.

          BTW, does the reference genome mm9 I used is right?

          Comment


          • #6
            Originally posted by camelbbs View Post
            I am running tophat index/mm9 1.fastq

            to get accepted_hist.bam

            The running process has no error. But 1.fastq has 1Gb, accepted_hist.bam only has 470K.

            So What does it mean? It means only a few reads has mapped to reference genome?

            I downloaded mm9 from bowtie website and contructed a .fa file using bowtie-inspect mm9 > mm9.fa

            Anyone know this problem?

            thanks!
            This is absolutely normal. accepted_hits.bam is in the binary BAM format which uses much less space than fastq (which is in text format).

            Comment


            • #7
              Originally posted by shurjo
              This is absolutely normal. accepted_hits.bam is in the binary BAM format which uses much less space than fastq (which is in text format).
              True, but 470K is a very small bam file, even for a small data set of only 1GB.

              Comment


              • #8
                Originally posted by chadn737 View Post
                True, but 470K is a very small bam file, even for a small data set of only 1GB.
                Agreed. Maybe we will get a better insight after camelbbs looks at the actual BAM file to see what's in it.

                Another option would be to run FastQC on the input file to check the quality of the reads.

                Comment

                Working...
                X