| A simple FDR peak caller. The script is designed for processing DamID-seq datasets such as those generated by the damidseq_pipeline software, but will work on any DNA binding track in bedgraph or GFF format (for example, background-subtracted ChIP-seq data). The output is a GFF file of all peaks with an FDR less than an assigned value (by default, FDR < 0.01). The mean binding intensity of the peak is represented by the score column of the GFF (column 6) and the FDR is present in the attributes column (column 9). The peak file can be loaded into genome browser software such as IGV for viewing, and associated genes can be called with the included peaks2genes script. |
hg clone https://toolshed.g2.bx.psu.edu/repos/mvdbeek/damidseq_findpeaks
| Name | Description | Version | Minimum Galaxy Version |
|---|---|---|---|
| Simple FDR random permutation peak caller | 0.1.3 | 16.01 | |