src/hg/tcga/pBamBam/asBamBam.c 1.5

1.5 2010/03/07 03:28:00 jsanborn
updated
Index: src/hg/tcga/pBamBam/asBamBam.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/tcga/pBamBam/asBamBam.c,v
retrieving revision 1.4
retrieving revision 1.5
diff -b -B -U 4 -r1.4 -r1.5
--- src/hg/tcga/pBamBam/asBamBam.c	7 Mar 2010 03:26:33 -0000	1.4
+++ src/hg/tcga/pBamBam/asBamBam.c	7 Mar 2010 03:28:00 -0000	1.5
@@ -26,15 +26,15 @@
   "asBamBam - Allele-Specific BamBam.\n"
   "usage:\n"
   "   asBamBam [options] <left.bam> <right.bam>\n"
   "options:\n"
-  "   -position=chrX:1-100 - Position to do analysis. [None]\n" 
+  "   -position=chr1:1-100 - Position to do analysis. [NULL=process all reads]\n" 
   "   -minSupport=INT   - Minimum Support in Germline. [4]\n"
   "   -minQ=INT         - Minimum acceptable base quality. [10]\n"
   "   -minMapQ=INT      - Minimum acceptable mapping quality. [20]\n"
   "   -avgMapQ=INT      - Minimum acceptable average mapping quality. [30]\n"
   "   -minChiSq=INT     - Minimum acceptable chi square of nucleotide freqs. [250]\n"
-  "   -maxISize=INT     - Maximum Insert Size. [10000]\n"
+  "   -maxISize=INT     - Maximum Insert Size. [2000]\n"
   "   -readsPerBin=INT  - Reads per CNV bin [2500]\n"
   "\n"
   );
 }
@@ -58,9 +58,9 @@
 int minSupport  = 4;
 int minQ        = 10;
 int minMapQ     = 20;
 int avgMapQ     = 30;
-int maxISize    = 10000;
+int maxISize    = 2000;
 int minChiSq    = 250;
 int readsPerBin = 2500;
 
 typedef struct {
@@ -1442,9 +1442,9 @@
     minMapQ = optionInt("minMapQ", 20);
 if (optionExists("avgMapQ"))
     avgMapQ = optionInt("avgMapQ", 30);
 if (optionExists("maxISize"))
-    maxISize = optionInt("maxISize", 10000);
+    maxISize = optionInt("maxISize", 2000);
 if (optionExists("minChiSq"))
     minChiSq = optionInt("minChiSq", 250);
 if (optionExists("readsPerBin"))
     readsPerBin = optionInt("readsPerBin", 2500);