library(LTRpred)
# de novo LTR transposon prediction of 'O. sativa'
LTRpred(
genome.file = "Osativa.fa",
cluster = TRUE,
cores = 4,
copy.number.est = FALSE,
minlenltr = 100,
maxlenltr = 5000,
mindistltr = 4000,
maxdistltr = 30000,
mintsd = 3,
maxtsd = 20,
vic = 80,
overlaps = "no",
xdrop = 7,
motifmis = 1,
pbsradius = 60,
pbsalilen = c(8,40),
pbsoffset = c(0,10),
quality.filter = TRUE,
n.orfs = 0
)
# import LTRpred output
Osativa_LTRpred <- read.ltrpred("Osativa_ltrpred/Osativa_LTRpred_DataSheet.tsv")
LTR转座子的重新注释
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
推荐阅读更多精彩内容
- =========================================================...