#安装TEsorter和seqkit
conda install -c bioconda tesorter seqkit
#经笔者测试TEsorter貌似只对LTR分类较好,对DNA转座子分类比较一般,从转座子库中提取LTR和未分类的TE
cat AST.fa.mod.EDTA.TElib.fa | seqkit grep --by-name --use-regexp --ignore-case --pattern Unknown --pattern LTR > AST.LTRandUnknown.fa
#除了LTR和未分类的TE,剩下的另存为一个文件
cat AST.fa.mod.EDTA.TElib.fa | seqkit grep -v --by-name --use-regexp --ignore-case --pattern Unknown --pattern LTR > AST.other.fa
#分类
TEsorter -db rexdb-plant -p 10 AST.LTRandUnknown.fa
#分类完了可以再用RepeatMasker比对回基因组，统计各重复序列比例
RepeatMasker -pa 40 -lib hap2.after_sorter.fa -no_is -xsmall -poly -html -gff -dir masker2 hap2_rename.fasta
perl /opt/conda/share/RepeatMasker/util/buildSummary.pl hap2_rename.fasta.out > repeatSummary.txt
#生成gff文件
grep TE *out | awk -F" *" '{print $6"\tEDTA\t"$12"\t"$7"\t"$8"\t.\t"$10"\t.\tID=TE_"$16";Name="$11";Classification="$12}' > TEsorter.gff
sed 's/\tC\t/\t-\t/' TEsorter.gff > TEsorter.MLA_h2.gff
使用TEsorter对转座子分类
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- transposon属于二型转座子，使用“剪切-粘贴”机制。在转座酶（Transposase）的作用下，将DN...
- 基因组中重复序列大体分为两类：串联重复（Tandem repeats，Tandem Duplication） （T...
