Skip to content

Commit 007a98a

Browse files
Merge branch 'master' of https://github.com/algbio/GraphChainer
2 parents f986f37 + babcb40 commit 007a98a

File tree

1 file changed

+15
-3
lines changed

1 file changed

+15
-3
lines changed

README.md

Lines changed: 15 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
# GraphChainer
22

3-
Here we will version the main project consisting on building a sequence to graph aligner based on Co-linear Chaining.
3+
GraphChainer is an accurate aligner of long reads to a variation graph, based on co-linear chaining.
44

55
#### Compilation
66

@@ -19,7 +19,7 @@ Quickstart: `./bin/GraphChainer -t 4 -f reads.fastq -g graph.gfa -a out.gam`
1919
Parameters inherited from GraphAligner:
2020
- `-t` Number of threads.
2121
- `-f` Input reads. Format .fasta / .fastq / .fasta.gz / .fastq.gz. You can input multiple files with `-f file1 -f file2 ...` or `-f file1 file2 ...`.
22-
- `-g` Input graph. Format .gfa / .vg.
22+
- `-g` Input graph. **This graph must be acyclic**. Format .gfa / .vg.
2323
- `-a` Output file name. Format .gam or .json.
2424

2525
Parameters related to colinear chaining:
@@ -28,6 +28,18 @@ Parameters related to colinear chaining:
2828
- `--colinear-split-gap <int>` Default 35. The distance between consecutive fragments. If `--speed` is set, then always `--colinear-split-gap = --speed * --colinear-split-len`.
2929
- `--colinear-gap 1000` Default 10000. When converting an optimal chain of anchors into an alignment path, split the path if the distance between consecutive anchors is greater than this value.
3030

31+
### Constructing an (acyclic) variation graph
32+
33+
Use [vg](https://github.com/vgteam/vg) and run:
34+
35+
`vg construct -t 30 -a -r {ref} -v {vcf} -R 22 -p -m 3000000`
36+
37+
### Citation
38+
39+
If you use GraphChainer, please cite as:
40+
41+
Jun Ma, Manuel Cáceres, Leena Salmela, Veli Mäkinen, Alexandru I. Tomescu. GraphChainer: Co-linear Chaining for Accurate Sequence-to-Variation-Graph Alignment. Submitted, 2021
42+
3143
### Credits
3244

33-
GraphChainer is built on the excellent code base of [GraphAligner](https://github.com/maickrau/GraphAligner), which is released under [MIT License](https://github.com/maickrau/GraphAligner/blob/master/LICENSE.md). GraphAligner is described in the paper [GraphAligner: Rapid and Versatile Sequence-to-Graph Alignment](GraphAligner: rapid and versatile sequence-to-graph alignment) by Mikko Rautiainen and Tobias Marschall.
45+
GraphChainer is built on the excellent code base of [GraphAligner](https://github.com/maickrau/GraphAligner), which is released under [MIT License](https://github.com/maickrau/GraphAligner/blob/master/LICENSE.md). GraphAligner is described in the paper [GraphAligner: Rapid and Versatile Sequence-to-Graph Alignment](https://doi.org/10.1186/s13059-020-02157-2) by Mikko Rautiainen and Tobias Marschall.

0 commit comments

Comments
 (0)