#加载sangerseqR包
library(sangerseqR)
#读入数据
seq = readsangerseq('input.ab1')
#读取碱基数据,0.33指的是将达到主峰0.33的次峰定义为杂合子峰
bc = makeBaseCalls(seq, ratio = 0.33)
#读主峰
primarySeq(seq)
#读次峰
secondarySeq(seq)
#输出可视化图像,也可以选择输出到pdf
chromatogram(bc, showcalls = 'both')
R包SangerSeqR处理ab1数据
猜你喜欢
转载自blog.csdn.net/Cassiel60/article/details/89396259
今日推荐
周排行