Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • counting junction reads in TopHat

    Hi,
    I am running TopHat and using the accepted_hits.bam file to generate counts for genes.
    How does TopHat identify junctions in this bam file, as exons?
    How are users generating read counts for genes which include the junction and exon reads?
    Thanks,
    Lana

  • #2
    If you just want to use the junctions directly, then, in addition to bam file, you should also have a junctions, insertions and deletions .bed file.

    If you are wanting to understand, then, you should look at Sam Format Specification. Your bam file can be viewed with samtools. If necessary, you can convert to sam with picard tools.
    In RNA-Seq data, basically, when tophat finds a read that splices across previously identified exon regions, that is, for one read R (80bp), say, of gene G (with 5 exons, say), between E2 and E3 (the intron between E2 and E3 is 500bp, say); lets say R1 = 30bp maps to E2 and R2 = 50bp maps to E3, then tophat writes this as 30M500N50M. This is a CIGAR string (from the sam format specification). In addition to that there is a "start" position in the SAM file format and using this you can find out the junction (check out other possible options for CIGAR string). Is this what you asked for?

    Comment


    • #3
      exon join frequencies from RNAseq data

      Hi,

      Don't know whether this will help - or whether I went about it the correct way - but this is my related experience.

      I was faced with a similar problem - determining the read evidence for particular exon joins within a certain genomic region.

      Basically what I did was take the RNAseq read data and use tophat to align to a particular genomic region for which I knew the exon positions.

      After getting a sam/bam file, I wrote a python script to parse the tophat match data i.e. the cigar string info: 30M500N50M to determine the exon join read frequencies.

      Shoot me any questions you like.

      Good luck,

      Alex

      Comment


      • #4
        After going thro' your question again, the simple picture as far as I understand is this: tophat first uses bowtie to align reads against the reference genome. Bowtie aligns reads without identifying splice sites. It tries to map entire read directly to a matching region of your genome (and it is fast). From here, tophat "knows" potential exons and it tries to map the other unmapped reads by splitting them in to smaller parts and finding a region.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Latest Developments in Precision Medicine
          by seqadmin



          Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

          Somatic Genomics
          “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
          05-24-2024, 01:16 PM
        • seqadmin
          Recent Advances in Sequencing Analysis Tools
          by seqadmin


          The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
          05-06-2024, 07:48 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 01:32 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-24-2024, 07:15 AM
        0 responses
        199 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-23-2024, 10:28 AM
        0 responses
        221 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-23-2024, 07:35 AM
        0 responses
        232 views
        0 likes
        Last Post seqadmin  
        Working...
        X