Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Hi Boetsie,

    I've tried to run the SSPACE example (E. Coli, contigs_abyss.fasta, etc.) in order to test the program.
    I follow all the steps, but the final result is a file with 111 scaffolds. Why is that? I was supposed to get only a few scaffolds, isn't?

    Could you please give me some advices? I don't know what I'm doing wrong.

    Thank you very much.
    Germán.

    Comment


    • Originally posted by boetsie View Post
      yes, this is possible. The file should be in a TAB delimited format like:

      <contig1> <startpos_on_contig1> <endpos_on_contig1> <contig2> <startpos_on_contig2> <endpos_on_contig2>

      E.g.
      contig1 100 150 contig1 350 300
      contig1 4000 4050 contig2 110 60

      There is a script in the 'tools' directory of the package to convert SAM/BAM to a tab format.

      Regards,
      Boetsie
      Where should we report SSPACE bugs? The tools/sam_bam2Tab.pl script included in SSPACE basic v2.0 is extremely brittle - it cannot handle properly formatted paired end SAM/BAM (using the FLAG field), and instead expects fake single read data with read suffices (and only looks at the reverse complement bit of the FLAG). Thanks.

      Update: I've started hacking together a Python replacement https://github.com/peterjc/picobio/b..._sspace_tab.py
      Last edited by maubp; 02-20-2014, 02:54 AM. Reason: Adding link to alternative

      Comment


      • Hi Anna,

        How did you solve the problem with "Can't locate getopts.pl in @INC (@INC contains"
        I can't figure it out

        Thanks

        Comment


        • Hello fellows,

          I've been trying to assemble pre-assembled contigs into scaffolds using SSPACE. Everything goes ok up to a specific point where the software stops and gives me the following message:

          Bowtie-build error; -1 at /usr/local/src/sspace-basic/bin/PairingAndScaffolding.pl line 829.

          I tried running with -v option but no extra messages are output.
          Could please someone give me some help on how to debug that issue?

          I'm using a single fastq file of pre-assembled contigs and enabling contig extension ( -x 1 ) option, which should use a set of mate-pair reads I'm passing in on libraries.txt file.


          Thanks a lot.

          Cheers,
          George Condomitti.

          Comment


          • Hi Boetsie,

            Do you have a clue on how I could find out about that error?

            Thank you very much!

            Best rgrds,
            Conomitti.

            Comment


            • Sorry for the late reply. It seems that I do not get an e-mail notification anymore when a reply is given in this thread

              Please have a look at this post where colindaven suggests how to fix the problem;

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


              Boetsie

              Originally posted by condomitti View Post
              Hi Boetsie,

              Do you have a clue on how I could find out about that error?

              Thank you very much!

              Best rgrds,
              Conomitti.

              Comment


              • Originally posted by boetsie View Post
                Sorry for the late reply. It seems that I do not get an e-mail notification anymore when a reply is given in this thread

                Please have a look at this post where colindaven suggests how to fix the problem;

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


                Boetsie

                Thanks Boetsie! That seems to solve!!

                Cheers,
                Condomitti.

                Comment


                • I'd still like to use 454 PE data to scaffold my assembly. Is there some intermediate file or Bowtie like PE format that we can feed to SSPACE?

                  Comment


                  • Hi Boetsie

                    Hi hope you can help me. I have been trying to use sspace to scaffold a old assembly with a new paired end library, several times, using different parameters and keeps giving the following error:
                    "Thread 26 terminated abnormally: Can't open bwa output -- fatal"

                    Although in my library file I have bowtie as selected aligner

                    "lib1 bowtie ../corrected_trimmomatic_rubellus_1.fastq.gz ../corrected_trimmomatic_rubellus_2.fastq.gz 225 0.60 FR"

                    The last log message is that is doing the mapping stage

                    "2014: Mapping reads to contigs with Bowtie"

                    and that's it,

                    Can you please advise me where to look in order to try track the issue?

                    Thanking in advance

                    Luis

                    Comment


                    • Originally posted by luiscunhamx View Post
                      Hi Boetsie

                      Hi hope you can help me. I have been trying to use sspace to scaffold a old assembly with a new paired end library, several times, using different parameters and keeps giving the following error:
                      "Thread 26 terminated abnormally: Can't open bwa output -- fatal"

                      Although in my library file I have bowtie as selected aligner

                      "lib1 bowtie ../corrected_trimmomatic_rubellus_1.fastq.gz ../corrected_trimmomatic_rubellus_2.fastq.gz 225 0.60 FR"

                      The last log message is that is doing the mapping stage

                      "2014: Mapping reads to contigs with Bowtie"

                      and that's it,

                      Can you please advise me where to look in order to try track the issue?

                      Thanking in advance

                      Luis
                      I have same problem. Have you solved the problem, Luis?

                      Amery

                      Comment


                      • Hi Amery, sorry just saw your question now, the good news is yes I solve it but the bad news is that the solution was running it with only one processor. I would love to know if you could also solve it by other way, because it takes ages running it with only one thread.

                        Cheers

                        L

                        Comment


                        • hi boetsie,

                          i just stumbled across your program and would like to use it.
                          for the read libraries, how are they build? could you give an example of the steps to get there? i have read files in fastq format let's say R1.fastq and R2.fastq

                          thanks, julia

                          Comment


                          • Originally posted by luiscunhamx View Post
                            Hi Amery, sorry just saw your question now, the good news is yes I solve it but the bad news is that the solution was running it with only one processor. I would love to know if you could also solve it by other way, because it takes ages running it with only one thread.

                            Cheers

                            L
                            Just reporting the same problem with SSPACE_Standard_v3.0, but only when "-x 1". "-T 1" didn't fix it for me. Could the problem be in ExtendOrFormatContigs.pl if statements?

                            Code:
                            if($#filesbowtie >[B]=[/B] 0)
                            ..
                            if($#filesbwa >[B]=[/B] 0 || $#filesbwasw >[B]=[/B] 0){
                            Shouldn't these be ">=1" or ">0"? The program finishes just fine if I use bwa or bwasw instead of bowtie..
                            Last edited by rhinoceros; 05-29-2014, 09:32 PM.
                            savetherhino.org

                            Comment


                            • Is that possible to use already paired library for SSPACE

                              Hi all,

                              We normally use CLC bio for assembly, and CLC bio always pair the two Reads into one (i.e. R1.fastq + R2.fastq = pair.fastq).

                              However, in the SSPACE tutorial, the two Reads should be listed separately in the "library file" as: Lib1 bwa R1.fastq R2.fastq 4000 0.5 RF

                              Therefore, I would like to ask whether I could change the "library file" as:
                              Lib1 bwa pair.fastq 4000 0.5 RF

                              In addition, it seems that untrimmed reads could give better results than trimmed results. Do you trim your reads before running SSPACE?

                              Thank you for your helps in advance!

                              Cheers,

                              Kai

                              Comment


                              • waiting................

                                Comment

                                Latest Articles

                                Collapse

                                • 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
                                • seqadmin
                                  The Impact of AI in Genomic Medicine
                                  by seqadmin



                                  Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                                  02-26-2024, 02:07 PM

                                ad_right_rmr

                                Collapse

                                News

                                Collapse

                                Topics Statistics Last Post
                                Started by seqadmin, 03-14-2024, 06:13 AM
                                0 responses
                                32 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-08-2024, 08:03 AM
                                0 responses
                                72 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-07-2024, 08:13 AM
                                0 responses
                                80 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-06-2024, 09:51 AM
                                0 responses
                                68 views
                                0 likes
                                Last Post seqadmin  
                                Working...
                                X