1. Shifu Chen, Yanqing Zhou, Yaru Chen, Jia Gu; fastp: an ultra-fast all-in-one FASTQ preprocessor, Bioinformatics, Volume 34, Issue 17, 1 September 2018, Pages i884–i890
    2. Wang, K., Liang, C., Liu, J. et al. Prediction of piRNAs using transposon interaction and a support vector machine. BMC Bioinformatics 15, 419 (2014). https://doi.org/10.1186/s12859-014-0419-6
    3. novoalign [http://www.novocraft.com/products/novoalign/]
    4. Sai Lakshmi, S. and S. Agrawal, piRNABank: a web resource on classified and clustered Piwi-interacting RNAs. Nucleic Acids Res, 2008. 36(Database issue): p. D173-7.
    5. Love M.I., Huber W., Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. doi: 10.1186/s13059-014-0550-8.
    6. Enright AJ, John B, Gaul U, Tuschl T, Sander C and Marks DS; (2003) Genome Biology; 5(1):R1.
    7. Huang HY, Lin YC, Li J, Huang KY, Shrestha S, Hong HC, Tang Y, Chen YG, Jin CN, Yu Y, Xu JT, Li YM, Cai XX, Zhou ZY, Chen XH, Pei YY, Hu L, Su JJ, Cui SD, Wang F, Xie YY, Ding SY, Luo MF, Chou CH, Chang NW, Chen KW, Cheng YH, Wan XH, Hsu WL, Lee TY, Wei FX, Huang HD*. (2020) Nucleic Acids Research.
    8. Guangchuang Yu, Li-Gen Wang, Yanyan Han, Qing-Yu He. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS: A Journal of Integrative Biology. 2012, 16(5):284-287.
    9. https://cran.r-project.org/web/packages/pheatmap/