Description

Computational framework for tracking and quantifying DNA damage patterns among ancient DNA sequencing reads generated by Next-Generation Sequencing platforms.

Input

name:type
description
pattern

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

bam:file

BAM file

*.{bam}

fasta:file

Fasta file, the reference the input BAM was mapped against

*.{fasta}

Output

name:type
description
pattern

runtime_log

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Runtime_log.txt:file

Log file with a summary of command lines used and timestamps.

Runtime_log.txt

fragmisincorporation_plot

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Fragmisincorporation_plot.pdf:file

A pdf file that displays both fragmentation and misincorporation patterns.

Fragmisincorporation_plot.pdf

length_plot

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Length_plot.pdf:file

A pdf file that displays length distribution of singleton reads per strand and cumulative frequencies of C->T at 5’-end and G->A at 3’-end are also displayed per strand.

Length_plot.pdf

misincorporation

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/misincorporation.txt:file

Contains a table with occurrences for each type of mutations and relative positions from the reads ends.

misincorporation.txt

lgdistribution

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/lgdistribution.txt:file

Contains a table with read length distributions per strand.

lgdistribution.txt

dnacomp

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/dnacomp.txt:file

Contains a table of the reference genome base composition per position, inside reads and adjacent regions.

dnacomp.txt

stats_out_mcmc_hist

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Stats_out_MCMC_hist.pdf:file

A MCMC histogram for the damage parameters and log likelihood.

Stats_out_MCMC_hist.pdf

stats_out_mcmc_iter

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Stats_out_MCMC_iter.csv:file

Values for the damage parameters and log likelihood in each MCMC iteration.

Stats_out_MCMC_iter.csv

stats_out_mcmc_trace

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Stats_out_MCMC_trace.pdf:file

A MCMC trace plot for the damage parameters and log likelihood.

Stats_out_MCMC_trace.pdf

stats_out_mcmc_iter_summ_stat

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Stats_out_MCMC_iter_summ_stat.csv:file

Summary statistics for the damage parameters estimated posterior distributions.

Stats_out_MCMC_iter_summ_stat.csv

stats_out_mcmc_post_pred

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Stats_out_MCMC_post_pred.pdf:file

Empirical misincorporation frequency and posterior predictive intervals from the fitted model.

Stats_out_MCMC_post_pred.pdf

stats_out_mcmc_correct_prob

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/Stats_out_MCMC_correct_prob.csv:file

Position specific probability of a C->T and G->A misincorporation is due to damage.

Stats_out_MCMC_correct_prob.csv

dnacomp_genome

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/dnacomp_genome.csv:file

Contains the global reference genome base composition (computed by seqtk).

dnacomp_genome.csv

rescaled

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/*rescaled.bam:file

Rescaled BAM file, where likely post-mortem damaged bases have downscaled quality scores.

*.{bam}

pctot_freq

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/5pCtoT_freq.txt:file

Contains frequencies of Cytosine to Thymine mutations per position from the 5’-ends.

5pCtoT_freq.txt

pgtoa_freq

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/3pGtoA_freq.txt:file

Contains frequencies of Guanine to Adenine mutations per position from the 3’-ends.

3pGtoA_freq.txt

fasta

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

results_*/*.fasta:file

Alignments in a FASTA file, only if flagged by -d.

*.{fasta}

folder

meta:map

Groovy Map containing sample information e.g. [ id:‘test’, single_end:false ]

*/:directory

Folder created when —plot-only, —rescale and —stats-only flags are passed.

*/

versions

versions.yml:file

File containing software versions

versions.yml