Skip to content

Commit 92b7d57

Browse files
authored
Update README.md
1 parent 100aead commit 92b7d57

File tree

1 file changed

+4
-2
lines changed

1 file changed

+4
-2
lines changed

README.md

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -30,8 +30,10 @@ cargo build --release
3030
```
3131
### (Optional) Prune the Himito graph using paired srWGS data
3232
The graph served as the foundation for downstream assembly, variant calling and methylation analysis.
33-
If you have paired short reads data, you can optional refine the graph by using ```Himito correct```.
34-
(Optional) trim graph paths using short reads, short reads aligned to chrM should be firstly compressed into a [msBWT](https://github.com/HudsonAlpha/rust-msbwt).
33+
34+
If you have paired short reads data, you can refine the graph by using ```Himito correct```.
35+
36+
```Himito correct``` trims graph paths according to the occurrence counts of path kmers in the paired srWGS data. We recommend to first subset the srWGS BAM file to reads aligned to chrM, then compressed these reads into a [msBWT](https://github.com/HudsonAlpha/rust-msbwt).
3537
```
3638
msbwt2-build -o sr_msbwt.npy <srWGS.chrM.fasta.gz>
3739
./target/release/Himito correct -g <output.gfa> -b <bwt_file, e.g. sr_msbwt.npy> -o <corrected.gfa> -m <minimal_supporting_sr> -q <query_length, should be less than short read length>

0 commit comments

Comments
 (0)