Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Need some suggestion for overlap assembler

    Recently, i faced with a problem, I want to overlap the sequences in order to have a longer one, but i can't find out a assembler which could be suitable to my data.
    The length of my data is wide range, from 102bp to 87249bp, with N50 length of 738bp, so it count not be assembled using 454 data assembler, such as Newbler.
    What's more, the total bases number of my data is about 1.9G, so it seems to be impossible to use Sanger data assembler.
    What i want to do is to extend the sequences based on the overlap relationship between them, if someone have any suggestion, please inform me.
    thank you

  • #2
    hello,

    I also want to do the same but could not find any software. Did you find any overlaping assembler?
    ~Adnan~

    Comment


    • #3
      Have a look at Celera Assembler, http://sourceforge.net/apps/mediawik...itle=Main_Page

      Comment


      • #4
        Hello,
        I used Pear to merge my Illumina (MiSeq sequencer) pair-end reads.
        I'd like to know the overlap length in each of the merged reads in order to calculate min length, max, average, mode, median etc...
        Can you suggest my how to do?
        I tried with a perl script (length R1 + length R2 - length MergedR1R2) but I am not very good in programming....
        Can anybody help me and tell me how to do this?
        Thanks!!

        Comment


        • #5
          Code:
          #!/usr/bin/perl
          #!/usr/bin/perl
          # debugging left as exercise for student :-)
          use strict;
          
          foreach my $arg(@ARGV)
          {
             my %lengths=();
             # read sequence file for lengths
             my $rdr=new Bio::SeqIO(-file=>$arg,-format=>'fastq');
             my $sum=0; my $cnt=0;
             while (my $rec=$rdr->next_seq)
             {
                $lengths{$rec->length}++; $sum+=$rec->length; $cnt++;
             }
             my @sortedLengths=sort {$a<=>$b} keys %lengths;
             my $minLen=$sortedLengths[0];
             my $maxLen=$sortedLengths[$#sortedLengths];
             my $meanLen=$sum/$cnt;
             my ($modeLen)=sort {$lengths{$b}<=>$lengths{$a}} keys %lengths;
          
             my $medLen="error";
             # calculating median left as exercise
          
             print join("\t",$arg,$minLen,$maxLen,$meanLen,$modeLen,$medLen),"\n";
          }

          Comment


          • #6
            Thank you very much!

            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
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            51 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            68 views
            0 likes
            Last Post seqadmin  
            Working...
            X