samtools过滤多重比对reads的命令如下：
samtools view -h -q 1 -F 4 -F 256 big.sort.dedup.bam |grep -v XA:Z |grep -v SA:Z |samtools view -b - > big.sort.dedup.unique.bam
其中参数：
-q 1：过滤掉比对质量小于1的reads；
-F 4：过滤掉没有比对上的reads；
-F 256：过滤掉比对上的多次的reads。
此外，我们在网易云课堂上有各种教学视频，有兴趣可以了解一下：
1. 文章越来越难发？是你没发现新思路，基因家族分析发2-4分文章简单快速，学习链接：基因家族分析实操课程
2. 转录组数据理解不深入？图表看不懂？点击链接学习深入解读数据结果文件，学习链接：转录组（有参）结果解读；转录组（无参）结果解读
3. 转录组数据深入挖掘技能-WGCNA，提升你的文章档次，学习链接：WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘？学习链接：转录组标准分析后的数据挖掘
5. 微生物16S/ITS/18S分析原理及结果解读
6. 更多学习内容：linux、perl、R语言画图，更多免费课程请点击以下链接：
https://study.omicsclass.com/
