R脚本:N50.N90.sort.R #!/usr/bin/env Rscriptsetwd("Bin_all/Bin_summary/")data=read.table("N50.N90.txt",header=F)BinID=vector()N50=vector()N90=vector()index=seq(from=1,to=length(data[,1]),by=3)a=1for(iinindex){BinID[a]=as.character(data[i,"V1"])N50[a]=as.character(da...
cat n50_n90_oneline.pl #!/usr/bin/perl -w my ($len,$total)=(0,0);my @x;while(<>){if(/^[\>\@]/){if($len>0){$total+=$len;push@x,$len;};$len=0;}else{s/\s//g;$len+=length($_);}}if ($len>0){$total+=$len;push @x,$len;}@x=sort{$b<=>$a}@x; my (...