Perform the new stitching in the denovo pipeline

cfe-lab · Nov 15, 2023 · 50c2ccc · 50c2ccc
1 parent f39195a
commit 50c2ccc
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/micall/core/denovo.py b/micall/core/denovo.py
@@ -19,7 +19,7 @@
 from Bio.SeqRecord import SeqRecord
 
 from micall.core.project_config import ProjectConfig
-from micall.core.contig_stitcher import GenotypedContig
+from micall.core.contig_stitcher import GenotypedContig, stitch_consensus
 
 IVA = "iva"
 DEFAULT_DATABASE = os.path.join(os.path.dirname(__file__),
@@ -84,6 +84,7 @@ def write_contig_refs(contigs_fasta_path,
                          group_refs=group_refs)
 
     contigs = list(read_assembled_contigs(group_refs, genotypes, contigs_fasta_path))
+    contigs = list(stitch_consensus(contigs))
 
     for contig in contigs:
         writer.writerow(dict(ref=contig.ref_name,