Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • KevinLam
    Senior Member
    • Nov 2009
    • 204

    SOAP alignment format convert to SAM/BAM

    Dear All,
    would you know if there's a convertor for SOAP alignment format to SAM?
    http://kevin-gattaca.blogspot.com/
  • gaffa
    Member
    • Oct 2010
    • 82

    #2
    There is a script called "soap2sam.pl" - find it here: http://soap.genomics.org.cn/soapaligner.html

    Comment

    • Aman Mahajan
      Member
      • Jan 2012
      • 22

      #3
      I assembled my illumina data using SOAP, now I want to carry out expression analysis using Rseq tool. it accepts only SAM format so I downloaded SAMTOOLS to convert my soap output to SAM. Can anyone tell me how to run it and convert, tutorial has been of no use yet!

      or if you can throw light on how to use this script.

      Comment

      • gringer
        David Eccles (gringer)
        • May 2011
        • 845

        #4
        From the source code from that script:

        Code:
        Usage: soap2sam.pl [-p] <aln.soap>
        This will output a SAM format file on standard output. If you want that output to be sent to a particular file, then put a redirection sequence ('> aln_file.sam') at the end of the command line.

        It looks like the '-p' option is used to indicate that your alignments are paired-end.

        Comment

        • Aman Mahajan
          Member
          • Jan 2012
          • 22

          #5
          Thanks. but I actually want to know how can I input my file to run through this.

          Comment

          • gringer
            David Eccles (gringer)
            • May 2011
            • 845

            #6
            Have you aligned your reads using SOAPaligner? That program produces output that can be fed into soap2sam as the 'aln.soap' command-line argument.

            Comment

            • Aman Mahajan
              Member
              • Jan 2012
              • 22

              #7
              My alignment output has generated a lot of files.. which file out of these has to be the input for soap2sam.pl.. Should the .contig file be renamed as .soap for input?

              Comment

              • ganygan25
                Junior Member
                • Dec 2011
                • 9

                #8
                While running soap, the output file name can be explicitly specified.
                soap -D <Reference Genome Index> <Read Files> -o <Output File>
                (Note:- no file extension required)

                This output File can then be converted to SAM format as:-
                soap2sam.pl -p <Output File> > <Sam File>.sam

                Comment

                • Aman Mahajan
                  Member
                  • Jan 2012
                  • 22

                  #9
                  @ganygany - my output file name is abc.soap.

                  now I am putting this command

                  soap2sam.pl abc.soap abc.sam

                  It generates a sam file but it gives an error message.

                  Comment

                  • gringer
                    David Eccles (gringer)
                    • May 2011
                    • 845

                    #10
                    you missed out the > symbol

                    Comment

                    • Aman Mahajan
                      Member
                      • Jan 2012
                      • 22

                      #11
                      I have put this symbol now

                      soap2sam.pl abc.soap > abc.sam..

                      Error of command not found. If i put this command :-

                      perl soap2sam.pl abc.soap > abc.sam..

                      Still the same result, file generated of the name abc.soap but zero space..

                      Comment

                      • ganygan25
                        Junior Member
                        • Dec 2011
                        • 9

                        #12
                        "Error of command not found. "
                        Not a soap issue. Could your verify your installation and run them again?
                        Or post your sequence of commands here..

                        "file generated of the name abc.soap but zero space"
                        If your commands are fine, do you mean to say that there are no reads aligned ??

                        Comment

                        • Aman Mahajan
                          Member
                          • Jan 2012
                          • 22

                          #13
                          Okay. I aligned my data using SOAP so I have a file abc.soap which is aligned.

                          I have copied this file in the same folder with soap2sam.pl. Now I am trying to run soap2sam.pl on this file abc.soap.

                          My command is not showing any error but the output generated i.e abc.sam is of zero space.

                          Any idea if there is a need to compile soap2sam.pl first because I am using it through samtools.

                          Thank you.

                          Comment

                          • gringer
                            David Eccles (gringer)
                            • May 2011
                            • 845

                            #14
                            could you please show the result of running the following command (including the exact error, if any)? I suspect that soap2sam isn't in the directory that you're trying to run it from.

                            Code:
                            ls $(pwd)/soap2sam.pl $(pwd)/abc.soap

                            Comment

                            • Aman Mahajan
                              Member
                              • Jan 2012
                              • 22

                              #15
                              result is


                              /root/Desktop/samtools-0.1.18/misc/soap2sam.pl


                              in the same folder i have my file abc.soap

                              Comment

                              Latest Articles

                              Collapse

                              • SEQadmin2
                                From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
                                by SEQadmin2


                                Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


                                The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
                                ...
                                06-02-2026, 10:05 AM
                              • SEQadmin2
                                Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends
                                by SEQadmin2


                                With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.


                                Introduction

                                Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
                                05-22-2026, 06:42 AM
                              • SEQadmin2
                                Environmental Genomics in the Age of NGS: From Microbes to Conservation Strategies
                                by SEQadmin2

                                Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.


                                Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...
                                05-06-2026, 09:04 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by SEQadmin2, Yesterday, 08:59 AM
                              0 responses
                              14 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 06-02-2026, 12:03 PM
                              0 responses
                              22 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 06-02-2026, 11:40 AM
                              0 responses
                              19 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 05-28-2026, 11:40 AM
                              0 responses
                              32 views
                              0 reactions
                              Last Post SEQadmin2  
                              Working...