该问题已被锁定!
2
关注
1667
浏览

如何提取比对上某一个基因或者某一段序列的所有reads

为什么被折叠? 0 个回复被折叠
孟浩巍 超级管理员 用户来自于: 北京市
2018-10-05 13:41
samtools view 命令的使用说明如下: [code]Usage: samtools view [options] || [region ...] Options: -b output BAM -C output CRAM (requires -T) -1 use fast BAM compression (implies -b) -u uncompressed BAM output (implies -b) -h include header in SAM output -H print SAM header only (no alignments) -c print only the count of matching records -o FILE output file name [stdout] -U FILE output reads not selected by filters to FILE [null] -t FILE FILE listing reference names and lengths (see long help) [null] -L FILE only include reads overlapping this BED FILE [null] -r STR only include reads in read group STR [null] -R FILE only include reads with read group listed in FILE [null] -q INT only include reads with mapping quality >= INT [0] -l STR only include reads in library STR [null] -m INT only include reads with number of CIGAR operations consuming query sequence >= INT [0] -f INT only include reads with all of the FLAGs in INT present [0] -F INT only include reads with none of the FLAGS in INT present [0] -G INT only EXCLUDE reads with all of the FLAGs in INT present [0] -s FLOAT subsample reads (given INT.FRAC option value, 0.FRAC is the fraction of templates/read pairs to keep; INT part sets seed) -M use the multi-region iterator (increases the speed, removes duplicates and outputs the reads as they are ordered in the file) -x STR read tag to strip (repeatable) [null] -B collapse the backward CIGAR operation -? print long help, including note about region specification -S ignored (input format is auto-detected) --input-fmt-option OPT[=VAL] Specify a single input file format option in the form of OPTION or OPTION=VALUE -O, --output-fmt FORMAT[,OPT[=VAL]]... Specify output format (SAM, BAM, CRAM) --output-fmt-option OPT[=VAL] Specify a single output file format option in the form of OPTION or OPTION=VALUE -T, --reference FILE Reference sequence FASTA FILE [null] -@, --threads INT Number of additional threads to use [0][/code]其中,-L参数可以帮助你完成你的要求,你只需要把你需要的区间存成BED格式就好。  

关于作者

Czc 注册会员

这家伙很懒,还没有设置简介

问题动态

发布时间
2018-10-04 10:32
更新时间
2018-10-05 13:41
关注人数
2 人关注

相关问题

如何批量绘制多条折线
seurat包,如何调节大热图左侧的基因字体大小
从多个.fa文件中提取以“poptri.”开头的蛋白序列
有台Dell工作站,如何配置成服务器
如何方便快捷地批量重命名qiime2做出的ASV号?
x <- x[keep.exprs,, keep.lib.sizes=FALSE] 请问一下这条命令该如何解读呢?
如何获得cas9蛋白的高效低效的sgRNA
SRR1924430的sra数据与hdf5提取的fastq不一致
如何提取可变剪切位点?
求助一个shell脚本问题,如何批量处理下面这种情况?

推荐内容

群体内同源基因的所有变异
VCONTACT2的结果文件genome_by_genome_overview如何进行统计分类
去除批次效应选定靶标细胞群后如何进行亚群分析?
ASR祖先序列重建,最后一步使用PAML时出现一些问题
使用Tracking Tumor ImmunoPhenotype(TIP)网站分析TCGA的BLCA_tpm数据
infercnv运行报错
如何提取可变剪切位点?
使用tophat2和bowtie1寻找环形RNA时报错
如何对特征数量少的空间蛋白组数据进行细胞聚类?
如何按一个列表对基因型文件进行过滤,剔除不需要的样本?
All Rights Reserved Powered BY WeCenter V4.1.0 © 2025