src/hg/encode/encodeValidate/doEncodeValidate.pl 1.163
1.163 2009/03/06 20:47:31 mikep
diff header coming from GIS
Index: src/hg/encode/encodeValidate/doEncodeValidate.pl
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/encode/encodeValidate/doEncodeValidate.pl,v
retrieving revision 1.162
retrieving revision 1.163
diff -b -B -U 4 -r1.162 -r1.163
--- src/hg/encode/encodeValidate/doEncodeValidate.pl 6 Mar 2009 20:36:18 -0000 1.162
+++ src/hg/encode/encodeValidate/doEncodeValidate.pl 6 Mar 2009 20:47:31 -0000 1.163
@@ -811,15 +811,20 @@
# >461_19_90_F3
# T203033330010111011221200302001
# >461_19_209_F3
# T022213002230311203200200322000
+
+ # Files from GIS have this header:
+ # >920_22_656_F3,1.-152654094.1.35.35.0###,19.43558664.1.35.35.0###
+ # T01301010111200210102321210100112312
+
my ($path, $file, $type) = @_;
doTime("beginning validateCsfasta") if $opt_timing;
my $fh = openUtil($path, $file);
my $line = 0;
my $state = 'header';
my $seqName;
- my $states = {header => {REGEX => "^>\\d+_\\d+_\\d+_\.\\d+", NEXT => 'seq'},
+ my $states = {header => {REGEX => "^>\\d+_\\d+_\\d+_\.\\d+.*", NEXT => 'seq'},
seq => {REGEX => "^T\\d+", NEXT => 'header'},
};
while(<$fh>) {
chomp;