Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SOAP: short oligonucleotide alignment program.

    New paper out describing a new alignment tool. I don't have a subscription to Bioinformatics, so I must be content with the abstract until I can dig up the paper...

    MOTIVATION: We have developed a program SOAP for efficient gapped and ungapped alignment of short oligonucleotides onto reference sequences. The program is designed to handle the huge amounts of short reads generated by parallel sequencing using the new generation Illumina-Solexa sequencing technology. SOAP is compatible with numerous applications, including single-read or pair-end resequencing, small RNA discovery, and mRNA tag sequence mapping. SOAP is a command-driven program, which supports multithreaded parallel computing, and has a batch module for multiple query sets. AVAILABILITY: http://soap.genomics.org.cn CONTACT: [email protected].
    http://www.ncbi.nlm.nih.gov/pubmed/18227114

    edit: Due to intermittent downtime of the original SOAP site above, I've hosted the current tar file as an attachment here. If the author has any problem with this please let me know.
    Attached Files

  • #2
    cant find SOAP site in canada

    this program is supposed to be available at



    but cant access the site

    i was handed a file for a resequencing project for a microcystis strain that doesnt produce toxins ... i dont have the solexa software ... i ran a program called rmap using a fasta version i created of the solexa file against the published m.aer 6megbase genome and the sequences look promising ... so how do people who dont have the solexa software align their sequences? i have 5,130,912 sequences of 36 bases each ... rmap identifies about 120,000 of them as having nomismatches with the published genome ... about 360,000 have 1 or none mismatches ... the rmap output has start and stop positions and an indicator of the strand ...

    i am thinking it would be easier to do this with a program like SOAP if it were publicly available ... ?

    any suggestions ?

    thanks a bunch in advance

    rudy

    Comment


    • #3
      hey rudy. welcome. definitely sounds like you need to take a look at this thread:

      Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc


      ...and the SOAP site is down for me too... :/

      Comment


      • #4
        yuchang ::AT:: genomics.org.cn

        ECO

        kind fellow in canada sent me a tar file with the SOAP software with apologies for inconvenience accessing site ... will let you know what i find ... thanks for list of links
        rudy

        Comment


        • #5
          SOAP site down, still

          It is disappointing when a paper comes out and the download site for the software doesn't even last a month. If anyone on this board has the SOAP package, could they post a mirror (or send it to me)?

          Thanks,

          Ryan

          Comment


          • #6
            Originally posted by myrna View Post
            It is disappointing when a paper comes out and the download site for the software doesn't even last a month. If anyone on this board has the SOAP package, could they post a mirror (or send it to me)?

            Thanks,

            Ryan
            As usual, if someone PM's me, I will host it!

            Comment


            • #7
              SOAP site is back up

              The SOAP site is available again

              Tim

              Comment


              • #8
                It's back up for me too, but I've hosted the tar in the first post just in case! Thanks RudyS.

                Comment


                • #9
                  Any hints on compiling this thing on OSX?

                  Code:
                  eco@##:~/perl/soap/soap_1.10$ sudo make
                  Password:
                  g++ -static -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c align.cpp -o align.chr.o
                  /var/tmp//ccrPoVlo.s:20045:FATAL:incompatible feature used: section type non_lazy_symbol_pointers (must specify "-dynamic" to be used)
                  make: *** [align.chr.o] Error 1
                  eco@##:~/perl/soap/soap_1.10$

                  Comment


                  • #10
                    Originally posted by ECO View Post
                    Any hints on compiling this thing on OSX?

                    Code:
                    eco@##:~/perl/soap/soap_1.10$ sudo make
                    Password:
                    g++ -static -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c align.cpp -o align.chr.o
                    /var/tmp//ccrPoVlo.s:20045:FATAL:incompatible feature used: section type non_lazy_symbol_pointers (must specify "-dynamic" to be used)
                    make: *** [align.chr.o] Error 1
                    eco@##:~/perl/soap/soap_1.10$
                    You can remove "-static" from the FLAGS= line in the makefile

                    Comment


                    • #11
                      Originally posted by jjcook View Post
                      You can remove "-static" from the FLAGS= line in the makefile
                      Thanks JJcook...that worked a treat. I have another problem with this one, this time on openSUSE11:

                      Code:
                      ericbox:/home/eco/dl/soap_1.10 # make
                      g++ -static -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c main.cpp -o main.chr.o
                      In file included from main.cpp:11:
                      align.h: In member function ‘void SingleAlign::Reverse_Seq()’:
                      align.h:286: error: ‘reverse’ was not declared in this scope
                      align.h: In member function ‘void SingleAlign::Reverse_Qual()’:
                      align.h:294: error: ‘reverse’ was not declared in this scope
                      main.cpp: In function ‘void usage()’:
                      main.cpp:294: error: ‘exit’ was not declared in this scope
                      main.cpp: In function ‘int mGetOptions(int, char**)’:
                      main.cpp:308: error: ‘atoi’ was not declared in this scope
                      main.cpp: In function ‘void RunProcess()’:
                      main.cpp:348: error: ‘exit’ was not declared in this scope
                      main.cpp:353: error: ‘exit’ was not declared in this scope
                      main.cpp:358: error: ‘exit’ was not declared in this scope
                      main.cpp:363: error: ‘exit’ was not declared in this scope
                      main.cpp:385: error: ‘exit’ was not declared in this scope
                      main.cpp:391: error: ‘exit’ was not declared in this scope
                      main.cpp:398: error: ‘exit’ was not declared in this scope
                      main.cpp: In function ‘int main(int, char**)’:
                      main.cpp:428: error: ‘exit’ was not declared in this scope
                      make: *** [main.chr.o] Error 1
                      Last edited by ECO; 07-15-2008, 09:05 PM. Reason: to follow up on the OSX change working...

                      Comment


                      • #12
                        In align.h, you may need to change the 'reverse' call to 'std::reverse'

                        In main.cpp, you should add
                        Code:
                        #include <stdlib.h>
                        to the top of the file.

                        Given I can't duplicate your error with my version of gcc, I can't test if those fixes will work -- but hopefully they do

                        Comment


                        • #13
                          Originally posted by jjcook View Post
                          In align.h, you may need to change the 'reverse' call to 'std::reverse'

                          In main.cpp, you should add
                          Code:
                          #include <stdlib.h>
                          to the top of the file.

                          Given I can't duplicate your error with my version of gcc, I can't test if those fixes will work -- but hopefully they do
                          We're gettin' there...thanks again for your help. The #include removed most of the errors:

                          Code:
                          g++ -static -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c main.cpp -o main.chr.o
                          In file included from main.cpp:12:
                          align.h: In member function ‘void SingleAlign::Reverse_Seq()’:
                          align.h:286: error: ‘reverse’ was not declared in this scope
                          align.h: In member function ‘void SingleAlign::Reverse_Qual()’:
                          align.h:294: error: ‘reverse’ was not declared in this scope
                          make: *** [main.chr.o] Error 1
                          ...then changing the align.h didn't appear to work...

                          Code:
                          g++ -static -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c main.cpp -o main.chr.o
                          In file included from main.cpp:12:
                          align.h: In member function ‘void SingleAlign::Reverse_Seq()’:
                          align.h:286: error: ‘reverse’ is not a member of ‘std’
                          align.h: In member function ‘void SingleAlign::Reverse_Qual()’:
                          align.h:294: error: ‘reverse’ is not a member of ‘std’
                          make: *** [main.chr.o] Error 1

                          Comment


                          • #14
                            Ah -- you need to add this to align.h:

                            Code:
                            #include<algorithm>

                            Comment


                            • #15
                              Heh...error keeps getting shorter....I also removed the -static flag on advice received from the SOAP team.

                              Code:
                              g++ -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c main.cpp -o main.chr.o
                              g++ -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c pairs.cpp -o pairs.chr.o
                              g++ -DMAXGAP=3 -DMAXHITS=10000 -DTHREAD -O3   -DDB_CHR -DREAD_60  -c param.cpp -o param.chr.o
                              param.cpp: In member function ‘void Param::SetMrnaTag(int)’:
                              param.cpp:98: error: ‘exit’ was not declared in this scope

                              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