Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • visualizing scripture output BED

    Hi,

    I'm working through an analysis with Scripture for the first time.

    I've been able to generate a bed and dot file. UCSC does not seem like this BED file.

    I was able to find this suggestion that scripture includes a gff converter:



    java -Xmx1000m -jar <path_to_scripture>/scripture.jar -task toGFF -cufflinks -in scripture.t0.bed -source SCRIPTURE -out scripture.t0.gtf -prefix t0

    but... this toGFF function does not seem to be described in the scripture documentation, and does not seem to work for me either.

    Any suggestions? All I am really trying to do is to get this into UCSC so I can quickly assess how things look. Scripture suggests IGV browser which I suppose I could also use.

    thanks!

  • #2
    What happens when you try upload the .bed file to the UCSC browser?

    Comment


    • #3
      error message

      I get this:

      Code:
      File 'allScores.bed' - Unrecognized format line 1 of custom track: chr1	199	249	chr1:199-249	0.0	*	199	249	0,0,0	1	50,	0,	1.0	0.3684297339584846	1.0	0.02 3.1075660687968516	0.0542844351489113	50.0	1.0 (note: chrom names are case sensitive)
      thoughts? Maybe I need some sort of header?

      generally the scripture documentation seems somewhat half written, unless I am missing something?

      Comment


      • #4
        fyi, the chromosome name is fine... (no caps problem there)

        Comment


        • #5
          Hi,

          I have been able to open bed files in UCSC Genome browser using Scripture-beta2 version and this type of command:

          java -Xmx2000m -jar scripture-beta2.jar -task extractDot -in Segmentation_Results/chr1.scriptureESTest.segments.graph.post.add.pairs.dot -chr chr1 -start 6235241 -end 6239548 -out chr1_6235241-6239548.dot

          It gives a dot file and a bed file for the indicated region of the chromosome.
          Hope this help !

          Emilie
          Last edited by Emilie; 01-26-2011, 05:57 AM.

          Comment


          • #6
            .bed files don't require a header, but you can include one. Bed format instructions:
            http://genome.ucsc.edu/FAQ/FAQformat#format1
            I think you have too many columns in your file. Can you post the first few lines from your file, exactly as written in the file?

            Comment


            • #7
              Originally posted by chrisbala View Post
              Hi,

              I'm working through an analysis with Scripture for the first time.

              I've been able to generate a bed and dot file. UCSC does not seem like this BED file.

              I was able to find this suggestion that scripture includes a gff converter:



              java -Xmx1000m -jar <path_to_scripture>/scripture.jar -task toGFF -cufflinks -in scripture.t0.bed -source SCRIPTURE -out scripture.t0.gtf -prefix t0

              but... this toGFF function does not seem to be described in the scripture documentation, and does not seem to work for me either.

              Any suggestions? All I am really trying to do is to get this into UCSC so I can quickly assess how things look. Scripture suggests IGV browser which I suppose I could also use.

              thanks!
              Would you tell me that how can you transfer the bed file from scripture to gtf file please?

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 06:37 PM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              49 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X