GPatch 0.3.6__tar.gz → 0.3.8__tar.gz

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: GPatch
-Version: 0.3.6
+Version: 0.3.8
 Summary: Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.
 Home-page: https://github.com/adadiehl/GPatch
 Author: Adam Diehl
@@ -61,7 +61,7 @@ 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. |
-        
+        | __-t, --no_trim__ | 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. |
         
         ## Output
         
@@ -76,7 +76,10 @@ Description: # GPatch
         ## 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.3.6 → GPatch-0.3.8}/README.md

@@ -53,7 +53,7 @@ Starting with alignments of contigs to a reference genome, produce a chromosome-
 | __-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. |
-
+| __-t, --no_trim__ | 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. |
 
 ## Output
 
@@ -68,4 +68,6 @@ GPatch produces three output files:
 ## 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

{GPatch-0.3.6 → GPatch-0.3.8}/setup.py

@@ -29,7 +29,7 @@ def main():
 
     metadata = dict(
         name="GPatch",
-        version="0.3.6",
+        version="0.3.8",
         author="Adam Diehl",
         author_email="adadiehl@umich.edu",
         description="Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.",

{GPatch-0.3.6 → GPatch-0.3.8}/src/GPatch/GPatch.py

@@ -7,7 +7,7 @@ from Bio.Seq import Seq
 import pysam
 import random
 
-__version__ = "0.3.6"
+__version__ = "0.3.8"
 
 """
 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
 
 
@@ -171,6 +180,9 @@ 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.')
+    parser.add_argument('-t', '--no_trim',
+                        required=False, default=False, 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.')
     
     args = parser.parse_args()
 
@@ -207,6 +219,7 @@ def main():
     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.
@@ -289,7 +302,15 @@ def main():
                 # 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
+                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)

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: GPatch
-Version: 0.3.6
+Version: 0.3.8
 Summary: Assemble contigs into a chromosome-scalse pseudo-assembly using alignments to a reference sequence.
 Home-page: https://github.com/adadiehl/GPatch
 Author: Adam Diehl
@@ -61,7 +61,7 @@ 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. |
-        
+        | __-t, --no_trim__ | 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. |
         
         ## Output
         
@@ -76,7 +76,10 @@ Description: # GPatch
         ## 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