Skip to content

Commit 6b536e5

Browse files
committed
.
1 parent 3e650b0 commit 6b536e5

1 file changed

Lines changed: 12 additions & 12 deletions

File tree

README.md

Lines changed: 12 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -146,6 +146,17 @@ Example HTML output for one sample:
146146

147147
![HTML Output Example](html_example.png)
148148

149+
## Design and how Virasign works
150+
151+
Virasign works in two steps: **(1)** map reads to a large viral DB (e.g. RVDB, RefSeq), **(2)** pick one best reference per viral species and remap all reads to that curated set.
152+
153+
- **Best-reference selection**: References are deduplicated by viral species and one representative ref per species is chosen (by read count, identity, coverage), giving a small set of best refs instead of many near-duplicates.
154+
- **Why remap**: In one pass to the full DB, minimap2 reports one primary alignment per read, so the best ref “wins” and other species get no count. Remapping to the curated set gives correct counts per virus and meaningful BAM/FASTQ for consensus and downstream analysis, and increases specificity (fewer spurious hits).
155+
- **Why not secondary mapping**: Reporting multiple alignments per read in the first pass would create many false positives with a large DB and lower specificity. Virasign uses primary-only for the first pass, then deduplicates and remaps.
156+
- **Scale and flexibility**: Large viral DBs give sensitivity and completeness (e.g. catching more true viral hits, including low-abundance or divergent viruses). Custom databases and accessions are supported, and the tool is expandable to future reference sets.
157+
158+
Also on [**Zenodo**](https://zenodo.org/records/18387009) and [**Docker**](https://hub.docker.com/repository/docker/daanjansen94/virasign/general).
159+
149160
## Citation
150161

151162
If you use Virasign in your research, please cite:
@@ -164,15 +175,4 @@ Contributions are welcome! Please feel free to submit a Pull Request.
164175

165176
## Support
166177

167-
If you encounter any problems or have questions, please open an issue on GitHub.
168-
169-
## Design and how Virasign works
170-
171-
Virasign works in two steps: **(1)** map reads to a large viral DB (e.g. RVDB, RefSeq), **(2)** pick one best reference per viral species and remap all reads to that curated set.
172-
173-
- **Best-reference selection**: References are deduplicated by viral species and one representative ref per species is chosen (by read count, identity, coverage), giving a small set of best refs instead of many near-duplicates.
174-
- **Why remap**: In one pass to the full DB, minimap2 reports one primary alignment per read, so the best ref “wins” and other species get no count. Remapping to the curated set gives correct counts per virus and meaningful BAM/FASTQ for consensus and downstream analysis, and increases specificity (fewer spurious hits).
175-
- **Why not secondary mapping**: Reporting multiple alignments per read in the first pass would create many false positives with a large DB and lower specificity. Virasign uses primary-only for the first pass, then deduplicates and remaps.
176-
- **Scale and flexibility**: Large viral DBs give sensitivity and completeness (e.g. catching more true viral hits, including low-abundance or divergent viruses). Custom databases and accessions are supported, and the tool is expandable to future reference sets.
177-
178-
Also on [**Zenodo**](https://zenodo.org/records/18387009) and [**Docker**](https://hub.docker.com/repository/docker/daanjansen94/virasign/general).
178+
If you encounter any problems or have questions, please open an issue on GitHub.

0 commit comments

Comments
 (0)