F-UT-URE
Thursday, January 14, 2010
Find nnnn in genome
perl -e'open(FILE,"chr1.fa");my $content=""; while(
){chomp $_; $content.=$_;} while($content=~/(NNNNNNN.*?N[A|C|G|T|a|c|g|t])/g){print length($1),"\t",pos($content),"\n";}'
No comments:
Post a Comment
Newer Post
Older Post
Home
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment