R包SangerSeqR处理ab1数据

#加载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')