sqz is a tool for compressing and decompressing path annotations in GFA files. It uses a modified version of the GFA-format containing Q- and Z-lines.
Use
sqz compress <GFA-FILE> > <COMPRESSED-GFA-FILE>
to compress a GFA file and
sqz decompress <COMPRESSED-GFA-FILE> > <GFA-FILE>
sqz is written in RUST and requires a working RUST build system (version >= 1.74.1) for installation. See here for more details.
git clone git@github.com:codialab/sqz.git
cd sqz
cargo build --release
The format of sqz is mostly based on the GFA format, but extended to include two new types of lines: Q- and Z-lines.
A Q-line defines part of a compressed walk that can be used as part of other compressed walks.
| Column | Field | Type | Regexp | Description |
|---|---|---|---|---|
| 1 | RecordType |
Character | Q |
Record type |
| 2 | Name |
String | [!-)+-<>-~][!-~]* |
Rule name |
| 3 | CompressedWalk |
String | ([><][!-;=?-~]+)+ |
Compressed Walk |
A Walk is defined as
<walk> ::= ( `>' | `<' <segId> )+where <segId> corresponds either to the identifier of a segment or the
identifier oft a Q-line. A valid walk must exist in the graph.
A walk line describes an oriented walk in the graph. It is only intended for a graph without overlaps between segments. Note that Z-lines can not use jump connections (introduced in v1.2).
| Column | Field | Type | Regexp | Description |
|---|---|---|---|---|
| 1 | RecordType |
Character | W |
Record type |
| 2 | SampleId |
String | [!-)+-<>-~][!-~]* |
Sample identifier |
| 3 | HapIndex |
Integer | [0-9]+ |
Haplotype index |
| 4 | SeqId |
String | [!-)+-<>-~][!-~]* |
Sequence identifier |
| 5 | SeqStart |
Integer | \*|[0-9]+ |
Optional Start position |
| 6 | SeqEnd |
Integer | \*|[0-9]+ |
Optional End position (BED-like half-close-half-open) |
| 7 | CompressedWalk |
String | ([><][!-;=?-~]+)+ |
Compressed Walk |
For a haploid sample, HapIndex takes 0. For a diploid or polyploid sample,
HapIndex starts with 1. For two W-lines with the same
(SampleId,HapIndex,SeqId), their [SeqSart,SeqEnd) should have no
overlaps. A Walk is defined as
<walk> ::= ( `>' | `<' <segId> )+where <segId> corresponds either to the identifier of a segment or the
identifier oft a Q-line. A valid walk must exist in the graph.
S s11 ACCTT
S s12 TC
S s13 GATT
L s11 + s12 - 0M
L s12 - s13 + 0M
L s11 + s13 + 0M
Q q1 >s11<s12
Z NA12878 1 chr1 0 11 >q1>s13