Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • anyone1985
    Member
    • Mar 2009
    • 68

    #1

    how convert ace to phd without consed?

    I have an ace file generated by amos2ace, but it can not be opened by consed. what i want is a phd file coming from the ace file. So can any tools convert the ace file to phd file?
  • TJK-OHSU
    Junior Member
    • Jan 2009
    • 9

    #2
    RE: how convert ace to phd without consed?

    Consed comes with several utility programs:
    ace2Fasta.perl generates a fasta file (actually I think it is fastq)
    then
    fasta2Phd.perl would give you the phd file you want.
    You need to get these from the good people at U. Washington though.
    Tom K
    OHSU

    Comment

    • anyone1985
      Member
      • Mar 2009
      • 68

      #3
      I got consed v19, and mapped the Solexa reads to my reference sequence. It work well. But I had more two hundred contigs and how can i convert the ace file to phd files seperately by consed/phrap at one time? I saw that every contig could consert to phd file by file->export consesus sequence with options.

      Comment

      • cabauser
        Junior Member
        • Aug 2008
        • 3

        #4
        Viewing Ace Files

        There are several open source software to view the alignments in the assembled contigs: Eagleview (Marth Lab, Boston College). More are listed in other posts. If you want to edit and manipulate the contigs try the commercial software Lasergene SeqMan from DNAStar. Both programs have trouble with large numbers of reads. I don't know how Consed works with 5 Mio reads (or more) aligned to a reference sequence.
        Last edited by cabauser; 03-22-2009, 02:01 AM. Reason: correct inaccuracy

        Comment

        • anyone1985
          Member
          • Mar 2009
          • 68

          #5
          I want to convert ace file which contains thousands contigs to phd file at one time by consed, is there any good way?

          Comment

          • Rao
            Member
            • Oct 2008
            • 36

            #6
            I've no idea about consed... but we can use BioPerl for this...
            see this link for more on this... http://www.bioperl.org/wiki/Module:Bio::SeqIO

            Comment

            • sklages
              Senior Member
              • May 2008
              • 628

              #7
              @anyone1985,

              ace format is an alignment format, phd isn't. What do you want to achieve?

              You have mapped your solexa reads against a reference sequence and get a few thousand contigs? Optimise your reference

              If you just want to reassemble your contigs, convert ace to fasta and use phrap and/or cap3 directly?

              Sven

              Comment

              • anyone1985
                Member
                • Mar 2009
                • 68

                #8
                I just assemble my solexa data with velvet, and get about 240 contigs larger than 500. Then, I map the solexa data to the contigs using the consed/addSolexaReads.perl to get a contigs with quality. Now, I want to convert the ace file to phd file individually. The consed can export the phd file only one at a time, I just want to know how to convert it all at a time.

                Comment

                • sklages
                  Senior Member
                  • May 2008
                  • 628

                  #9
                  Well, I think you need to write a script on your own, which parses the ACE files and dumps PHD. That's probably the fastest way. No other idea, sorry.

                  Sven

                  Comment

                  • anyone1985
                    Member
                    • Mar 2009
                    • 68

                    #10
                    thanks^^^^^^^^^^^^^^

                    Comment

                    • baohua100
                      Senior Member
                      • Jun 2008
                      • 103

                      #11
                      What's Consed's memory usage when open ACE file contains 10 million solexa reads ?

                      Who knowns the memory usage about Consed?

                      Comment

                      • anyone1985
                        Member
                        • Mar 2009
                        • 68

                        #12
                        To be honest, I maped my reference sequence with about 7 million reads, which is about half of my data. I divided the reference into six files, and the server that i use worked well. If i put them together, the program would be killed. My server is inspur with 8 G menory, about 4 G swap.

                        Comment

                        • sklages
                          Senior Member
                          • May 2008
                          • 628

                          #13
                          @baohua100,

                          opening an ACE file with 6.5Mio reads (2Mb genome) used 5.5GB physical memory.

                          Comment

                          • baohua100
                            Senior Member
                            • Jun 2008
                            • 103

                            #14
                            Originally posted by sklages View Post
                            @baohua100,

                            opening an ACE file with 6.5Mio reads (2Mb genome) used 5.5GB physical memory.
                            what's the version of consed ? 19.0 ?

                            what's the length of your read?

                            what's the cpu time load this ACE file?

                            thanks!
                            Last edited by baohua100; 03-25-2009, 03:28 AM.

                            Comment

                            • sklages
                              Senior Member
                              • May 2008
                              • 628

                              #15
                              Originally posted by baohua100 View Post
                              what's the version of consed ? 19.0 ?

                              what's the length of your read?

                              what's the cpu time load this ACE file?

                              thanks!
                              v19.00, 32b PE, 6min to open, didn't check cpu time

                              Comment

                              Latest Articles

                              Collapse

                              • SEQadmin2
                                Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                                by SEQadmin2



                                Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                                Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                                ...
                                10-14-2026, 09:55 AM
                              • SEQadmin2
                                New Genomics Technologies Take Aim at Long-Standing Limits
                                by SEQadmin2


                                Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

                                We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
                                ...
                                09-28-2026, 10:25 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by SEQadmin2, Yesterday, 11:13 AM
                              0 responses
                              21 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, Yesterday, 10:48 AM
                              0 responses
                              18 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 10-15-2026, 10:52 AM
                              0 responses
                              36 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 10-14-2026, 05:39 AM
                              0 responses
                              42 views
                              0 reactions
                              Last Post SEQadmin2  
                              Working...