00001 #include "log.h"
00002 #include "format.h"
00003 #include "gfr.h"
00004
00005
00006
00007 int main (int argc, char *argv[])
00008 {
00009 GfrEntry *currGE;
00010 int count;
00011 int countRemoved;
00012
00013 count = 0;
00014 countRemoved = 0;
00015 gfr_init ("-");
00016 puts (gfr_writeHeader ());
00017 while (currGE = gfr_nextEntry ()) {
00018 if (strEqual( currGE->chromosomeTranscript1, "chrM") ||
00019 strEqual( currGE->chromosomeTranscript2, "chrM") ) {
00020 countRemoved++;
00021 continue;
00022 } else {
00023 puts (gfr_writeGfrEntry (currGE));
00024 count++;
00025 }
00026 }
00027 gfr_deInit ();
00028 warn ("%s_numRemoved: %d",argv[0],countRemoved);
00029 warn ("%s_numGfrEntries: %d",argv[0],count);
00030 return 0;
00031 }
00032
