Skip to content

Commit a2db4cb

Browse files
Update authentication.qmd
1 parent 23d8f9b commit a2db4cb

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

authentication.qmd

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -396,7 +396,7 @@ Typically, reads with PMD scores greater than 3 are considered to be reliably an
396396
As mapDamage, PMDtools can also compute deamination profile. However, the advantage of PMDtools that it can compute deamination profile for UDG / USER treated samples (with the flag *--CpG*). For this purpose, PMDtools uses only CpG sites which escape the treatment, so deamination is not gone completely and there is a chance to authenticate treated samples. Computing deamination pattern with PMDtools can be achieved with the following command line (please note that the scripts *pmdtools.0.60.py* and *plotPMD.v2.R* can be downloaded from the github repository here https://github.com/pontussk/PMDtools):
397397

398398
```bash
399-
samtools view Y.pestis_sample10.bam | pmdtools --platypus > PMDtemp.txt
399+
samtools view Y.pestis_sample10.bam | pmdtools --platypus > PMD_temp.txt
400400
```
401401

402402
We can then run simple R commands directly from the terminal (without loading R itself with) the following.

0 commit comments

Comments
 (0)