famCircle 基因家族分析，tandem可视化分析

分析流程

基因组下载以及处理 路线： 前期处理序列比对、共线性扫描、ks计算 ks展示 全基因组 ks展示 block ks展示 circle全基因组共线性展示 line局部共线性绘制 hmmer 隐马尔可夫模型搜索筛选 screen tandem数量分部 断开处： pfam结构域筛选# http://pfam.xfam.org/ α/β筛选# http://cb.csail.mit.edu/cb/paircoil2/ 主程序选择： inner typing 数据格式化 inner绘图 outer typing 数据格式化 outer绘图 part局部tandem绘制

配置文件解释

总体ks分布可视化

run ks ks = ks file vertical = False bins = 100 model = YN00/NG86 savefile = save file (*.png, *.pdf)

基因块ks分布可视化

run Ks_allocation ks = ks file area = 0,2 vertical = False model = YN00/NG86 blockfile = block file blocklength = 6 pvalue = 0.05 savefile = save file (*.png, *.pdf)

共线性可视化

run circle lens = lens file gff = gff file chrolist = Genome name ks = ks file genepairs = genepairs file peripheral = False block = 6 Ks_concern = 0,0.15 bridge = 1 radius = 0.3 savefile = save file (*.png, *.pdf)

基因家族查找

run hmmer pep = pep file# 蛋白质文件 hmmmoldpath = hmm file# 模型地址 format_conversion = Fales# 格式转换 comparison = muscle/clustal# 比对软件 e_value1 = value1# 初筛阈值 e_value2 = value2# 复筛阈值

结构域筛选

α/β筛选

结构域分布情况

run screen domainpath = domain file path lens = lens file gff = gff file chrolist = Genome name series = 25# 串联数 outpath = out file path

文件格式化

run typing domainpath = domain file# 文件路径 domainlist = Genome name# 文件名列表 position = inner/outer# 目标程序格式 savefile = out file# 保存文件

内卷型的tandem可视化

run inner lens = lens file gff = gff file chrolist = Genome name ks = ks file genefamily = famliy file Ks_concern = 0.1,0.2,1.4,1.5# ks分割参数可调 bridge = 0.05# 最远连接 radius = 0.3 peripheral = False savecsv = outer file (.csv) savefile = save file (.png, *.pdf)

放射型的tandem可视化

run outer lens = lens file gff = gff file chrolist = Genome name ks = ks file genepairs = genepairs file Ks_concern = 0,0.15# ks分割参数可调 series = 25# 串联数 clusters = None# 集团参数可调 peripheral = False bridge = 1# 跨区域连线参数可调 radius = 0.3 savecsv = outer file (.csv) savefile = save file (.png, *.pdf)

[circle] lens = lens file gff = gff file species_list = Genome name ks = ks file genepairs = genepairs file block = 6 radius = 0.45 savefile = save file (*.png, *.pdf)

[circle_all] lens = lens file gff = gff file species_list = Genome name blockfile = block file radius = 0.3 savefile = save file (*.png, *.pdf)

[hmmer] pep = pep file cds = cds file hmmmoldpath = hmm file format_conversion = Fales comparison = clustal e_value1 = value1 e_value2 = value2

[inner] lens = lens file gff = gff file chrolist = Genome name ks = ks file genefamily = famliy file Ks_concern = 0,1.5 radius = 0.3 space = 0.005 clusters = True peripheral = False savecsv = outer file (.csv) savefile = save file (.png, *.pdf)

[Ks] cds_file = cds file pep_file = pep file align_software = muscle pairs_file = gene pairs file ks_file = ks result

[Ks_allocation] ks = ks file species_list = specise name area = 0,2.5 model = YN00/NG86 savefile = save file (*.png, *.pdf)

[Ks_block] species_list = specise name ks = ks file area = 0,2.5 model = YN00/NG86 blockfile = block file blocklength = 6 pvalue = 1e-5 savecsv = save csv savefile = save file (*.png, *.pdf)

[line] pairs_file = pairs file gff1 = gff1 file gff2 = gff2 file lens1 = lens1 file lens2 = lens2 file chr1_name = chr1 name chr2_name = chr2 name savefile = savefile(.png,.pdf)

[outer] lens = lens file gff = gff file chrolist = Genome name ks = ks file genefamily = famliy file Ks_concern = 0,0.15 radius = 0.3 space = 0.005 clusters = True peripheral = False savecsv = outer file (.csv) savefile = save file (.png, *.pdf)

[part] lens = lens file gff = gff file chrolist = Genome name ks = ks file genefamily = famliy file Ks_concern = 0,1.5 radius = 0.3 space = 0.005 clusters = True peripheral = False savecsv = outer file (.csv) savefile = save file (.png, *.pdf)

[screen] domainpath = domain file lens = lens file gff = gff file chrolist = Genome name series = 25 outpath = out file

[typing] domainpath = domain file domainlist = Genome name savefile = out file