Skip contigs not in range requested by the user

Schaudge · Dec 20, 2021 · fe6d364 · fe6d364
1 parent dd10d72
commit fe6d364
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/utgcns/utgcns.C b/src/utgcns/utgcns.C
@@ -160,6 +160,10 @@ processImportedTigs(cnsParameters  &params) {
 
   while (tig->importData(importFile, reads, importedLayouts, importedReads) == true) {
 
+    if (tig->tigID() < params.tigBgn || tig->tigID() > params.tigEnd) {
+       fprintf(stdout, "Skipping tig %d\n", tig->tigID());
+       continue;
+    }
     //  Log that we're processing.
 
     if (tig->numberOfChildren() > 1) {