Merge pull request #264 from ncsa/fix/generate-reads-contig-shorter-than-read-len

Fix infinite loop when contig is shorter than read_len

Author: joshfactorial

neat/read_simulator/utils/generate_reads.py

@@ -183,6 +183,13 @@ def generate_reads(
     # _LOG.info(f'Sampling reads for thread {thread_index}...')
     start_time = time.time()
 
+    if len(reference) < options.read_len:
+        _LOG.warning(
+            f"Contig '{contig_name}' (length {len(reference)}) is shorter than read_len "
+            f"({options.read_len}). Skipping contig."
+        )
+        return []
+
     # _LOG.debug("Covering dataset.")
     t = time.time()
     reads = cover_dataset(