Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • coldturkey
    Member
    • Nov 2008
    • 51

    #1

    newbler ace to amos bnk

    does anyone know how to convert a newbler derived ace file to amos bnk so i can use amosvalidate. I know the standard command to do this using toAmos, however, we i do this lose mate pair information.
    I tried using toAmos with a traceinfo file derived from sff_extract which seemed to work but when i ran amosvalidate i encountered the following error

    Command: /usr/local/bioinf/Amos/bin/asmQC -b 454Contigs.bnk -scaff -recompute -update -numsd 2 exited with status:1


    can anyone help?

    regards

    Brian
  • coldturkey
    Member
    • Nov 2008
    • 51

    #2
    in case anyone is interested



    Unfortunately, we don't support 454 mates at this time, so I recommend you skip
    the mate analysis of amosvalidate. Copy the amosvalidate script to a new
    location:

    $ cp /path/to/amos/bin/amosvalidate amosvalidate-nomates

    Edit the script to remove the mate analysis by commenting out steps 300-340

    ## Analyzing mate-pairs
    ## 300: $(BINDIR)/asmQC -b $(PREF).bnk -scaff -recompute -update -numsd 2
    ## 310: $(BINDIR)/asmQC -b $(PREF).bnk -scaff -recompute -update -numsd 2
    ## 320: $(BINDIR)/asmQC -b $(PREF).bnk -scaff -feat -numsd 3 -shortcvg -1 -longcvg
    -1
    ## 330: $(BINDIR)/cestat-cov -i -f 4 $(PREF).bnk > $(PREF).ce.feat
    ## 340: $(BINDIR)/loadFeatures -i $(PREF).bnk $(PREF).ce.feat

    Then try the new script
    $ amosvalidate-nomates 454Contigs

    Good luck!

    Michael Schatz


    Which leads me to ask if there is any program which can be used to validate 454 mate pair assemblies

    Comment

    • AppleInformatics
      Junior Member
      • Nov 2009
      • 5

      #3
      Hi All,
      I am attempting to use amosalidate to perform assembly forensics. I am stuck, however, as everytime I follow the instructions on the amos wiki to convert my newbler .ace file to a .bnk, I get the following error:

      $ toAmos -ace 454Contigs.ace.1 -o - | bank-transact -m - -o 454Contigs.bnk -c
      bank-transact: invalid option -- o
      ERROR: The -b option is mandatory
      ERROR: Bank directory is not accessible, No such file or directory


      .USAGE.
      bank-transact [options] -b <bank path> -m <message path>
      Try 'bank-transact -h' for more information.




      Any suggestions? Thanks a bunch!

      Comment

      • shahid.manzoor
        Junior Member
        • Jun 2009
        • 8

        #4
        hi me shahid plz can any body help me i m using velvet for denovo and producing afg file and from that afg file creating amos bank file but when i run amos validate command on that bnk folder it give error like this:
        Doing step 100: Creating FastA
        Command: /usr/local/amos/bin/bank2fasta -b /home/shahid/denovo_run2/bank/assembly.bnk/.bnk > /home/shahid/denovo_run2/bank/assembly.bnk/.fasta exited with status: 1

        so any body can help me what this error come and any suggestion for solution.thanks

        Comment

        • ZhigangLi
          Member
          • Nov 2010
          • 11

          #5
          Originally posted by shahid.manzoor View Post
          hi me shahid plz can any body help me i m using velvet for denovo and producing afg file and from that afg file creating amos bank file but when i run amos validate command on that bnk folder it give error like this:
          Doing step 100: Creating FastA
          Command: /usr/local/amos/bin/bank2fasta -b /home/shahid/denovo_run2/bank/assembly.bnk/.bnk > /home/shahid/denovo_run2/bank/assembly.bnk/.fasta exited with status: 1

          so any body can help me what this error come and any suggestion for solution.thanks
          A similar error is returned when i use amosvalidate!
          github:
          https://github.com/Bioinformatics-and-Genomics

          Comment

          • Irsan_Kooi
            Member
            • Mar 2011
            • 10

            #6
            I have the same problem exactly the same problem as shahid.manzoor. I do velvet assembly with create .afg file option, put it in the AMOS bank, no problem. When I do amosvalidate I get the same response unfortunately...

            Does anyone know how to do it the right way?

            Comment

            • nathanhaigh
              Member
              • Jul 2011
              • 11

              #7
              Originally posted by AppleInformatics View Post
              I am attempting to use amosalidate to perform assembly forensics. I am stuck, however, as everytime I follow the instructions on the amos wiki to convert my newbler .ace file to a .bnk, I get the following error:

              [I]$ toAmos -ace 454Contigs.ace.1 -o - | bank-transact -m - -o 454Contigs.bnk -c
              bank-transact: invalid option -- o
              ERROR: The -b option is mandatory
              ERROR: Bank directory is not accessible, No such file or directory
              I know this is an old thread, but thought I would add what I've done recently to convert Newbler ACE files to AMOS bank:

              Code:
              toAmos -ace newbler_assembly.ace -o - | bank-transact -m - -b newbler_assembly.bnk -c

              Comment

              • nangillala
                Junior Member
                • Aug 2011
                • 9

                #8
                Originally posted by Irsan_Kooi View Post
                I have the same problem exactly the same problem as shahid.manzoor. I do velvet assembly with create .afg file option, put it in the AMOS bank, no problem. When I do amosvalidate I get the same response unfortunately...

                Does anyone know how to do it the right way?
                Hi,
                I know it's an old thread but I ran into the same problem. It seems as if the solution is the following:
                I guess you used tab-complete to write the bank name?
                So you would have something like this as an input:
                Code:
                amosvalidate assembly.bnk/
                just cut off the "/". So:
                Code:
                amosvalidate assembly.bnk
                seems to work for me.

                Maybe this helps someone running into similar problems.

                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, 10-22-2026, 11:13 AM
                0 responses
                28 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 10-22-2026, 10:48 AM
                0 responses
                20 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 10-15-2026, 10:52 AM
                0 responses
                38 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 10-14-2026, 05:39 AM
                0 responses
                44 views
                0 reactions
                Last Post SEQadmin2  
                Working...