标签：help exons 五花八门 bed sample 质量 bamdst 软件 gene

1.bamdst -- a BAM Depth Stat Tool  对bam文件进行统计

安装bamdst  https://github.com/shiquan/bamdst

git clone https://github.com/shiquan/bamdst.git
cd bamdst 
make

make后直接help后，查看用法，有时候可能会提示bamdst不存在，本来已经存在了，还会这样报错。把其路径加入环境变量后，就可以了。

bamdst --help    # ./bamdst --help 

使用bamdst  对bed格式目标区域的捕获统计，查看结果中的region.tsv.gz

bamdst -p region.bed -o ./result in.bam

如果没有时候，可以用CCDS  ，参考：https://www.jianshu.com/p/5a548c8df669

cat CCDS.20160908.txt |perl -alne '{/\[(.*?)\]/;next unless $1;$gene=$F[2];$exons=$1;$exons=~s/\s//g;$exons=~s/-/\t/g;print "$F[0]\t$_\t$gene" foreach split/,/,$exons;}'|sort -u |bedtools sort -i >exon_probe.hg38.gene.bed

ls ../alignment/*.bam|while read id;
do 
file=$(basename $id )
sample=${file%%.*} 
echo $sample
mkdir $sample
~/biosoft/bamdst/bamdst -p exon_probe.GRCh38.gene.bed  -o $sample $id 
done

2.samtools depth  计算每一个位点或者区域的测序深度并在标准显示设备中显示

samtools  还有很多其他用法，可以深入研究

3.用bedtools genomecov统计深度和覆盖度  

bedtools其他用法

标签：help,exons,五花八门,bed,sample,质量,bamdst,软件,gene
来源： https://blog.csdn.net/Cassiel60/article/details/89026250

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。