Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Is anyone having trouble with tophat's SAM file -> cufflinks?

    I am having trouble running cufflinks on my own SAM file generated from tophat. I used many SAM files generated from real data as well as simulated data, but I can't get cufflinks to assemble the transcripts using any of my SAM files.

    I tried downloading the file: test_data.sam and it works fine. It identifies three exons and notes that they are part of one transcript. When I try this on my SAM file, it identifies all the exons as their own transcripts, even though I see many reads that span the exons.

    Any help would be appreciated. Thanks

  • #2
    Did you run Tophat on more than one sequence reads at the same time? If so you might want to run it one at a time. If not, then well, at the moment I can't think of an answer.

    Siva

    Comment


    • #3
      Would this matter if I specify a different ouuput file for each run?

      I am trying to figure out if the SAM format is a problem or if the actual reads do not cover the junction enough for the transcript to be called.

      Comment


      • #4
        Actually it is really strange. I cannot get any of my tophat-generated SAM outputs to successfully run using cufflinks (it runs but it lists each exon as its own transcript). I've done this with real data and simulated data.

        Comment


        • #5
          Originally posted by peterlchang View Post
          Actually it is really strange. I cannot get any of my tophat-generated SAM outputs to successfully run using cufflinks (it runs but it lists each exon as its own transcript). I've done this with real data and simulated data.
          Hi are u working with an organism whose typical intron size is very large? If so you can change Tophat's default parameters. Cufflinks typically gives you a GTF file detailing FPKM etc, a comma delimited transcript.expr file and a list of gene coordinates (genes.expr).

          Comment


          • #6
            The gene I am testing is 34K long (ie. the reference sequence is 34K long). The default intron length for tophat is 500K and I changed that parameter to 50K. The default intron length for tophat is 300K, and that was left unchanged.

            In addition, tophat is able to find reads that span the junction. I don't know what the issue could be. I already emailed Cole so I hope he can get back to me asap.

            Thanks

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Exploring the Dynamics of the Tumor Microenvironment
              by seqadmin




              The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
              07-08-2024, 03:19 PM
            • seqadmin
              Exploring Human Diversity Through Large-Scale Omics
              by seqadmin


              In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
              06-25-2024, 06:43 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 07-10-2024, 07:30 AM
            0 responses
            25 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 07-03-2024, 09:45 AM
            0 responses
            201 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 07-03-2024, 08:54 AM
            0 responses
            211 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 07-02-2024, 03:00 PM
            0 responses
            193 views
            0 likes
            Last Post seqadmin  
            Working...
            X