I am pleased to share with you the R package I have developed for analysis of ctDNA sequencing data (ctDNAtools).
You can do two kinds of analyses with ctDNAtools:
- Minimal residual disease testing: this is useful if you have a list of pre-detected mutations in the tumor (for example detected at diagnosis), and you want to test a follow-up sample whether it has cancer (i.e. ctDNA) or not.
- Analysis of cell-free DNA fragmentation: you can study the fragment size histograms, profiles and fragmentation patterns along the genome.
The ctDNAtools package is mainly built on GenomicAlignments and Rsamtools packages, and most functions require inputs including one or more of a bam file, mutations data frame, targets data frame and a reference genome.
You can check the full documentation and the examples in the vignette, and I am happy to answer questions about the usage.