# #——消息= FALSE ---------------------------------------------------------- 库(SICtools ) ## ----------------------------------------------------------------------------- bam1 < -执行(包=‘SICtools’,‘extdata’,‘example1.bam’)bam2 < -执行(包=‘SICtools’,‘extdata’,‘example2.bam’)refFsa < -执行(包=‘SICtools’,‘extdata’,‘example.ref.fasta’)snpDiffDf < - snpDiff (bam1、bam2 refFsa, chr04, 962501, 1026983, pValueCutOff = 1, baseDistCutOff = 0) snpDiffDf # #----------------------------------------------------------------------------- snpDiffDfSort < - snpDiffDf[秩序(snpDiffDf p.value美元,snpDiffDf d.value美元),]snpDiffDfSort ## ----------------------------------------------------------------------------- indelDiffDf < - indelDiff (bam1、bam2 refFsa,“chr07”,828514年,828914年,pValueCutOff = 1, gtDistCutOff = 0) indelDiffDfSort < - indelDiffDf[秩序(indelDiffDf p.value美元,indelDiffDf d.value美元)]indelDiffDfSort # #----------------------------------------------------------------------------- sessionInfo ()