Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Tophat different outputs according to the different number of threads

    Hi! All,

    I don't know anyone who had the same issue like I had, while using tophat software. I have a server which has 16 cores.

    I tried to run tophat with different thread number for one of my sample. Unfortunately I had different outputs. (single-end RNA-seq data)

    1. tophat -p 16 --no-coverage-search --solexa-quals --keep-tmp --b2-very-sensitive -G Gene_GTF --no-novel-juncs -o Tophat_out bowtie_index 1st_sample.fastq

    2. tophat -p 4 --no-coverage-search --solexa-quals --keep-tmp --b2-very-sensitive -G Gene_GTF --no-novel-juncs -o Tophat_out bowtie_index 1st_sample.fastq

    Any idea?
    Thanks in advance...

  • #2
    When you say you got different outputs do you mean that you got the same hits, but in a different order (which wouldn't be surprising), or that you got different hits (which would)?

    Comment


    • #3
      simonandrews, thanks for your inquiry

      Of course, I know about different order of outputs.
      In my case, the result hits were different. Simply, the read numbers in accepted_hits.bam and unmapped.bam were slightly different. (-p 4 vs -p 16)

      I don't know why this problem was happened. I tried the other number of threads (like, 1, 2, and 4) and their results were all same, except 16 (max num of available threads).

      --Kumtl

      Comment


      • #4
        If this is reproducible then I think it's one to send to the tophat authors.

        Comment


        • #5
          That's not bad idea. I'll let them know, if this problem happened again.
          Thanks

          Comment


          • #6
            Simon, I have another question about tophat.

            As you know from my post, I had different outputs because of different number of threads.
            So, I tried to remap reads from unmapped.bam, and I got another hits from them. Of course, I applied exactly same option on tophat command line.

            Do you have any idea?

            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