Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Input file "intree"

    Hi Bas,

    I just installed the stand alone version on my mac and have been trying to run it.
    It seems that drawtree is looking for the input tree file "intree" which was not made.
    Is there some configuration settings I missing for drawtree to allow it to accept cladogram.ph files?

    Many thanks,
    Martin

    Comment


    • #17
      Originally posted by martintay View Post
      Hi Bas,

      I just installed the stand alone version on my mac and have been trying to run it.
      It seems that drawtree is looking for the input tree file "intree" which was not made.
      Is there some configuration settings I missing for drawtree to allow it to accept cladogram.ph files?

      Many thanks,
      Martin
      Hi Martin,
      The drawtree program is used to visualize the cladogram, taking a Phylip bracketnotation as input, and outputting a PNG image. It is called from within crAss (v1.2) with the following system call (with branch lengths):

      echo -e "output.$size_correction.cladogram.ph\n$crAss_dir/fontfile\nV\nN\nL\nA\nD\nM\n0\n0\\nF\nHelvetica\nC\n0.5\nY\n" | $drawtree_exe

      Something is probably going wrong with this system call. There are three variables, as indicated with the dollar sign (basic if you know some Perl...):
      - $size_correction is either "minimum", "reads", "shot", or "wootters"
      - $crAss_dir/fontfile is the location of the file "fontfile" that drawtree uses (it should be in the directory where the crAss.pl script is)
      - $drawtree_exe is the location of the drawtree executable in your shell

      By default, drawtree looks for a Phylip cladogram file called "infile", but as you can see above, the input file should be called "output.$size_correction.cladogram.ph". This means the location of the input file in the system call was not passed correctly, possibly as a result of different use of "echo" on a Mac.

      As a solution, you could try running the drawtree command on your command line. Alternatively, and this may be the easiest solution, is not to use drawtree but another program to visualize the cladogram files (the four "output.$size_correction.cladogram.ph" files in the output directory). You could use any tree plotting program for this, such as FigTree, NJPlot, or the website http://itol.embl.de/ are just a few examples. Those programs take the bracketnotation in the *.ph files as input and show you a cladogram (tree) as output. This is all that drawtree does as well.

      Good luck, let me know if/where you get stuck!
      Bas

      Comment


      • #18
        Hi Bas,

        Thanks for the advice, it is indeed the echo command that was giving the problem.
        I continued on manually and got what I needed.

        Cheers,
        Martin

        Comment


        • #19
          Hi,

          I would like to try crAss.pl out on a metagenome study.

          Two treatments
          Three time points
          Three technical replicates per time point

          I have a cross assembly of all the samples, though I do not have an ace file. The assembler that I used produces only a fasta file of assembled contigs.

          I'm wondering how to use crass.pl properly on this study. Any guidance would be greatly appreciated. Thanks.

          Comment


          • #20
            Hi Brettin, the new version of crAss can also take a SAM file as input. This allows you to first make cross-contigs, and then do the read mapping back to those contigs for example with Bowtie2, which creates a SAM file which you could use as input. You can find the latest crAss version on SourceForge (https://sourceforge.net/projects/crass/). Best! Bas


            Originally posted by brettin View Post
            Hi,

            I would like to try crAss.pl out on a metagenome study.

            Two treatments
            Three time points
            Three technical replicates per time point

            I have a cross assembly of all the samples, though I do not have an ace file. The assembler that I used produces only a fasta file of assembled contigs.

            I'm wondering how to use crass.pl properly on this study. Any guidance would be greatly appreciated. Thanks.

            Comment

            Latest Articles

            Collapse

            • 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...
              Yesterday, 07:48 AM
            • seqadmin
              Essential Discoveries and Tools in Epitranscriptomics
              by seqadmin




              The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
              04-22-2024, 07:01 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Today, 06:57 AM
            0 responses
            9 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 07:17 AM
            0 responses
            13 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 05-02-2024, 08:06 AM
            0 responses
            19 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-30-2024, 12:17 PM
            0 responses
            22 views
            0 likes
            Last Post seqadmin  
            Working...
            X