Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Problems with TopHat and Paired End RNASeq

    Hi, I'm very new at RNASeq Analysis and I think that I'm doing something wrong.
    I've tried to align that data: http://trace.ncbi.nlm.nih.gov/Traces...mple=SRS167019.

    First I execute sratoolkit to extract the fastq files with this command:

    fastq-dump --split-3 ./1_SRR097897.sra
    Then I have two files: 1_SRR097897_1.fastq and 1_SRR097897_2.fastq.

    The next step that I've done is to execute TopHat to generate the bam file with the reads aligned. The command was:

    ../tophat2 -N 0 Pombe 1_SRR097897_1.fastq 1_SRR097897_2.fastq
    Then I use the IGB to visualize the results and than shows the + and the - strands are the same



    If I only use one of the _n files to generate the output I have two strands with polarity, but depending on the file to select the polarity is reversed.



    I have also tried running TopHat with different parameters enabling --library-type or a GFF annotation file with identical results.
    I could stay only with the results obtained through _2.fastq file but I wonder to know what I'm doing wrong, for other occasions.

    Thanks for your time!!

  • #2
    I think it makes sense that the different paired ends will be on opposite strands (since you basically read the second end forwards, but on the other strand). At the level of zoom you're at in the included figures, I'm not sure that you could count on seeing the different spacing of the read1s vs the read2s (in my experiments, something like 100-300 bp). It doesn't seem like the best solution, but you could perhaps go through the BAM file and change all of the "SEQ being reverse complemented" bits in the FLAG field to get both ends pointing the same way.

    Comment


    • #3
      Thanks for your answer.
      Finaly I change the reads in one of the FASTQ files to their reverse complementary form and I can align them with polarity but I haven't understood yet why it works :S

      Comment

      Latest Articles

      Collapse

      • seqadmin
        New Genomics Tools and Methods Shared at AGBT 2025
        by seqadmin


        This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

        The Headliner
        The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
        03-03-2025, 01:39 PM
      • seqadmin
        Investigating the Gut Microbiome Through Diet and Spatial Biology
        by seqadmin




        The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
        02-24-2025, 06:31 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Today, 12:50 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-03-2025, 01:15 PM
      0 responses
      181 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-28-2025, 12:58 PM
      0 responses
      275 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-24-2025, 02:48 PM
      0 responses
      663 views
      0 likes
      Last Post seqadmin  
      Working...
      X