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";}'
Posted by Xusheng at 2:30 PM

No comments:

Post a Comment

Newer Post Older Post Home
Subscribe to: Post Comments (Atom)

About Me

Xusheng
View my complete profile

Blog Archive

  • ►  2011 (8)
    • ►  September (1)
    • ►  August (1)
    • ►  July (4)
    • ►  May (1)
    • ►  February (1)
  • ▼  2010 (44)
    • ►  December (5)
    • ►  November (3)
    • ►  October (7)
    • ►  September (11)
    • ►  August (2)
    • ►  July (3)
    • ►  June (4)
    • ►  May (4)
    • ►  April (2)
    • ►  February (1)
    • ▼  January (2)
      • Find nnnn in genome
      • maq
  • ►  2009 (46)
    • ►  December (3)
    • ►  November (2)
    • ►  October (2)
    • ►  September (19)
    • ►  August (6)
    • ►  July (2)
    • ►  June (5)
    • ►  May (5)
    • ►  February (1)
    • ►  January (1)
  • ►  2008 (27)
    • ►  November (9)
    • ►  October (4)
    • ►  September (14)
Xusheng Wang. Simple theme. Powered by Blogger.