当前位置: 首页 > news >正文

北京南站列车时刻表商务网站模板

北京南站列车时刻表,商务网站模板,怎么用思维导图做网站结构图,扬中网站建设哪家好在进行变异检测时#xff0c;以群体基因组重测序数据为例#xff0c;涉及到的个体基本都是上百个#xff0c;而其中大多数流程均是重复的步骤。 本文将基于GATK进行SNP calling的流程写入循环#xff0c;便于批量分析。 1 涉及变量 1.工作目录work_dir/ 2.参考基因组ref…在进行变异检测时以群体基因组重测序数据为例涉及到的个体基本都是上百个而其中大多数流程均是重复的步骤。 本文将基于GATK进行SNP calling的流程写入循环便于批量分析。 1 涉及变量 1.工作目录work_dir/ 2.参考基因组ref_genome.fa 3.Reads列表read_list.txt 4.测序平台Illumina 5.调用线程数 2 调用数据 1.参考基因组ref_genome.fa 2.重测序数据sample1/sample1_1.fq.gz、sample1/sample1_2.fq.gz…… 3.Reads列表read_list.txt 生成方法预先将存放各个个体Reads的文件夹放入一个文件夹work_dir/然后使用下列命令生成 ls work_dir/ read_list.txt3 主要脚本 usage: bash GATK_pipeline.sh work_dir/ ref_genome.fa read_list.txt Illumina 10GATK_pipeline.sh #---------------------------------------------------------------# # objection defined by user # #---------------------------------------------------------------#set -au# 1. # Master dir.: WORK_dir$1# 2. # Reference genome: REF$2# 3. # Read list: READ_list$3# 4. # Seqencing platform: PL$4# 5. # number of threads: NT$5#---------------------------------------------------------------# # main loop for SNPs calling by gatk pipeline # #---------------------------------------------------------------##READ_list.txt is a list of read groups. while read -r READdoSAMPLESM_${READ} ID${READ} READ1${WORK_dir}${READ}_1.fq READ2${WORK_dir}${READ}_2.fq OUT${READ}#1. #Alignning reads to reference genome by BWA-MEM2-mem, producing a .sam data bwa-mem2 \mem \-M \-t ${NT} \-R RG\tID:${ID}\tSM:${SAMPLE}\tPL:${PL} \${REF} \${READ1} \${READ2} \ ${OUT}.sam#2. #Sorting .sam by gatk-SortSam, producing a .bam data gatk \SortSam \-I ${OUT}.sam \-O ${OUT}.bam \-SO coordinate \-VALIDATION_STRINGENCY LENIENT \-CREATE_INDEX true \-TMP_DIR ./${OUT}tmp.sort #3. #Marking dupulications in .bam by gatk-MarkDuplicates #producing a .dup.bam and .dup.txt data gatk \MarkDuplicates \-I ${OUT}.bam \-O ${OUT}.dup.bam \-M ${OUT}.dup.txt \-REMOVE_DUPLICATES true \-VALIDATION_STRINGENCY LENIENT \-CREATE_INDEX true \-TMP_DIR ${OUT}tmp.dup#4. #QC by samtools-flagstat, producing a .dup.bam.stat data samtools \flagstat \${OUT}.dup.bam \ ${OUT}.dup.bam.stat#5. #Calling SNPs by gatk-HaplotypeCaller, producing a .dup.vcf data gatk \HaplotypeCaller \-R ${REF} \-I ${OUT}.dup.bam \-O ${OUT}.dup.vcfdone $READ_list ##
http://wiki.neutronadmin.com/news/173588/

相关文章:

  • 免费做网站软件下载做模板网站赚钱吗
  • 同声传译公司网站建设东莞好的网站国外站建设价格
  • 做程序开发的网站东莞百度seo在哪
  • 怎么做一个购物平台网站平面排版网站
  • 网站开发公司外包wordpress个人站主题
  • VPS做镜像网站国外做网站公司能赚钱
  • 网站怎么做邮箱广州网站优化价格
  • 什么网站可以做旅行行程单wordpress自定义函数
  • 做网站的傻瓜软件手机网站设计思路
  • 不备案的网站能打开吗网站建设与网络营销
  • 网站平台 生态建设长沙网开亿面做网站多少钱
  • 南涧县城乡建设局网站高清图片免费下载
  • 义乌 网站建设推广营销厦门app网站设计
  • 西安学校网站建设多少钱推广策略是什么意思
  • wordpress图片站主题河北邢台特色美食
  • 游戏攻略网站开发小型电子商务网站规划
  • 网站建设怎样避免犯法网站制作工作室专业公司
  • 可烦网站建设南阳网站优化费用
  • p2p的网站开发织梦模板首页修改
  • 浙江建设职业技术学院塘栖校区网站个人网站域名申请
  • 西安专业网站建设服务好国外优秀的网站建设公司
  • 2018爱情动做网站苏州建站公司认准苏州聚尚网络
  • 旅游网站大全施工企业奖惩制度范本
  • 做个公司展示网站多少钱 后期有什么费用网络游戏陪玩
  • 建设肯德基网站的好处葫芦岛建设网站
  • 怎样手机网站建设正保建设工程网站
  • 提升网站速度电销外包怎么收费
  • 十八把网站做扇子wordpress角色权限管理
  • 聊城建网站哪家好电子印章手机在线制作软件
  • 丰城做网站凡科自助建站平台