GPatch 0.3.6__tar.gz → 0.4.0__tar.gz

{GPatch-0.3.6/src/GPatch.egg-info → gpatch-0.4.0}/PKG-INFO

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
-Name: GPatch
-Version: 0.3.6
-Summary: Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.
+Name: gpatch
+Version: 0.4.0
+Summary: Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence. Download the github repository for helper scripts to automate GPatch workflows, identify and correct misjoins in the contig assembly, produce dot-plots of patched pseudoassemblies to a reference assembly, and generate chrom.sizes and liftover chains for patched pseudoassemblies.
 Home-page: https://github.com/adadiehl/GPatch
 Author: Adam Diehl
 Author-email: adadiehl@umich.edu
@@ -40,11 +40,10 @@ Description: # GPatch
         
         ## Usage
         ```
-        usage: GPatch.py [-h] -q SAM/BAM -r FASTA [-x STR] [-b FILENAME] [-m N]
-                         [-w PATH] [-d]
+        GPatch.py [-h] -q SAM/BAM -r FASTA [-x STR] [-b FILENAME] [-m N] [-w PATH] [-d] [-s] [-l N] [-e] [-k]
         ```
         
-        Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. Reference chromosomes with no mapped contigs are printed to output unchanged.
+        Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. By default, reference chromosomes with no mapped contigs are printed to output unchanged. Use the --drop_missing option to disable this behavior. By default, patches are applied to the 5' and 3' telomere ends of pseudochromsomes if the first and last contig alignments do not extend to the start/end of the reference chromsome. In some cases, this may cause spurious duplications. Use the --no_extend option if this is a concern. Note that GPatch is designed to be run on single-haplotype or unphased genome assemblies. For phased assemblies, each haplotype should be separated into its own input FASTA file prior to alignment. GPatch can then be run separately on the BAM files for each haplotype to obtain phased pseudoassemblies, otherwise results will be unpredictable and likely incorrect.
         
         #### Required Arguments
         | Argument | Description |
@@ -61,7 +60,10 @@ Description: # GPatch
         | __-m N, --min_qual_score N__ | Minimum mapping quality score to retain an alignment. Default=30 |
         | __-w PATH, --whitelist PATH__ | Path to BED file containing whitelist regions: i.e., the inverse of blacklist regions. Supplying this will have the effect of excluding alignments that fall entirely within blacklist regions. Default=None |
         | __-d, --drop_missing__ | Omit unpatched reference chromosome records from the output if no contigs map to them. Default: Unpatched chromosomes are printed to output unchanged. |
-        
+        | __-s, --scaffold_only__ | Pad gaps between placed contigs with strings of N characters instead of patching with sequence from the reference assembly. Effectively turns GPatch into a reference-guided scaffolding tool. Note that patches.bed will still be generated to document (inverse) mapped contig boundaries in reference frame. |
+        | __-l N, --gap_length N__ | Length of "N" gaps separating placed gontigs when using --scaffold_only. Has no effect when in default patching mode. Default=Estimate gap length from alignment. |
+        | __-e, --no_extend__ | Do not patch telomere ends of pseudochromosomes with reference sequence upstream of the first mapped contig and downstream of the last mapped contig. Default is to include 5' and 3' patches to extend telomeres to the ends implied by the alignment. |
+        | __-k, --keep_nested__ | Do not drop contigs with mapped positions nested entirely inside other mapped contigs. Instead, these will be bookended after the contig in which they are nested. Default is to drop contigs with mapped positions nested entirely within other mapped contigs. This option should be used with caution as these mappings cannot be placed unambigiously relative to other mapped contigs, thus including them is likely to lead to unpredictable and possibly incorrect results. Do not use this unless you are sure you know what you are doing! |
         
         ## Output
         
@@ -72,11 +74,17 @@ Description: # GPatch
         | __contigs.bed__ | Location of contigs in the coordinate frame of the patched genome. |
         | __patches.bed__ | Location of patches in the coordinate frame of the reference genome. |
         
+        ## Helper Scripts
+        
+        The scripts directory contains helper scripts for running GPatch and working with its output, including a shell script to automate a two-stage patching process, including initial alignment and patching steps, misjoin breakpoint prediction and contig-breaking, and subsequent realignment and patching of the split contigs, with dot-plots against the reference assembly created after both patching stages. In addition, scripts are provided to generate a chrom.sizes file for a patched pseudoassembly and a set of liftover chains that can be used to translate features mapped to the unpatched contigs to the patched pseudoassembly. Please see the README.md within the scripts for detailed usage information.
         
         ## Citing GPatch
         Please use the following citation if you use this software in your work:
         
-        CITATION_HERE
+        Fast and Accurate Draft Genome Patching with GPatch
+        Adam Diehl, Alan Boyle
+        bioRxiv 2025.05.22.655567; doi: https://doi.org/10.1101/2025.05.22.655567
+        
 Keywords: genomics,genome assembly
 Platform: UNKNOWN
 Classifier: Development Status :: 4 - Beta

gpatch-0.4.0/README.md

@@ -0,0 +1,78 @@
+# GPatch
+## Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.
+
+Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference.
+
+## Dependencies
+* Python >= v3.7
+* samtools (https://github.com/samtools/samtools)
+* biopython (https://biopython.org/)
+* pysam (https://github.com/pysam-developers/pysam)
+* minimap2 (https://github.com/lh3/minimap2)
+
+We recommend using minimap2 for alignment, using the -a option to generate SAM output.
+
+## Installation
+
+We recommend installing with conda, into a new environment:
+```
+conda create -n GPatch -c conda-forge -c bioconda Bio pysam minimap2 samtools GPatch
+```
+
+Install with pip:
+```
+pip install GPatch
+```
+
+Installation from the github repository is not recommended. However, if you must, follow the steps below:
+1) git clone https://github.com/adadiehl/GPatch
+2) cd GPatch/
+3) python3 -m pip install -e .
+
+
+## Usage
+```
+GPatch.py [-h] -q SAM/BAM -r FASTA [-x STR] [-b FILENAME] [-m N] [-w PATH] [-d] [-s] [-l N] [-e] [-k]
+```
+
+Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. By default, reference chromosomes with no mapped contigs are printed to output unchanged. Use the --drop_missing option to disable this behavior. By default, patches are applied to the 5' and 3' telomere ends of pseudochromsomes if the first and last contig alignments do not extend to the start/end of the reference chromsome. In some cases, this may cause spurious duplications. Use the --no_extend option if this is a concern. Note that GPatch is designed to be run on single-haplotype or unphased genome assemblies. For phased assemblies, each haplotype should be separated into its own input FASTA file prior to alignment. GPatch can then be run separately on the BAM files for each haplotype to obtain phased pseudoassemblies, otherwise results will be unpredictable and likely incorrect.
+
+#### Required Arguments
+| Argument | Description |
+|---|---|
+| __-q SAM/BAM, --query_bam SAM/BAM__ | Path to SAM/BAM file containing non-overlapping contig mappings to the reference genome. |
+| __-r FASTA, --reference_fasta FASTA__ | Path to reference genome fasta. |
+
+#### Optional Arguments:
+| Argument | Description |
+|---|---|
+| __-h, --help__ | Show this help message and exit. |
+| __-x STR, --prefix STR__ | Prefix to add to output file names. Default=None |
+| __-b FILENAME, --store_final_bam FILENAME__ | Store the final set of primary contig alignments to the given file name. Default: Do not store the final BAM. |
+| __-m N, --min_qual_score N__ | Minimum mapping quality score to retain an alignment. Default=30 |
+| __-w PATH, --whitelist PATH__ | Path to BED file containing whitelist regions: i.e., the inverse of blacklist regions. Supplying this will have the effect of excluding alignments that fall entirely within blacklist regions. Default=None |
+| __-d, --drop_missing__ | Omit unpatched reference chromosome records from the output if no contigs map to them. Default: Unpatched chromosomes are printed to output unchanged. |
+| __-s, --scaffold_only__ | Pad gaps between placed contigs with strings of N characters instead of patching with sequence from the reference assembly. Effectively turns GPatch into a reference-guided scaffolding tool. Note that patches.bed will still be generated to document (inverse) mapped contig boundaries in reference frame. |
+| __-l N, --gap_length N__ | Length of "N" gaps separating placed gontigs when using --scaffold_only. Has no effect when in default patching mode. Default=Estimate gap length from alignment. |
+| __-e, --no_extend__ | Do not patch telomere ends of pseudochromosomes with reference sequence upstream of the first mapped contig and downstream of the last mapped contig. Default is to include 5' and 3' patches to extend telomeres to the ends implied by the alignment. |
+| __-k, --keep_nested__ | Do not drop contigs with mapped positions nested entirely inside other mapped contigs. Instead, these will be bookended after the contig in which they are nested. Default is to drop contigs with mapped positions nested entirely within other mapped contigs. This option should be used with caution as these mappings cannot be placed unambigiously relative to other mapped contigs, thus including them is likely to lead to unpredictable and possibly incorrect results. Do not use this unless you are sure you know what you are doing! |
+
+## Output
+
+GPatch produces three output files:
+| File | Description |
+|---|---|
+| __patched.fasta__ | The final patched genome. |
+| __contigs.bed__ | Location of contigs in the coordinate frame of the patched genome. |
+| __patches.bed__ | Location of patches in the coordinate frame of the reference genome. |
+
+## Helper Scripts
+
+The scripts directory contains helper scripts for running GPatch and working with its output, including a shell script to automate a two-stage patching process, including initial alignment and patching steps, misjoin breakpoint prediction and contig-breaking, and subsequent realignment and patching of the split contigs, with dot-plots against the reference assembly created after both patching stages. In addition, scripts are provided to generate a chrom.sizes file for a patched pseudoassembly and a set of liftover chains that can be used to translate features mapped to the unpatched contigs to the patched pseudoassembly. Please see the README.md within the scripts for detailed usage information.
+
+## Citing GPatch
+Please use the following citation if you use this software in your work:
+
+Fast and Accurate Draft Genome Patching with GPatch
+Adam Diehl, Alan Boyle
+bioRxiv 2025.05.22.655567; doi: https://doi.org/10.1101/2025.05.22.655567

{GPatch-0.3.6 → gpatch-0.4.0}/setup.py

@@ -28,11 +28,11 @@ def readme():
 def main():
 
     metadata = dict(
-        name="GPatch",
-        version="0.3.6",
+        name="gpatch",
+        version="0.4.0",
         author="Adam Diehl",
         author_email="adadiehl@umich.edu",
-        description="Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.",
+        description="Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence. Download the github repository for helper scripts to automate GPatch workflows, identify and correct misjoins in the contig assembly, produce dot-plots of patched pseudoassemblies to a reference assembly, and generate chrom.sizes and liftover chains for patched pseudoassemblies.",
         long_description=readme(),
         long_description_content_type="text/markdown",
         url="https://github.com/adadiehl/GPatch",

{GPatch-0.3.6 → gpatch-0.4.0}/src/GPatch/GPatch.py

@@ -1,13 +1,13 @@
 #!/usr/bin/env python3
 
 import sys, os, re
-from argparse import ArgumentParser
+from argparse import ArgumentParser, SUPPRESS
 from Bio import SeqIO
 from Bio.Seq import Seq
 import pysam
 import random
 
-__version__ = "0.3.6"
+__version__ = "0.4.0"
 
 """
 Given a name-sorted BAM file, cluster mapped reads by name and
@@ -21,16 +21,25 @@ def create_primary_alignments_list(query_bam, min_qual_score):
     Extract all primary alignments from the BAM file and return
     as a list.
     """
+    #i = 0
+    #j = 0
+    #k = 0
+    #l = 0
     ret = []
     for aln in query_bam:
+        #i += 1
         if aln.is_unmapped:
             # Skip unmapped contigs
+            #j += 1
             continue
         if aln.mapping_quality < min_qual_score:
             # Skip low-quality mappings
+            #k += 1
             continue
         if not (aln.is_supplementary or aln.is_secondary):
+            #l += 1
             ret.append(aln)
+    #sys.stderr.write("%d\t%d\t%d\t%d\t%d\n" % (i, j, k, l, len(ret)))
     return ret
 
 
@@ -139,6 +148,7 @@ def cluster_contig_positions(sorted_contig_intervals, drop_nested=False):
                         cluster.append(contig)
                     else:
                         #sys.stderr.write("\tnested\n")
+                        #sys.stderr.write("\t%s\n\t%s\n" % (cluster[-1], contig))
                         pass
                 else:
                     cluster.append(contig)
@@ -151,7 +161,7 @@ def cluster_contig_positions(sorted_contig_intervals, drop_nested=False):
 
 
 def main():
-    parser = ArgumentParser(description="Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference.")
+    parser = ArgumentParser(description="Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. By default, reference chromosomes with no mapped contigs are printed to output unchanged. Use the --drop_missing option to disable this behavior. By default, patches are applied to the 5' and 3' telomere ends of pseudochromsomes if the first and last contig alignments do not extend to the start/end of the reference chromsome. In some cases, this may cause spurious duplications. Use the --no_extend option if this is a concern. Note that GPatch is designed to be run on single-haplotype or unphased genome assemblies. For phased assemblies, each haplotype should be separated into its own input FASTA file prior to alignment. GPatch can then be run separately on the BAM files for each haplotype to obtain phased pseudoassemblies, otherwise results will be unpredictable and likely incorrect.")
     parser.add_argument('-q', '--query_bam', metavar='SAM/BAM', type=str,
                         required=True, help='Path to SAM/BAM file containing non-overlapping contig mappings to the reference genome.')
     parser.add_argument('-r', '--reference_fasta', metavar='FASTA', type=str,
@@ -171,6 +181,25 @@ def main():
     parser.add_argument('-d', '--drop_missing',
                         required=False, default=False, action="store_true",
                         help='Omit unpatched reference chromosome records from the output if no contigs map to them. Default: Unpatched chromosomes are printed to output unchanged.')
+    # Note the --no_trim option is deprecated in versions 0.4.0. Trimming
+    # is disabled as it had universally negative effects in our testing
+    # and should not be used.
+    parser.add_argument('-t', '--no_trim',
+                        required=False, default=True, action="store_true",
+                        #help='Do not trim the 5-prime end of contigs whose mappings overlap the previously-placed contig. Default: Overlapping contig sequence will be trimmed at the previous 3-prime contig breakpoint.')
+                        help=SUPPRESS)
+    parser.add_argument('-s', '--scaffold_only',
+                        required=False, default=False, action="store_true",
+                        help='Pad gaps between placed contigs with strings of N characters instead of patching with sequence from the reference assembly. Effectively turns GPatch into a reference-guided scaffolding tool. Note that patches.bed will still be generated to document (inverse) mapped contig boundaries in reference frame.')
+    parser.add_argument('-l', '--gap_length', metavar='N', type=int,
+                        required=False, default=-1,
+                        help='Length of "N" gaps separating placed gontigs when using --scaffold_only. Has no effect when in default patching mode. Default=Estimate gap length from alignment.')
+    parser.add_argument('-e', '--no_extend',
+                        required=False, default=False, action="store_true",
+                        help='Do not patch telomere ends of pseudochromosomes with reference sequence upstream of the first mapped contig and downstream of the last mapped contig. Default is to include 5\' and 3\' patches to extend telomeres to the ends implied by the alignment.')
+    parser.add_argument('-k', '--keep_nested',
+                        required=False, default=False, action="store_true",
+                        help='Do not drop contigs with mapped positions nested entirely inside other mapped contigs. Instead, these will be bookended after the contig in which they are nested. Default is to drop contigs with mapped positions nested entirely within other mapped contigs. This option should be used with caution as these mappings cannot be placed unambigiously relative to other mapped contigs, thus including them is likely to lead to unpredictable and possibly incorrect results. Do not use this unless you are sure you know what you are doing!')
     
     args = parser.parse_args()
 
@@ -201,19 +230,24 @@ def main():
         os.rename(filtered_bam + '.bai', sorted_primary_alignments_bam + ".bai")
 
     # Check for nested/overlapping alignment spaces among the remaining contigs.
-    sys.stderr.write("Checking for nested or overlapping contigs...\n")
+    if not args.keep_nested:
+        sys.stderr.write("Checking for nested contigs...\n")
     sorted_primary_alignments = pysam.AlignmentFile(sorted_primary_alignments_bam, "rb")
     final_alignments = []
     for ref_seq in SeqIO.parse(args.reference_fasta, "fasta"):
         # Select all contigs mapped to this sequence.
         contigs = list(sorted_primary_alignments.fetch(ref_seq.id))
+        #sys.stderr.write("%s\t%s\n" % (ref_seq.id, contigs))
 
         # Convert the list of contigs into a list of tuples: (start, end, name),
         # representing the inferred contig breakpoints, sorted by start, then end position.
         sorted_contig_intervals = contigs_to_sorted_interval_list(contigs, contig_breakpoints)
         
         # Cluster the linear intervals on overlap, leaving out nested contigs
-        clusters = cluster_contig_positions(sorted_contig_intervals, drop_nested=True)
+        drop_nested = True
+        if args.keep_nested:
+            drop_nested = False
+        clusters = cluster_contig_positions(sorted_contig_intervals, drop_nested=drop_nested)
 
         # Next thing to do is to push all retained contigs across clusters into
         # a final list of useful contigs. Since pysam.fetch does not allow
@@ -229,9 +263,13 @@ def main():
                 final_alignments.append(contig)
 
     sorted_primary_alignments.close()
-                    
+
+    # Clean up the initial sorted BAM and index files.
+    os.remove(sorted_primary_alignments_bam)
+    os.remove(sorted_primary_alignments_bam + '.bai')
+    
     # Write the final set of alignments to a sorted BAM file
-    sorted_primary_alignments_bam = sorted_bam_from_aln_list(final_alignments, query_bam.header, sorted_bam_name=args.store_final_bam)                
+    sorted_primary_alignments_bam = sorted_bam_from_aln_list(final_alignments, query_bam.header, sorted_bam_name=args.store_final_bam)
     pysam.index(sorted_primary_alignments_bam)
         
     # Set up output streams for fasta, patches.bed, and contigs.bed
@@ -281,15 +319,39 @@ def main():
             if rstart > pos:
                 # No overlap. Append the patch, in all lower-case
                 # for easy identification.
-                patched_seq = patched_seq + ref_seq.seq[pos:rstart].lower()
-                # Write patch coordinates in reference frame to patches_bed
-                patches_bed.write("%s\t%d\t%d\n" % (ref_seq.id, pos, rstart))
+                patch = ""
+                if args.scaffold_only:
+                    # Use a string of N characters the length of the patch instead
+                    # of an actual patch.
+                    if args.gap_length >= 0:
+                        patch = "n" * args.gap_length
+                    else:
+                        patch = "n" * (rstart-pos)
+                else:
+                    patch = ref_seq.seq[pos:rstart].lower()
+                if pos == 0 and args.no_extend:
+                    pass
+                else:
+                    patched_seq = patched_seq + patch
+                    # Write patch coordinates in reference frame to patches_bed
+                    patches_bed.write("%s\t%d\t%d\n" % (ref_seq.id, pos, rstart))
             else:
                 # Handle overlapping contig ends by trimming the 5' end of
                 # this contig sequence by the length of the overlap. Note this
                 # evaluates to zero if contigs are bookended (i.e., pos and
                 # rstart are equal.)
-                qstart = pos - rstart
+                #
+                # NOTE THAT THIS BEHAVIOR IS DEPRECATED AND THE OPTION TO TRIM
+                # IS DISABLED IN VERSIONS 0.4.0 AND GREATER.
+                if not args.no_trim:
+                    qstart = pos - rstart
+                    if qstart > len(Seq(contig.query_sequence)):
+                        # This should not happen, as it is indicative of a nested
+                        # mapping, which should be filtered out prior to sequence-
+                        # building. However, I have seen instances of nested mappings
+                        # that somehow slip through the cracks. This is a hack to fix
+                        # those instances until I can figure out why it happens.
+                        continue
                 
             # Append the contig sequence to the patched sequence string
             contig_start = len(patched_seq)
@@ -316,7 +378,7 @@ def main():
                 
         # Once the above loop finishes, we need to add the terminal segment
         # from the reference genome.
-        if pos < len(ref_seq.seq):
+        if pos < len(ref_seq.seq) and not args.no_extend:
             #sys.stderr.write("%s\n" % (ref_seq.seq[pos:len(ref_seq.seq)]))
             patched_seq = patched_seq + ref_seq.seq[pos:len(ref_seq.seq)]
             patches_bed.write("%s\t%d\t%d\n" % (ref_seq.id, pos, len(ref_seq.seq)))

{GPatch-0.3.6 → gpatch-0.4.0/src/gpatch.egg-info}/PKG-INFO

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
-Name: GPatch
-Version: 0.3.6
-Summary: Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.
+Name: gpatch
+Version: 0.4.0
+Summary: Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence. Download the github repository for helper scripts to automate GPatch workflows, identify and correct misjoins in the contig assembly, produce dot-plots of patched pseudoassemblies to a reference assembly, and generate chrom.sizes and liftover chains for patched pseudoassemblies.
 Home-page: https://github.com/adadiehl/GPatch
 Author: Adam Diehl
 Author-email: adadiehl@umich.edu
@@ -40,11 +40,10 @@ Description: # GPatch
         
         ## Usage
         ```
-        usage: GPatch.py [-h] -q SAM/BAM -r FASTA [-x STR] [-b FILENAME] [-m N]
-                         [-w PATH] [-d]
+        GPatch.py [-h] -q SAM/BAM -r FASTA [-x STR] [-b FILENAME] [-m N] [-w PATH] [-d] [-s] [-l N] [-e] [-k]
         ```
         
-        Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. Reference chromosomes with no mapped contigs are printed to output unchanged.
+        Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. By default, reference chromosomes with no mapped contigs are printed to output unchanged. Use the --drop_missing option to disable this behavior. By default, patches are applied to the 5' and 3' telomere ends of pseudochromsomes if the first and last contig alignments do not extend to the start/end of the reference chromsome. In some cases, this may cause spurious duplications. Use the --no_extend option if this is a concern. Note that GPatch is designed to be run on single-haplotype or unphased genome assemblies. For phased assemblies, each haplotype should be separated into its own input FASTA file prior to alignment. GPatch can then be run separately on the BAM files for each haplotype to obtain phased pseudoassemblies, otherwise results will be unpredictable and likely incorrect.
         
         #### Required Arguments
         | Argument | Description |
@@ -61,7 +60,10 @@ Description: # GPatch
         | __-m N, --min_qual_score N__ | Minimum mapping quality score to retain an alignment. Default=30 |
         | __-w PATH, --whitelist PATH__ | Path to BED file containing whitelist regions: i.e., the inverse of blacklist regions. Supplying this will have the effect of excluding alignments that fall entirely within blacklist regions. Default=None |
         | __-d, --drop_missing__ | Omit unpatched reference chromosome records from the output if no contigs map to them. Default: Unpatched chromosomes are printed to output unchanged. |
-        
+        | __-s, --scaffold_only__ | Pad gaps between placed contigs with strings of N characters instead of patching with sequence from the reference assembly. Effectively turns GPatch into a reference-guided scaffolding tool. Note that patches.bed will still be generated to document (inverse) mapped contig boundaries in reference frame. |
+        | __-l N, --gap_length N__ | Length of "N" gaps separating placed gontigs when using --scaffold_only. Has no effect when in default patching mode. Default=Estimate gap length from alignment. |
+        | __-e, --no_extend__ | Do not patch telomere ends of pseudochromosomes with reference sequence upstream of the first mapped contig and downstream of the last mapped contig. Default is to include 5' and 3' patches to extend telomeres to the ends implied by the alignment. |
+        | __-k, --keep_nested__ | Do not drop contigs with mapped positions nested entirely inside other mapped contigs. Instead, these will be bookended after the contig in which they are nested. Default is to drop contigs with mapped positions nested entirely within other mapped contigs. This option should be used with caution as these mappings cannot be placed unambigiously relative to other mapped contigs, thus including them is likely to lead to unpredictable and possibly incorrect results. Do not use this unless you are sure you know what you are doing! |
         
         ## Output
         
@@ -72,11 +74,17 @@ Description: # GPatch
         | __contigs.bed__ | Location of contigs in the coordinate frame of the patched genome. |
         | __patches.bed__ | Location of patches in the coordinate frame of the reference genome. |
         
+        ## Helper Scripts
+        
+        The scripts directory contains helper scripts for running GPatch and working with its output, including a shell script to automate a two-stage patching process, including initial alignment and patching steps, misjoin breakpoint prediction and contig-breaking, and subsequent realignment and patching of the split contigs, with dot-plots against the reference assembly created after both patching stages. In addition, scripts are provided to generate a chrom.sizes file for a patched pseudoassembly and a set of liftover chains that can be used to translate features mapped to the unpatched contigs to the patched pseudoassembly. Please see the README.md within the scripts for detailed usage information.
         
         ## Citing GPatch
         Please use the following citation if you use this software in your work:
         
-        CITATION_HERE
+        Fast and Accurate Draft Genome Patching with GPatch
+        Adam Diehl, Alan Boyle
+        bioRxiv 2025.05.22.655567; doi: https://doi.org/10.1101/2025.05.22.655567
+        
 Keywords: genomics,genome assembly
 Platform: UNKNOWN
 Classifier: Development Status :: 4 - Beta

gpatch-0.4.0/src/gpatch.egg-info/SOURCES.txt

@@ -0,0 +1,17 @@
+MANIFEST.in
+README.md
+pyproject.toml
+setup.py
+./CODE_OF_CONDUCT.md
+./LICENSE
+./requirements.txt
+src/GPatch/GPatch.py
+src/GPatch/__init__.py
+src/GPatch/__main__.py
+src/gpatch.egg-info/PKG-INFO
+src/gpatch.egg-info/SOURCES.txt
+src/gpatch.egg-info/dependency_links.txt
+src/gpatch.egg-info/entry_points.txt
+src/gpatch.egg-info/not-zip-safe
+src/gpatch.egg-info/requires.txt
+src/gpatch.egg-info/top_level.txt

GPatch-0.3.6/README.md

@@ -1,71 +0,0 @@
-# GPatch
-## Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.
-
-Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference.
-
-## Dependencies
-* Python >= v3.7
-* samtools (https://github.com/samtools/samtools)
-* biopython (https://biopython.org/)
-* pysam (https://github.com/pysam-developers/pysam)
-* minimap2 (https://github.com/lh3/minimap2)
-
-We recommend using minimap2 for alignment, using the -a option to generate SAM output.
-
-## Installation
-
-We recommend installing with conda, into a new environment:
-```
-conda create -n GPatch -c conda-forge -c bioconda Bio pysam minimap2 samtools GPatch
-```
-
-Install with pip:
-```
-pip install GPatch
-```
-
-Installation from the github repository is not recommended. However, if you must, follow the steps below:
-1) git clone https://github.com/adadiehl/GPatch
-2) cd GPatch/
-3) python3 -m pip install -e .
-
-
-## Usage
-```
-usage: GPatch.py [-h] -q SAM/BAM -r FASTA [-x STR] [-b FILENAME] [-m N]
-                 [-w PATH] [-d]
-```
-
-Starting with alignments of contigs to a reference genome, produce a chromosome-scale pseudoassembly by patching gaps between mapped contigs with sequences from the reference. Reference chromosomes with no mapped contigs are printed to output unchanged.
-
-#### Required Arguments
-| Argument | Description |
-|---|---|
-| __-q SAM/BAM, --query_bam SAM/BAM__ | Path to SAM/BAM file containing non-overlapping contig mappings to the reference genome. |
-| __-r FASTA, --reference_fasta FASTA__ | Path to reference genome fasta. |
-
-#### Optional Arguments:
-| Argument | Description |
-|---|---|
-| __-h, --help__ | Show this help message and exit. |
-| __-x STR, --prefix STR__ | Prefix to add to output file names. Default=None |
-| __-b FILENAME, --store_final_bam FILENAME__ | Store the final set of primary contig alignments to the given file name. Default: Do not store the final BAM. |
-| __-m N, --min_qual_score N__ | Minimum mapping quality score to retain an alignment. Default=30 |
-| __-w PATH, --whitelist PATH__ | Path to BED file containing whitelist regions: i.e., the inverse of blacklist regions. Supplying this will have the effect of excluding alignments that fall entirely within blacklist regions. Default=None |
-| __-d, --drop_missing__ | Omit unpatched reference chromosome records from the output if no contigs map to them. Default: Unpatched chromosomes are printed to output unchanged. |
-
-
-## Output
-
-GPatch produces three output files:
-| File | Description |
-|---|---|
-| __patched.fasta__ | The final patched genome. |
-| __contigs.bed__ | Location of contigs in the coordinate frame of the patched genome. |
-| __patches.bed__ | Location of patches in the coordinate frame of the reference genome. |
-
-
-## Citing GPatch
-Please use the following citation if you use this software in your work:
-
-CITATION_HERE

GPatch-0.3.6/src/GPatch.egg-info/SOURCES.txt

@@ -1,17 +0,0 @@
-MANIFEST.in
-README.md
-pyproject.toml
-setup.py
-./CODE_OF_CONDUCT.md
-./LICENSE
-./requirements.txt
-src/GPatch/GPatch.py
-src/GPatch/__init__.py
-src/GPatch/__main__.py
-src/GPatch.egg-info/PKG-INFO
-src/GPatch.egg-info/SOURCES.txt
-src/GPatch.egg-info/dependency_links.txt
-src/GPatch.egg-info/entry_points.txt
-src/GPatch.egg-info/not-zip-safe
-src/GPatch.egg-info/requires.txt
-src/GPatch.egg-info/top_level.txt