The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

MatchAlign

An R package for parsing GAF and GFA files, allowing easier interpretation and pangenomic analysis.

Installation

# install.packages("remotes")
remotes::install_github("cromazurek/MatchAlign")

Usage

Parse GFA files

GFA (Graphical Fragment Assembly) files encode the variation captured by pangenome graphs. parse_paths_gfa() extracts the segment path for each accession from P (path) and W (walk) lines.

library(MatchAlign)

accession_paths <- parse_paths_gfa("path/to/example.gfa")
# Returns named list: each accession represented as vector of segment names
accession_paths[["CM294333.1"]]
#> [1] "2+" "3+" "4+" "6+" "7+" "9+" "10+"

Parse GAF files

GAF (Graph Alignment Format) files store alignments of sequences to a pangenome graph. parse_paths_gaf() extracts the segments traversed by each alignment.

alignment_paths <- parse_paths_gaf("path/to/example.gaf")
# Returns a named list: read name -> vector of segment names
alignment_paths[["A01127:121:HCLF5DRX5:2:2203:13579:24680"]]
#> [1] "5001" "5002" "5003" "5004" "5005" "5006" "5007"

Input formats

• GFA: plain-text, tab-delimited. P lines (column 2 = accession name, column 3 = comma-separated segments) and W lines (key = sample#hap_index#seq_id, column 7 = >/<-delimited walk) are supported.
• GAF: plain-text, tab-delimited. Column 1 = query name, column 6 = path string (e.g. >seg1>seg2<seg3).

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.