http://pgi.re.kr/index.php?action=history&feed=atom&title=N50_and_N90 2026-04-10T22:44:48Z Revision history for this page on the wiki MediaWiki 1.31.3 http://pgi.re.kr/index.php?title=N50_and_N90&diff=3174&oldid=prev 2010-11-08T10:41:31Z

<p></p> <table class="diff diff-contentalign-left" data-mw="interface"> <col class="diff-marker" /> <col class="diff-content" /> <col class="diff-marker" /> <col class="diff-content" /> <tr class="diff-title" lang="en"> <td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td> <td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 10:41, 8 November 2010</td> </tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l6" >Line 6:</td> <td colspan="2" class="diff-lineno">Line 6:</td></tr> <tr><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td></tr> <tr><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;perl -e 'my ($len,$total)=(0,0);my @x;while(&amp;lt;&amp;gt;){if(/^[\&amp;gt;\@]/){if($len&amp;gt;0){$total+=$len;push@x,$len;};$len=0;}else{s/\s//g;$len+=length($_);}}if ($len&amp;gt;0){$total+=$len;push @x,$len;}@x=sort{$b&amp;lt;=&amp;gt;$a}@x; my ($count,$half)=(0,0);for (my $j=0;$j&amp;lt;@x;$j++){$count+=$x[$j];if(($count&amp;gt;=$total/2)&amp;amp;&amp;amp;($half==0)){print &amp;quot;N50: $x[$j]\n&amp;quot;;$half=$x[$j]}elsif($count&amp;gt;=$total*0.9){print &amp;quot;N90: $x[$j]\n&amp;quot;;exit;}}'&amp;nbsp; contigs.fahttp://seqanswers.com/forums/showthread.php?t=2766&lt;/p&gt;</div></td><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;perl -e 'my ($len,$total)=(0,0);my @x;while(&amp;lt;&amp;gt;){if(/^[\&amp;gt;\@]/){if($len&amp;gt;0){$total+=$len;push@x,$len;};$len=0;}else{s/\s//g;$len+=length($_);}}if ($len&amp;gt;0){$total+=$len;push @x,$len;}@x=sort{$b&amp;lt;=&amp;gt;$a}@x; my ($count,$half)=(0,0);for (my $j=0;$j&amp;lt;@x;$j++){$count+=$x[$j];if(($count&amp;gt;=$total/2)&amp;amp;&amp;amp;($half==0)){print &amp;quot;N50: $x[$j]\n&amp;quot;;$half=$x[$j]}elsif($count&amp;gt;=$total*0.9){print &amp;quot;N90: $x[$j]\n&amp;quot;;exit;}}'&amp;nbsp; contigs.fahttp://seqanswers.com/forums/showthread.php?t=2766&lt;/p&gt;</div></td></tr> <tr><td class='diff-marker'>−</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td><td class='diff-marker'>+</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;<ins class="diffchange diffchange-inline">[[N50 N90 calculation using Perl code]]</ins>&amp;nbsp;&lt;/p&gt;</div></td></tr> <tr><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td></tr> <tr><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&lt;a href=&quot;http://seqanswers.com/forums/showthread.php?t=2766&quot;&gt;http://seqanswers.com/forums/showthread.php?t=2766&lt;/a&gt;&lt;/p&gt;</div></td><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&lt;a href=&quot;http://seqanswers.com/forums/showthread.php?t=2766&quot;&gt;http://seqanswers.com/forums/showthread.php?t=2766&lt;/a&gt;&lt;/p&gt;</div></td></tr> <tr><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td><td class='diff-marker'>&#160;</td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>&lt;p&gt;&amp;nbsp;&lt;/p&gt;</div></td></tr> <!-- diff cache key pgi_re_kr:diff::1.12:old-3173:rev-3174 --> </table>

WikiSysop http://pgi.re.kr/index.php?title=N50_and_N90&diff=3173&oldid=prev 2010-11-08T10:40:56Z

<p>Created page with &quot;&lt;p&gt;The N50 contig size is a weighted median value and defined as &lt;br /&gt; the length of the smallest contig S in the sorted list of all &lt;br /&gt; contigs where the cumulative length f...&quot;</p> <p><b>New page</b></p><div>&lt;p&gt;The N50 contig size is a weighted median value and defined as &lt;br /&gt;<br /> the length of the smallest contig S in the sorted list of all &lt;br /&gt;<br /> contigs where the cumulative length from the largest contig to &lt;br /&gt;<br /> contig S is at least 50% of the total length.&lt;br /&gt;<br /> &amp;nbsp;&lt;/p&gt;<br /> &lt;p&gt;&amp;nbsp;&lt;/p&gt;<br /> &lt;p&gt;perl -e 'my ($len,$total)=(0,0);my @x;while(&amp;lt;&amp;gt;){if(/^[\&amp;gt;\@]/){if($len&amp;gt;0){$total+=$len;push@x,$len;};$len=0;}else{s/\s//g;$len+=length($_);}}if ($len&amp;gt;0){$total+=$len;push @x,$len;}@x=sort{$b&amp;lt;=&amp;gt;$a}@x; my ($count,$half)=(0,0);for (my $j=0;$j&amp;lt;@x;$j++){$count+=$x[$j];if(($count&amp;gt;=$total/2)&amp;amp;&amp;amp;($half==0)){print &amp;quot;N50: $x[$j]\n&amp;quot;;$half=$x[$j]}elsif($count&amp;gt;=$total*0.9){print &amp;quot;N90: $x[$j]\n&amp;quot;;exit;}}'&amp;nbsp; contigs.fahttp://seqanswers.com/forums/showthread.php?t=2766&lt;/p&gt;<br /> &lt;p&gt;&amp;nbsp;&lt;/p&gt;<br /> &lt;p&gt;&amp;nbsp;&lt;/p&gt;<br /> &lt;p&gt;&lt;a href=&quot;http://seqanswers.com/forums/showthread.php?t=2766&quot;&gt;http://seqanswers.com/forums/showthread.php?t=2766&lt;/a&gt;&lt;/p&gt;<br /> &lt;p&gt;&amp;nbsp;&lt;/p&gt;</div>

WikiSysop