Mercurial > repos > mkhan1980 > ctcf_analysis

comparison check2.pl @ 6:f1696b304b8d draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author mkhan1980
date Thu, 25 Apr 2013 11:22:23 -0400
parents 2cceb9398d33
children
comparison
equal deleted inserted replaced
5:1c09a4a56459 6:f1696b304b8d
104 $kount=0; 104 $kount=0;
105 $coun=0; 105 $coun=0;
106 106
107 # Define the pattern for CTCF. Because of pseudocount, a wildcard is allowed at 107 # Define the pattern for CTCF. Because of pseudocount, a wildcard is allowed at
108 #each position. 108 #each position.
109 my $pattern = "[ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC][ATGC]"; 109 my $pattern = "[ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN][ATGCN]";
110 110
111 # Compare the pattern with the 19 bp nucleotide segments. 111 # Compare the pattern with the 19 bp nucleotide segments.
112 112
113 while($var =~ m/$pattern/gi) 113 while($var =~ m/$pattern/gi)
114 { 114 {