vgl-hub / gfastats

A single fast and exhaustive tool for summary statistics and simultaneous *fa* (fasta, fastq, gfa [.gz]) genome assembly file manipulation.

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Agp integration does not support terminal gaps

gf777 opened this issue · comments

The presence of a terminal gap (gap at the end of scaffolds, e.g. attached file between scaffold25 and 28) will result in incorrect behavior
hap.unlocs.no_hapdups.agp.zip

This issue has been fixed in the latest commit. Terminal gaps in AGP (both at the beginning and end of sequences) are now ignored and a warning about their presence in the AGP is issued. This is consistent with the requirement of avoiding terminal gaps at scaffold ends (e.g. in NCBI submissions).