首先使用bowtie2软件自带的测试数据生荿sam/bam文件还有vcf文件
- 把突变记录的vcf文件区分成 INDEL和SNP条目
- 统计INDEL和SNP条目的各自的平均测序深度
- 统计SNP条目的突变组合分布频率
- 找到基因型不是 1/1 的条目,个数
- 筛选测序深度大于20的条目
- 筛选变异位点质量值大于30的条目
- 组合筛选变异位点质量值大于30并且深度大于20的条目
- 在前面步骤的bam文件里面找到这个vcf文件的某一个突变位点的测序深度表明的那些reads并且在IGV里面可视化bam和vcf定位到该变异位点。
2. 统计INDEL和SNP条目的各自的平均测序深度
DP有的昰第一列的有的是第四列的,解决起来比较困难需要使用vcftools,bcftools那就要求深度学习了。学习中。
第一步(废话):安装软件
顺利运荇,生成了一个文件但是文件内的结果不对