Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • TopHat:Splice Junctions

    Hello Everyone,

    I am trying to use TopHat's GFF feature so that genes are normalized with RPKM. However, I cannot seem to get that aspect of the program working. If I include this -G/GFF option then TopHat will not report junctions.bed, or islands.gff or islands.bed. I am using test data to make sure everything is working properly before I invest the time in running real solexa runs. If anyone has any suggestions as to how to get GFF function in tophat it would be greatly appreciated.

    Here is Sample Reads:
    @SRR013411.1874 :1:1:838:688 length=50
    aaaaagtcacagtctgtctgtctgtctctctcctctaatcttttttatcc
    +SRR013411.1874 :1:1:838:688 length=50
    I>AI8I6II1BI,EIB4+0;51-.7'*,#(+&((*#'&&#"++&"!!*%$
    @SRR013411.1875 :1:1:83:374 length=50
    gctagggttttgaagcaaggtttctcgcgattttctcgatatctctcgcc
    +SRR013411.1875 :1:1:83:374 length=50
    IIIGDGIIIII@DIIII3I/ICI:412B;=:5/+2/&503-39&2#1D,(
    @SRR013411.1876 :1:1:44:341 length=50
    cATGGCGAAACCAAGTCGTGGCCGTCGTTCCCCCTCCGTGTCTGGCTCGT
    +SRR013411.1876 :1:1:44:341 length=50
    22*4$'3%+1#)41&2%2$$%)#8'%#$&%(+%$'&!,"+!$*&,&!%%&
    @SRR013411.1877 :1:1:557:896 length=50
    CATCTCGTTCCAGTTCCAGATCTCGTTCGGGTTCGAGCCCCTCCAGGTCT
    +SRR013411.1877 :1:1:557:896 length=50
    I1II+GI8CI86I%A-I89;=40+-)A0.%)+.&&)(-&0%$!"!$#("%
    @SRR013411.1878 :1:1:866:720 length=50
    ATTTCCCGCTCACGCTCCCGTTCTAGATCGCTCTCTTCATCTTCATCTCC
    +SRR013411.1878 :1:1:866:720 length=50
    8"(-,<=(-+&8'&%*&-')"$)%"&"#%!!$)#$!!$!"""!!"%!!""

    Here is Sample GFF:
    chr1 TAIR8 chromosome 1 30432563 . . . ID=chr1;Name=chr1
    chr1 TAIR8 gene 3631 5899 . + . ID=AT1G01010;Note=protein_coding_gene;Name=AT1G01010
    chr1 TAIR8 mRNA 3631 5899 . + . ID=AT1G01010.1;Parent=AT1G01010;Name=AT1G01010.1;Index=1
    chr1 TAIR8 protein 3760 5630 . + . ID=AT1G01010.1-Protein;Name=AT1G01010.1;Derives_from=AT1G01010.1
    chr1 TAIR8 exon 3631 3913 . + . Parent=AT1G01010.1
    chr1 TAIR8 five_prime_UTR 3631 3759 . + . Parent=AT1G01010.1
    chr1 TAIR8 cDS 3760 3913 . + 0 Parent=AT1G01010.1,AT1G01010.1-Protein;
    chr1 TAIR8 exon 3996 4276 . + . Parent=AT1G01010.1
    chr1 TAIR8 cDS 3996 4276 . + 2 Parent=AT1G01010.1,AT1G01010.1-Protein;
    chr1 TAIR8 exon 4486 4605 . + . Parent=AT1G01010.1
    chr1 TAIR8 cDS 4486 4605 . + 0 Parent=AT1G01010.1,AT1G01010.1-Protein;
    chr1 TAIR8 exon 4706 5095 . + . Parent=AT1G01010.1
    chr1 TAIR8 cDS 4706 5095 . + 0 Parent=AT1G01010.1,AT1G01010.1-Protein;
    chr1 TAIR8 exon 5174 5326 . + . Parent=AT1G01010.1
    chr1 TAIR8 cDS 5174 5326 . + 0 Parent=AT1G01010.1,AT1G01010.1-Protein;
    chr1 TAIR8 exon 5439 5899 . + . Parent=AT1G01010.1
    chr1 TAIR8 cDS 5439 5630 . + 0 Parent=AT1G01010.1,AT1G01010.1-Protein;
    chr1 TAIR8 three_prime_UTR 5631 5899 . + . Parent=AT1G01010.1

    Thanks for any help or suggestions.

    cheers

  • #2
    Can you send me the logs directory from the run? Also, what Bowtie index are you using, and have you verified that the chromosome names in the GFF match the chromosome names in the index?

    Comment


    • #3
      I have been struggling with a similar problem.

      If you are getting "malformed exon" errors in the report log, try adding ID fields to the attributes of "exon" lines in the GFF file; that worked for me.

      Good luck!

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Choosing Between NGS and qPCR
        by seqadmin



        Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
        10-18-2024, 07:11 AM
      • seqadmin
        Non-Coding RNA Research and Technologies
        by seqadmin




        Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

        Nobel Prize for MicroRNA Discovery
        This week,...
        10-07-2024, 08:07 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 06:09 AM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 10-30-2024, 05:31 AM
      0 responses
      13 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 10-24-2024, 06:58 AM
      0 responses
      22 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 10-23-2024, 08:43 AM
      0 responses
      52 views
      0 likes
      Last Post seqadmin  
      Working...
      X