Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Maria_Lu
    replied
    Hi, Lifeng,

    Thank you so much for you kindness.
    I'll calculate these junctions again according to your remind.

    Leave a comment:


  • lifeng.tian
    replied
    Maria,

    Since you did not give the complete junctions.bed line, I attached two lines
    from my TopHat analysis:

    75nt run
    chr20 251862 256677 JUNC00000001 11 - 251862 256677 255,0,0 2 46,69 0,4746


    50nt run
    chr20 251879 257723 JUNC00000001 2 - 251879 257723 255,0,0 2 29,39 0,5805


    Here is the link to BED format at ucsc: http://genome.ucsc.edu/FAQ/FAQformat.html#format1

    You can calculate the (start,end) of the junction by:
    e.g.for 75nt run, find col 2, 10, 11, 12:
    col 2: feature start 251862
    col 10: block count 2
    col 11: block sizes 46,69
    col 12: block starts 0,4746
    You can find the splice junction coordinates by:
    start: col2+col11.first.item+1
    end: col2+col12.second.item,
    i.e.,
    start: 251862+46+1=251909, end: 251862+4746=256608
    chr20:251909-256608


    for 50nt run,
    junction is : chr20:251909-257684


    Now, in this case, TopHat finds two different junctions. You can also
    copy and paste the BED directly to ucsc custom track and visualize them.

    Hope it helps.

    Lifeng

    Leave a comment:


  • Maria_Lu
    started a topic TopHat: the results confused me

    TopHat: the results confused me

    I used TopHat to find exon-exon junctions.

    But the results in the output 'junctions.bed' files confused me.
    I seperated RNA-seq data into two datasets as one is 76bp*2 another is 40bp*2, then ran TopHat individually.

    However, the two 'junctions.bed' files gave different results.
    Here are examples of the two 'junctions.bed' result:
    One reports the following junctions:
    chromosome12 12302 12721 JUNC00000002 5 -
    chromosome12 33389 34997 JUNC00000003 6 +
    chromosome12 33688 34964 JUNC00000004 2 +
    chromosome12 35474 35675 JUNC00000005 5 +
    chromosome12 35718 35949 JUNC00000006 9 +

    Another reports the following junctions:
    chromosome12 12303 12723 JUNC00000005 26 -
    chromosome12 33679 34982 JUNC00000007 3 +
    chromosome12 35490 35674 JUNC00000008 6 +
    chromosome12 35711 35949 JUNC00000009 7 +

    These junction locations of each output file were similar but different.
    When I extracted the detailed sequences, no GT-AG was found.

    Does anybody know how to explain it?

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, 05-24-2024, 07:15 AM
0 responses
195 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-23-2024, 10:28 AM
0 responses
218 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-23-2024, 07:35 AM
0 responses
223 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-22-2024, 02:06 PM
0 responses
12 views
0 likes
Last Post seqadmin  
Working...
X