Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Scaffolding suggestion? qqsmallfrog Bioinformatics 13 05-29-2013 08:50 PM
Need suggestion about RNA-seq shw Bioinformatics 0 12-12-2011 09:01 AM
suggestion and price Dan_cham SOLiD 11 03-04-2011 06:07 AM
New topic suggestion Susanne Site Feedback/Suggestions 2 10-16-2009 04:47 AM
some suggestion Chien-Yuan Chen Site Feedback/Suggestions 1 08-24-2009 03:30 PM

Thread Tools
Old 11-30-2010, 08:36 PM   #1
Junior Member
Location: China

Join Date: Oct 2010
Posts: 3
Default 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
ljhwahaha is offline   Reply With Quote
Old 11-28-2012, 10:00 AM   #2
Location: Uppsala

Join Date: Jan 2010
Posts: 25


I also want to do the same but could not find any software. Did you find any overlaping assembler?
~Adnan~ is offline   Reply With Quote
Old 11-29-2012, 04:39 AM   #3
Senior Member
Location: Berlin, DE

Join Date: May 2008
Posts: 628

Have a look at Celera Assembler,
sklages is offline   Reply With Quote
Old 03-19-2014, 04:02 AM   #4
Junior Member
Location: Italy

Join Date: Sep 2012
Posts: 4

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?
amaurizio is offline   Reply With Quote
Old 03-24-2014, 06:11 AM   #5
Senior Member
Location: Boston area

Join Date: Nov 2007
Posts: 747

# 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";
krobison is offline   Reply With Quote
Old 03-25-2014, 08:27 AM   #6
Junior Member
Location: Italy

Join Date: Sep 2012
Posts: 4

Thank you very much!
amaurizio is offline   Reply With Quote


Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:00 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO