文章目录

  • 关于宏基因组分析
    • 什么是宏基因组分析测序?
  • 准备-----宏基因组集群上需要用到的环境和工具(在服务器对测序数据进行组装和质控,排除宿主的DNA,根据OTU获取菌落的丰度的表,或者进行基因功能注释。)
    • Anaconda------类似于虚拟机,能够快速构建环境,傻瓜式--即conda命令
    • snakemake (v5.14.0)----流程控制软件,调控软件合作运行
    • python3 (v3.6.10), pandas-----生物信息分析软件底层支持
    • bowtie2 (v2.3.5.1)----序列比对软件--主要用于去除宿主的序列
    • samtools (v1.9)----一款对于bam格式文件进行赛选,排序,操作的软件,可以转换文件类型。为metaphlan3做好数据准备
    • metaphlan3 (v3.0)---根据质控数据获取物种丰度
    • seqkit (v0.12.1)
  • **==__关于集群__==**
    • 关于节点篇
        • qstat命令内容展示
        • tail命令
        • awk命令
      • xargs
      • sort 排序命令
      • tr 替换命令
      • 考虑到服务器存在可能有无下载节点,即是否有网络节点,如果需要网络则要跳转至下载节点
    • 投递计算任务
    • 查看文件大小
      • du命令查看看要处理恶的文件的大小
    • 下载任务---scp或是使用xftp其他软件
  • 微生物宏基因组分析
    • ==profile==是什么?
    • 微生物宏基因组分析数据的获取
    • **实际操作**
      • 构建snakemake流程进行生物信息分析
    • 构建流程----非snakemake 的并运行,模仿snakemake的运行模式
    • 通常来说会使用多种方法对数据进行处理,得到不同的结果,通过对比几个方案的不同数据,来获取最佳的数据进行后续处理同时也是一个证明数据,为你自己所选择的方法进行作证。
  • 总结

宏基因组分析第一次分析,我也想成为大佬加油,今年实习一定很傻
今天我又意外听到了hadoops和sge框架的处理系统概念,暂时还没有去弄不过这个snakemake是sge框架下使用的一种调节的工具就可以控制资源这样就可以用同样的资源跑几个任务。

关于宏基因组分析

什么是宏基因组分析测序?

以特定的环境中的整个微生物群落,作为研究对象,不需要对环境中微生物进行分离培养利用高通量测序平台进行测序,可进行微生物的组成及种类进行鉴定,系统功能注释,样品间的物种或基因培养以及物种间的代谢网络分析。
##宏基因组技术路线
样品—提取基因组DNA—基因组DNA片段话—illumina测序–数据质控–高质量数据—生物信息分析
注意宏基因组测序特点为测序不进行PCR,PCR扩增容易引起碱基偏向性,使得菌落的物种DNA比例无法反应真实的环境的变换,所以很费钱。

准备-----宏基因组集群上需要用到的环境和工具(在服务器对测序数据进行组装和质控,排除宿主的DNA,根据OTU获取菌落的丰度的表,或者进行基因功能注释。)

Anaconda------类似于虚拟机,能够快速构建环境,傻瓜式–即conda命令

snakemake (v5.14.0)----流程控制软件,调控软件合作运行

python3 (v3.6.10), pandas-----生物信息分析软件底层支持

bowtie2 (v2.3.5.1)----序列比对软件–主要用于去除宿主的序列


Bowtie2 是将测序reads与长参考序列比对工具 (适用于将长度大约为50到100或1000字符的reads与相对较长的基因组, 如哺乳动物,进行比对)。
通常是比较基因组学(包括识别变体(variation calling),ChIP-seq,RNA-seq,BS-seq)管道的第一步。
可以处理非常长的读数(即10s或100s的千字节),但它针对近期测序仪产生的读数长度和误差模式进行了优化,如Illumina HiSeq 2000,Roche 454和Ion Torrent仪器。
Bowtie2使用FM索引(基于Burrows-Wheeler Transform 或 BWT)对基因组进行索引,以此来保持其占用较小内存。对于人类基因组来说,内存占用在3.2G左右。Bowtie2 支持间隔,局部和双端对齐模式。可以同时使用多个处理器来极大的提升比对速度。

https://blog.csdn.net/u011262253/article/details/79833969


samtools (v1.9)----一款对于bam格式文件进行赛选,排序,操作的软件,可以转换文件类型。为metaphlan3做好数据准备

将sam文件转换成bam文件;然后对bam文件进行各种操作,比如数据的排序(不属于本命令的功能)和提取(这些操作是对bam文件进行的,因而当输入为sam文件的时候,不能进行该操作);最后将排序或提取得到的数据输出为bam或sam(默认的)格式。
https://blog.csdn.net/weixin_30457465/article/details/96721961

metaphlan3 (v3.0)—根据质控数据获取物种丰度

MetaPhlAn 是一种计算工具,用于分析微生物群落(细菌、阿奇亚和真核生物)的组成,这些微生物群落的组成来自使用 StrainPhlAn 的元基因组猎枪测序数据,具有物种水平和应变级分辨率。
https://www.ershicimi.com/p/09b30cf1460901b17a153205d0573a94

seqkit (v0.12.1)

FASTA和FASTQ是用于存储核苷酸和蛋白质序列的基本且普遍存在的格式。FASTA / Q文件的常见操作包括转换,搜索,过滤,重复数据删除,拆分,混排和采样。现有工具仅实现这些操作中的某些操作,而并非特别有效,并且某些仅适用于某些操作系统。此外,所需软件包和运行环境的复杂安装过程会使这些程序的用户友好性降低。
该项目描述了用于FASTA / Q处理的跨平台超快速综合工具包。SeqKit为所有主要操作系统(包括Windows,Linux和Mac OS X)提供可执行的二进制文件,并且可以直接使用而无需任何依赖项或预先配置。与类似工具相比,SeqKit展示了在执行时间和内存使用方面的竞争性能。SeqKit的效率和可用性使研究人员能够快速完成常见的FASTA / Q文件操作。
https://github.com/shenwei356/seqkit

关于集群

将测序出来的数据构建一个prifile

关于节点篇

一般的服务器集群都会分为登录节点和计算节点,为了能够选择到计算节点首先了解节点的情况
使用SSH切换到计算节点

ssh cngb-9 #cngb-9为一个计算节点,切换之后用户名@后面的部分会发生改变
#如果不知道节点怎么办,可以使用qstat来掉看节点状态
qstat -u \*  | tail -n +3 | awk -F'[ @]+' '$5=="r" && $8=="st.q" {print $9}' | awk -F'.' '{print $1}' | sort -V | uniq | tr '\n' ' ' | xargs -I xxx echo qhost -h xxx | xargs -I xxx bash -c xxx
#上面的qstat为查看节点状态,后面通过使用|通道进行优化。

qstat命令内容展示


HOSTNAME------节点名
ARCH ----- 支持处理器架构
NCPU ----- 已经利用的CPU核心数
NSOC ----
NCOR -----可利用的核心总数
NTHR-------线程数
LOAD -------当前负载,负载和核心数一致,100%负载为112
MEMTOT----- 总内存
MEMUSE -----当前使用内存
SWAPTO
SWAPUS
其他的可能还要查下,自己还不太懂

tail命令

-f 循环读取
-q 不显示处理信息
-v 显示详细的处理信息
-c<数目> 显示的字节数
-n<行数> 显示文件的尾部 n 行内容
–pid=PID 与-f合用,表示在进程ID,PID死掉之后结束
-q, --quiet, --silent 从不输出给出文件名的首部
-s, --sleep-interval=S 与-f合用,表示在每次反复的间隔休眠S秒

awk命令

-F fs or --field-separator fs
指定输入文件折分隔符,fs是一个字符串或者是一个正则表达式,如-F:

awk -F'[ @]+' '$5=="r" && $8=="st.q" {print $9}'
#$n 中的n为提取第几列, -F是分隔符号的一个参数后面接分割符号,这里是@和 ' ',这个代码意思是将前面的数据以’ ‘和@分割成不同列,然后匹配出第五列等于'r'和第8列等于'st,q'的行的第九列输出出来。

xargs

是给命令传递参数的一个过滤器,也是组合多个命令的一个工具。可以将管道或标准输入(stdin)数据转换成命令行参数,也能够从文件的输中读取数据。也可以将单行或多行文本输入转换为其他格式,例如多行变单行,单行变多行。默认的命令是 echo,这意味着通过管道传递给 xargs 的输入将会包含换行和空白,不过通过 xargs 的处理,换行和空白将被空格取代。 是一个强有力的命令,它能够捕获一个命令的输出,然后传递给另外一个命令。
-a file 从文件中读入作为sdtin
-e flag ,注意有的时候可能会是-E,flag必须是一个以空格分隔的标志,当xargs分析到含有flag这个标志的时候就停止。
-p 当每次执行一个argument的时候询问一次用户。
-n num 后面加次数,表示命令在执行的时候一次用的argument的个数,默认是用所有的。
-t 表示先打印命令,然后再执行。
-i 或者是-I,这得看linux支持了,将xargs的每项名称,一般是一行一行赋值给 {},可以用 {} 代替。
-r no-run-if-empty 当xargs的输入为空的时候则停止xargs,不用再去执行了。
-s num 命令行的最大字符数,指的是 xargs 后面那个命令的最大命令行字符数。
-L num 从标准输入一次读取 num 行送给 command 命令。
-l 同 -L。
-d delim 分隔符,默认的xargs分隔符是回车,argument的分隔符是空格,这里修改的是xargs的分隔符。
-x exit的意思,主要是配合-s使用。。
-P 修改最大的进程数,默认是1,为0时候为as many as it can ,这个例子我没有想到,应该平时都用不到的吧。

sort 排序命令

 -V, --version-sort            在文本内进行自然版本排序

tr 替换命令

   将数据中的符号替换成你想要的符号

考虑到服务器存在可能有无下载节点,即是否有网络节点,如果需要网络则要跳转至下载节点

ssh software-install   # 只是切换节点而已每个可能名字不一样

投递计算任务

使用前记得调至计算节点
使用qsub提交出来任务

qsub -clear -cwd -l vf=4g,num_proc=4 -P P_项目编号 -binding linear:4  -q st.q work.sh

-cwd #指定当前路径为工作目录,sge的日志会输出到当前路径。
-S #指定远程计算节点的shell路径
-l #指定资源请求,多个请求用逗号(,)隔开
vf=1.5G #任务的预估内存,内存估计的值应稍微大于真实的内存,内存预估偏小可能会导致节点跑挂。
-h=compute-0-15 #指定任务跑在compute-0-15节点上
-p=8 #指定要申请的CPU核心数
-q #指定要投递到的队列,如果不指定的话,SGE会在用户可使用的队列中选择一个。
-P #参数指明任务所属的项目
-p #设置优先级,优先级高的优先执行
-N #指定任务名称
-o #指定标准输出路径
-e #指定标准错误路径
-binding #指定核心数

查看文件大小

du命令查看看要处理恶的文件的大小

```du -sh# 查看当前目录的大小du -sh #*查看当前目录下文件大小```

通过查看文件大小和数据的多少,来调节提交任务的时的参数

下载任务—scp或是使用xftp其他软件

cp   root@192.168.0.101:~/filenotdir  ~   #从远端传送filenotdir文件到本地~ scp   root@192.168.0.101:/libopen*   /usr  #从远端传送以libopen开头的多个文件到本地 /usr scp   -r root@192.168.0.101:/ /usr/   #从远端传送文件夹到本地/usr ,注意-r 在本地命令帮助里没有!

                                WDNM,NMSL,CNM,CCC,NTMGLZGUN

微生物宏基因组分析

profile是什么?

对于profile见解,profile是一类含有特征信息的数据集合,序列数据谱数据-来源于序列多重比对,包含已知关系的序列数据,位置位置特异性的份表以及空位罚分。谱数据中的每一个位置包括了所有可能氨基酸的得分,以及一个空位的罚分和另一个在特异位点维持空位罚分,通过优化过程来建立一个起始于给定多重比对的谱数据,人们试图增加数据的敏感性。
[我的理解就是出来的一些可以用来构建族谱的关键数据。]

微生物宏基因组分析数据的获取

流程解析:

#mermaid-svg-QXcmkdMWppRyBRts .label{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);fill:#333;color:#333}#mermaid-svg-QXcmkdMWppRyBRts .label text{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .node rect,#mermaid-svg-QXcmkdMWppRyBRts .node circle,#mermaid-svg-QXcmkdMWppRyBRts .node ellipse,#mermaid-svg-QXcmkdMWppRyBRts .node polygon,#mermaid-svg-QXcmkdMWppRyBRts .node path{fill:#ECECFF;stroke:#9370db;stroke-width:1px}#mermaid-svg-QXcmkdMWppRyBRts .node .label{text-align:center;fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .node.clickable{cursor:pointer}#mermaid-svg-QXcmkdMWppRyBRts .arrowheadPath{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .edgePath .path{stroke:#333;stroke-width:1.5px}#mermaid-svg-QXcmkdMWppRyBRts .flowchart-link{stroke:#333;fill:none}#mermaid-svg-QXcmkdMWppRyBRts .edgeLabel{background-color:#e8e8e8;text-align:center}#mermaid-svg-QXcmkdMWppRyBRts .edgeLabel rect{opacity:0.9}#mermaid-svg-QXcmkdMWppRyBRts .edgeLabel span{color:#333}#mermaid-svg-QXcmkdMWppRyBRts .cluster rect{fill:#ffffde;stroke:#aa3;stroke-width:1px}#mermaid-svg-QXcmkdMWppRyBRts .cluster text{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts div.mermaidTooltip{position:absolute;text-align:center;max-width:200px;padding:2px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);font-size:12px;background:#ffffde;border:1px solid #aa3;border-radius:2px;pointer-events:none;z-index:100}#mermaid-svg-QXcmkdMWppRyBRts .actor{stroke:#ccf;fill:#ECECFF}#mermaid-svg-QXcmkdMWppRyBRts text.actor>tspan{fill:#000;stroke:none}#mermaid-svg-QXcmkdMWppRyBRts .actor-line{stroke:grey}#mermaid-svg-QXcmkdMWppRyBRts .messageLine0{stroke-width:1.5;stroke-dasharray:none;stroke:#333}#mermaid-svg-QXcmkdMWppRyBRts .messageLine1{stroke-width:1.5;stroke-dasharray:2, 2;stroke:#333}#mermaid-svg-QXcmkdMWppRyBRts #arrowhead path{fill:#333;stroke:#333}#mermaid-svg-QXcmkdMWppRyBRts .sequenceNumber{fill:#fff}#mermaid-svg-QXcmkdMWppRyBRts #sequencenumber{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts #crosshead path{fill:#333;stroke:#333}#mermaid-svg-QXcmkdMWppRyBRts .messageText{fill:#333;stroke:#333}#mermaid-svg-QXcmkdMWppRyBRts .labelBox{stroke:#ccf;fill:#ECECFF}#mermaid-svg-QXcmkdMWppRyBRts .labelText,#mermaid-svg-QXcmkdMWppRyBRts .labelText>tspan{fill:#000;stroke:none}#mermaid-svg-QXcmkdMWppRyBRts .loopText,#mermaid-svg-QXcmkdMWppRyBRts .loopText>tspan{fill:#000;stroke:none}#mermaid-svg-QXcmkdMWppRyBRts .loopLine{stroke-width:2px;stroke-dasharray:2, 2;stroke:#ccf;fill:#ccf}#mermaid-svg-QXcmkdMWppRyBRts .note{stroke:#aa3;fill:#fff5ad}#mermaid-svg-QXcmkdMWppRyBRts .noteText,#mermaid-svg-QXcmkdMWppRyBRts .noteText>tspan{fill:#000;stroke:none}#mermaid-svg-QXcmkdMWppRyBRts .activation0{fill:#f4f4f4;stroke:#666}#mermaid-svg-QXcmkdMWppRyBRts .activation1{fill:#f4f4f4;stroke:#666}#mermaid-svg-QXcmkdMWppRyBRts .activation2{fill:#f4f4f4;stroke:#666}#mermaid-svg-QXcmkdMWppRyBRts .mermaid-main-font{font-family:"trebuchet ms", verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .section{stroke:none;opacity:0.2}#mermaid-svg-QXcmkdMWppRyBRts .section0{fill:rgba(102,102,255,0.49)}#mermaid-svg-QXcmkdMWppRyBRts .section2{fill:#fff400}#mermaid-svg-QXcmkdMWppRyBRts .section1,#mermaid-svg-QXcmkdMWppRyBRts .section3{fill:#fff;opacity:0.2}#mermaid-svg-QXcmkdMWppRyBRts .sectionTitle0{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .sectionTitle1{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .sectionTitle2{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .sectionTitle3{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .sectionTitle{text-anchor:start;font-size:11px;text-height:14px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .grid .tick{stroke:#d3d3d3;opacity:0.8;shape-rendering:crispEdges}#mermaid-svg-QXcmkdMWppRyBRts .grid .tick text{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .grid path{stroke-width:0}#mermaid-svg-QXcmkdMWppRyBRts .today{fill:none;stroke:red;stroke-width:2px}#mermaid-svg-QXcmkdMWppRyBRts .task{stroke-width:2}#mermaid-svg-QXcmkdMWppRyBRts .taskText{text-anchor:middle;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .taskText:not([font-size]){font-size:11px}#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutsideRight{fill:#000;text-anchor:start;font-size:11px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutsideLeft{fill:#000;text-anchor:end;font-size:11px}#mermaid-svg-QXcmkdMWppRyBRts .task.clickable{cursor:pointer}#mermaid-svg-QXcmkdMWppRyBRts .taskText.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutsideLeft.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutsideRight.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-QXcmkdMWppRyBRts .taskText0,#mermaid-svg-QXcmkdMWppRyBRts .taskText1,#mermaid-svg-QXcmkdMWppRyBRts .taskText2,#mermaid-svg-QXcmkdMWppRyBRts .taskText3{fill:#fff}#mermaid-svg-QXcmkdMWppRyBRts .task0,#mermaid-svg-QXcmkdMWppRyBRts .task1,#mermaid-svg-QXcmkdMWppRyBRts .task2,#mermaid-svg-QXcmkdMWppRyBRts .task3{fill:#8a90dd;stroke:#534fbc}#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutside0,#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutside2{fill:#000}#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutside1,#mermaid-svg-QXcmkdMWppRyBRts .taskTextOutside3{fill:#000}#mermaid-svg-QXcmkdMWppRyBRts .active0,#mermaid-svg-QXcmkdMWppRyBRts .active1,#mermaid-svg-QXcmkdMWppRyBRts .active2,#mermaid-svg-QXcmkdMWppRyBRts .active3{fill:#bfc7ff;stroke:#534fbc}#mermaid-svg-QXcmkdMWppRyBRts .activeText0,#mermaid-svg-QXcmkdMWppRyBRts .activeText1,#mermaid-svg-QXcmkdMWppRyBRts .activeText2,#mermaid-svg-QXcmkdMWppRyBRts .activeText3{fill:#000 !important}#mermaid-svg-QXcmkdMWppRyBRts .done0,#mermaid-svg-QXcmkdMWppRyBRts .done1,#mermaid-svg-QXcmkdMWppRyBRts .done2,#mermaid-svg-QXcmkdMWppRyBRts .done3{stroke:grey;fill:#d3d3d3;stroke-width:2}#mermaid-svg-QXcmkdMWppRyBRts .doneText0,#mermaid-svg-QXcmkdMWppRyBRts .doneText1,#mermaid-svg-QXcmkdMWppRyBRts .doneText2,#mermaid-svg-QXcmkdMWppRyBRts .doneText3{fill:#000 !important}#mermaid-svg-QXcmkdMWppRyBRts .crit0,#mermaid-svg-QXcmkdMWppRyBRts .crit1,#mermaid-svg-QXcmkdMWppRyBRts .crit2,#mermaid-svg-QXcmkdMWppRyBRts .crit3{stroke:#f88;fill:red;stroke-width:2}#mermaid-svg-QXcmkdMWppRyBRts .activeCrit0,#mermaid-svg-QXcmkdMWppRyBRts .activeCrit1,#mermaid-svg-QXcmkdMWppRyBRts .activeCrit2,#mermaid-svg-QXcmkdMWppRyBRts .activeCrit3{stroke:#f88;fill:#bfc7ff;stroke-width:2}#mermaid-svg-QXcmkdMWppRyBRts .doneCrit0,#mermaid-svg-QXcmkdMWppRyBRts .doneCrit1,#mermaid-svg-QXcmkdMWppRyBRts .doneCrit2,#mermaid-svg-QXcmkdMWppRyBRts .doneCrit3{stroke:#f88;fill:#d3d3d3;stroke-width:2;cursor:pointer;shape-rendering:crispEdges}#mermaid-svg-QXcmkdMWppRyBRts .milestone{transform:rotate(45deg) scale(0.8, 0.8)}#mermaid-svg-QXcmkdMWppRyBRts .milestoneText{font-style:italic}#mermaid-svg-QXcmkdMWppRyBRts .doneCritText0,#mermaid-svg-QXcmkdMWppRyBRts .doneCritText1,#mermaid-svg-QXcmkdMWppRyBRts .doneCritText2,#mermaid-svg-QXcmkdMWppRyBRts .doneCritText3{fill:#000 !important}#mermaid-svg-QXcmkdMWppRyBRts .activeCritText0,#mermaid-svg-QXcmkdMWppRyBRts .activeCritText1,#mermaid-svg-QXcmkdMWppRyBRts .activeCritText2,#mermaid-svg-QXcmkdMWppRyBRts .activeCritText3{fill:#000 !important}#mermaid-svg-QXcmkdMWppRyBRts .titleText{text-anchor:middle;font-size:18px;fill:#000;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts g.classGroup text{fill:#9370db;stroke:none;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);font-size:10px}#mermaid-svg-QXcmkdMWppRyBRts g.classGroup text .title{font-weight:bolder}#mermaid-svg-QXcmkdMWppRyBRts g.clickable{cursor:pointer}#mermaid-svg-QXcmkdMWppRyBRts g.classGroup rect{fill:#ECECFF;stroke:#9370db}#mermaid-svg-QXcmkdMWppRyBRts g.classGroup line{stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts .classLabel .box{stroke:none;stroke-width:0;fill:#ECECFF;opacity:0.5}#mermaid-svg-QXcmkdMWppRyBRts .classLabel .label{fill:#9370db;font-size:10px}#mermaid-svg-QXcmkdMWppRyBRts .relation{stroke:#9370db;stroke-width:1;fill:none}#mermaid-svg-QXcmkdMWppRyBRts .dashed-line{stroke-dasharray:3}#mermaid-svg-QXcmkdMWppRyBRts #compositionStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #compositionEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #aggregationStart{fill:#ECECFF;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #aggregationEnd{fill:#ECECFF;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #dependencyStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #dependencyEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #extensionStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts #extensionEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts .commit-id,#mermaid-svg-QXcmkdMWppRyBRts .commit-msg,#mermaid-svg-QXcmkdMWppRyBRts .branch-label{fill:lightgrey;color:lightgrey;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .pieTitleText{text-anchor:middle;font-size:25px;fill:#000;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .slice{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts g.stateGroup text{fill:#9370db;stroke:none;font-size:10px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts g.stateGroup text{fill:#9370db;fill:#333;stroke:none;font-size:10px}#mermaid-svg-QXcmkdMWppRyBRts g.statediagram-cluster .cluster-label text{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts g.stateGroup .state-title{font-weight:bolder;fill:#000}#mermaid-svg-QXcmkdMWppRyBRts g.stateGroup rect{fill:#ECECFF;stroke:#9370db}#mermaid-svg-QXcmkdMWppRyBRts g.stateGroup line{stroke:#9370db;stroke-width:1}#mermaid-svg-QXcmkdMWppRyBRts .transition{stroke:#9370db;stroke-width:1;fill:none}#mermaid-svg-QXcmkdMWppRyBRts .stateGroup .composit{fill:white;border-bottom:1px}#mermaid-svg-QXcmkdMWppRyBRts .stateGroup .alt-composit{fill:#e0e0e0;border-bottom:1px}#mermaid-svg-QXcmkdMWppRyBRts .state-note{stroke:#aa3;fill:#fff5ad}#mermaid-svg-QXcmkdMWppRyBRts .state-note text{fill:black;stroke:none;font-size:10px}#mermaid-svg-QXcmkdMWppRyBRts .stateLabel .box{stroke:none;stroke-width:0;fill:#ECECFF;opacity:0.7}#mermaid-svg-QXcmkdMWppRyBRts .edgeLabel text{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .stateLabel text{fill:#000;font-size:10px;font-weight:bold;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-QXcmkdMWppRyBRts .node circle.state-start{fill:black;stroke:black}#mermaid-svg-QXcmkdMWppRyBRts .node circle.state-end{fill:black;stroke:white;stroke-width:1.5}#mermaid-svg-QXcmkdMWppRyBRts #statediagram-barbEnd{fill:#9370db}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-cluster rect{fill:#ECECFF;stroke:#9370db;stroke-width:1px}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-cluster rect.outer{rx:5px;ry:5px}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-state .divider{stroke:#9370db}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-state .title-state{rx:5px;ry:5px}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-cluster.statediagram-cluster .inner{fill:white}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-cluster.statediagram-cluster-alt .inner{fill:#e0e0e0}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-cluster .inner{rx:0;ry:0}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-state rect.basic{rx:5px;ry:5px}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-state rect.divider{stroke-dasharray:10,10;fill:#efefef}#mermaid-svg-QXcmkdMWppRyBRts .note-edge{stroke-dasharray:5}#mermaid-svg-QXcmkdMWppRyBRts .statediagram-note rect{fill:#fff5ad;stroke:#aa3;stroke-width:1px;rx:0;ry:0}:root{--mermaid-font-family: '"trebuchet ms", verdana, arial';--mermaid-font-family: "Comic Sans MS", "Comic Sans", cursive}#mermaid-svg-QXcmkdMWppRyBRts .error-icon{fill:#522}#mermaid-svg-QXcmkdMWppRyBRts .error-text{fill:#522;stroke:#522}#mermaid-svg-QXcmkdMWppRyBRts .edge-thickness-normal{stroke-width:2px}#mermaid-svg-QXcmkdMWppRyBRts .edge-thickness-thick{stroke-width:3.5px}#mermaid-svg-QXcmkdMWppRyBRts .edge-pattern-solid{stroke-dasharray:0}#mermaid-svg-QXcmkdMWppRyBRts .edge-pattern-dashed{stroke-dasharray:3}#mermaid-svg-QXcmkdMWppRyBRts .edge-pattern-dotted{stroke-dasharray:2}#mermaid-svg-QXcmkdMWppRyBRts .marker{fill:#333}#mermaid-svg-QXcmkdMWppRyBRts .marker.cross{stroke:#333}:root { --mermaid-font-family: "trebuchet ms", verdana, arial;}#mermaid-svg-QXcmkdMWppRyBRts {color: rgba(0, 0, 0, 0.75);font: ;}

对样本进行测序
使用fastp对下机数据进行修剪
去除宿主的基因
生成样本的profile

1、对样本进行测序的步骤交给自动化测序仪器
2、trim ing然后将测序仪得到的下机数据进行修剪(这一步主要是为了去除低质量的碱基,通过对低质量碱基进行修剪,可以保留更多的reads,reads更多有利于profile,而但却会损失base(>10%),不利于组装)切两边应该是测序不知道哪边是尾巴。

观点:
因为测序的时候荧光标记的ddNTP是过量的,反应结束后一并沉淀下来干扰了小分zhi子量的引物后面几个核苷酸的信号。
至于后面一般800bp以上的区段也是不可信的,这也是目前测序公司技术优劣的一个指标,技术好的能读到1100bp,技术差的可能只有500bp。
原因就在于用的分离胶好不好,因为分子量越大,同样相差1个核苷酸所造成的影响越小。21个bp的片段和22个bp的片段可以分的比较开,而801bp和802bp就几乎贴在一起了,越到后面贴的越紧,最后就连成一片分不清谁是谁了。
具体从哪里到哪里可信要看测序结果的信号峰。强烈的单峰才是有效的。

3、将过滤后的数据去除宿主基因,在取样时是存在宿主基因的是存在所以有必要去除,同时要屏蔽掉一些与宿主高度同源的基因
原因:1、后期profiling分析不需要用到宿主DNA这样可以减少数据以及后期的处理。
2、做组装,会影响组装的准确度B。
3、减少之后步骤不想关数据量从而减少计算量。

4、生成profile-序列数据谱
样本数据质量统计情况表格来源MetaPhlAn3软件处理得到

project mining
Library ID 样本名
Raw_reads_count 总的reads数量
Raw_base_out 总的组装数
Read_length 读取序列的长度
Raw_Q30_ratio 原始质量
Dupliaction 重复数
clean_reads_count triming后的数据也就是剪去两端后
clean_reads_ratio 占总的比例
clean_Q30_ratio 质量
clean_bases_count 总组装数
clean_bases_ratio 占比
Nohomo_reads_count 去除宿主后的reads数量
Nohomo_reads_ratio 去除宿主后reads占总的比重
Nohomo_bases_ratio 没有人类的组装后占总比重
----------- -------------------------------------

metaphlan3.profile.merge

ID进行比后物质分类
test1.mps_unkbown,这一列表示这个物质的样本中每种物种丰度
UNKNOWN为未入库的比例
K_Bactria为细菌总占比
K_Bactria_或|者其他都是在细菌中的小分类,他们之和等于细菌占比。

-----------------------------------------------------------------------------------------------------------------------------------------我他的心态有点炸,这个程序不是一般的用了snakemake以及一些程序原理不是特别了解所以我并不能完美的理解在这里我决定通过3天来实现部分操作
-----------------------------------------------------------------------------------------------------------------------------------------

实际操作

构建snakemake流程进行生物信息分析

8.25-26没有进展由于环境问题还是怎么的发现无法运行了

9.11恢复写实


==我之前步骤主要是获取种群分布情况,该方法通过对照组比较主要用于寻找可能标志种群或差异种群有利于进一步研究快速找到有关菌种 ==
COG分析就是蛋白质功能分析

#mermaid-svg-GfjTDdNj1pb8n507 .label{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);fill:#333;color:#333}#mermaid-svg-GfjTDdNj1pb8n507 .label text{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .node rect,#mermaid-svg-GfjTDdNj1pb8n507 .node circle,#mermaid-svg-GfjTDdNj1pb8n507 .node ellipse,#mermaid-svg-GfjTDdNj1pb8n507 .node polygon,#mermaid-svg-GfjTDdNj1pb8n507 .node path{fill:#ECECFF;stroke:#9370db;stroke-width:1px}#mermaid-svg-GfjTDdNj1pb8n507 .node .label{text-align:center;fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .node.clickable{cursor:pointer}#mermaid-svg-GfjTDdNj1pb8n507 .arrowheadPath{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .edgePath .path{stroke:#333;stroke-width:1.5px}#mermaid-svg-GfjTDdNj1pb8n507 .flowchart-link{stroke:#333;fill:none}#mermaid-svg-GfjTDdNj1pb8n507 .edgeLabel{background-color:#e8e8e8;text-align:center}#mermaid-svg-GfjTDdNj1pb8n507 .edgeLabel rect{opacity:0.9}#mermaid-svg-GfjTDdNj1pb8n507 .edgeLabel span{color:#333}#mermaid-svg-GfjTDdNj1pb8n507 .cluster rect{fill:#ffffde;stroke:#aa3;stroke-width:1px}#mermaid-svg-GfjTDdNj1pb8n507 .cluster text{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 div.mermaidTooltip{position:absolute;text-align:center;max-width:200px;padding:2px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);font-size:12px;background:#ffffde;border:1px solid #aa3;border-radius:2px;pointer-events:none;z-index:100}#mermaid-svg-GfjTDdNj1pb8n507 .actor{stroke:#ccf;fill:#ECECFF}#mermaid-svg-GfjTDdNj1pb8n507 text.actor>tspan{fill:#000;stroke:none}#mermaid-svg-GfjTDdNj1pb8n507 .actor-line{stroke:grey}#mermaid-svg-GfjTDdNj1pb8n507 .messageLine0{stroke-width:1.5;stroke-dasharray:none;stroke:#333}#mermaid-svg-GfjTDdNj1pb8n507 .messageLine1{stroke-width:1.5;stroke-dasharray:2, 2;stroke:#333}#mermaid-svg-GfjTDdNj1pb8n507 #arrowhead path{fill:#333;stroke:#333}#mermaid-svg-GfjTDdNj1pb8n507 .sequenceNumber{fill:#fff}#mermaid-svg-GfjTDdNj1pb8n507 #sequencenumber{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 #crosshead path{fill:#333;stroke:#333}#mermaid-svg-GfjTDdNj1pb8n507 .messageText{fill:#333;stroke:#333}#mermaid-svg-GfjTDdNj1pb8n507 .labelBox{stroke:#ccf;fill:#ECECFF}#mermaid-svg-GfjTDdNj1pb8n507 .labelText,#mermaid-svg-GfjTDdNj1pb8n507 .labelText>tspan{fill:#000;stroke:none}#mermaid-svg-GfjTDdNj1pb8n507 .loopText,#mermaid-svg-GfjTDdNj1pb8n507 .loopText>tspan{fill:#000;stroke:none}#mermaid-svg-GfjTDdNj1pb8n507 .loopLine{stroke-width:2px;stroke-dasharray:2, 2;stroke:#ccf;fill:#ccf}#mermaid-svg-GfjTDdNj1pb8n507 .note{stroke:#aa3;fill:#fff5ad}#mermaid-svg-GfjTDdNj1pb8n507 .noteText,#mermaid-svg-GfjTDdNj1pb8n507 .noteText>tspan{fill:#000;stroke:none}#mermaid-svg-GfjTDdNj1pb8n507 .activation0{fill:#f4f4f4;stroke:#666}#mermaid-svg-GfjTDdNj1pb8n507 .activation1{fill:#f4f4f4;stroke:#666}#mermaid-svg-GfjTDdNj1pb8n507 .activation2{fill:#f4f4f4;stroke:#666}#mermaid-svg-GfjTDdNj1pb8n507 .mermaid-main-font{font-family:"trebuchet ms", verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .section{stroke:none;opacity:0.2}#mermaid-svg-GfjTDdNj1pb8n507 .section0{fill:rgba(102,102,255,0.49)}#mermaid-svg-GfjTDdNj1pb8n507 .section2{fill:#fff400}#mermaid-svg-GfjTDdNj1pb8n507 .section1,#mermaid-svg-GfjTDdNj1pb8n507 .section3{fill:#fff;opacity:0.2}#mermaid-svg-GfjTDdNj1pb8n507 .sectionTitle0{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .sectionTitle1{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .sectionTitle2{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .sectionTitle3{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .sectionTitle{text-anchor:start;font-size:11px;text-height:14px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .grid .tick{stroke:#d3d3d3;opacity:0.8;shape-rendering:crispEdges}#mermaid-svg-GfjTDdNj1pb8n507 .grid .tick text{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .grid path{stroke-width:0}#mermaid-svg-GfjTDdNj1pb8n507 .today{fill:none;stroke:red;stroke-width:2px}#mermaid-svg-GfjTDdNj1pb8n507 .task{stroke-width:2}#mermaid-svg-GfjTDdNj1pb8n507 .taskText{text-anchor:middle;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .taskText:not([font-size]){font-size:11px}#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutsideRight{fill:#000;text-anchor:start;font-size:11px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutsideLeft{fill:#000;text-anchor:end;font-size:11px}#mermaid-svg-GfjTDdNj1pb8n507 .task.clickable{cursor:pointer}#mermaid-svg-GfjTDdNj1pb8n507 .taskText.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutsideLeft.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutsideRight.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-GfjTDdNj1pb8n507 .taskText0,#mermaid-svg-GfjTDdNj1pb8n507 .taskText1,#mermaid-svg-GfjTDdNj1pb8n507 .taskText2,#mermaid-svg-GfjTDdNj1pb8n507 .taskText3{fill:#fff}#mermaid-svg-GfjTDdNj1pb8n507 .task0,#mermaid-svg-GfjTDdNj1pb8n507 .task1,#mermaid-svg-GfjTDdNj1pb8n507 .task2,#mermaid-svg-GfjTDdNj1pb8n507 .task3{fill:#8a90dd;stroke:#534fbc}#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutside0,#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutside2{fill:#000}#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutside1,#mermaid-svg-GfjTDdNj1pb8n507 .taskTextOutside3{fill:#000}#mermaid-svg-GfjTDdNj1pb8n507 .active0,#mermaid-svg-GfjTDdNj1pb8n507 .active1,#mermaid-svg-GfjTDdNj1pb8n507 .active2,#mermaid-svg-GfjTDdNj1pb8n507 .active3{fill:#bfc7ff;stroke:#534fbc}#mermaid-svg-GfjTDdNj1pb8n507 .activeText0,#mermaid-svg-GfjTDdNj1pb8n507 .activeText1,#mermaid-svg-GfjTDdNj1pb8n507 .activeText2,#mermaid-svg-GfjTDdNj1pb8n507 .activeText3{fill:#000 !important}#mermaid-svg-GfjTDdNj1pb8n507 .done0,#mermaid-svg-GfjTDdNj1pb8n507 .done1,#mermaid-svg-GfjTDdNj1pb8n507 .done2,#mermaid-svg-GfjTDdNj1pb8n507 .done3{stroke:grey;fill:#d3d3d3;stroke-width:2}#mermaid-svg-GfjTDdNj1pb8n507 .doneText0,#mermaid-svg-GfjTDdNj1pb8n507 .doneText1,#mermaid-svg-GfjTDdNj1pb8n507 .doneText2,#mermaid-svg-GfjTDdNj1pb8n507 .doneText3{fill:#000 !important}#mermaid-svg-GfjTDdNj1pb8n507 .crit0,#mermaid-svg-GfjTDdNj1pb8n507 .crit1,#mermaid-svg-GfjTDdNj1pb8n507 .crit2,#mermaid-svg-GfjTDdNj1pb8n507 .crit3{stroke:#f88;fill:red;stroke-width:2}#mermaid-svg-GfjTDdNj1pb8n507 .activeCrit0,#mermaid-svg-GfjTDdNj1pb8n507 .activeCrit1,#mermaid-svg-GfjTDdNj1pb8n507 .activeCrit2,#mermaid-svg-GfjTDdNj1pb8n507 .activeCrit3{stroke:#f88;fill:#bfc7ff;stroke-width:2}#mermaid-svg-GfjTDdNj1pb8n507 .doneCrit0,#mermaid-svg-GfjTDdNj1pb8n507 .doneCrit1,#mermaid-svg-GfjTDdNj1pb8n507 .doneCrit2,#mermaid-svg-GfjTDdNj1pb8n507 .doneCrit3{stroke:#f88;fill:#d3d3d3;stroke-width:2;cursor:pointer;shape-rendering:crispEdges}#mermaid-svg-GfjTDdNj1pb8n507 .milestone{transform:rotate(45deg) scale(0.8, 0.8)}#mermaid-svg-GfjTDdNj1pb8n507 .milestoneText{font-style:italic}#mermaid-svg-GfjTDdNj1pb8n507 .doneCritText0,#mermaid-svg-GfjTDdNj1pb8n507 .doneCritText1,#mermaid-svg-GfjTDdNj1pb8n507 .doneCritText2,#mermaid-svg-GfjTDdNj1pb8n507 .doneCritText3{fill:#000 !important}#mermaid-svg-GfjTDdNj1pb8n507 .activeCritText0,#mermaid-svg-GfjTDdNj1pb8n507 .activeCritText1,#mermaid-svg-GfjTDdNj1pb8n507 .activeCritText2,#mermaid-svg-GfjTDdNj1pb8n507 .activeCritText3{fill:#000 !important}#mermaid-svg-GfjTDdNj1pb8n507 .titleText{text-anchor:middle;font-size:18px;fill:#000;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 g.classGroup text{fill:#9370db;stroke:none;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);font-size:10px}#mermaid-svg-GfjTDdNj1pb8n507 g.classGroup text .title{font-weight:bolder}#mermaid-svg-GfjTDdNj1pb8n507 g.clickable{cursor:pointer}#mermaid-svg-GfjTDdNj1pb8n507 g.classGroup rect{fill:#ECECFF;stroke:#9370db}#mermaid-svg-GfjTDdNj1pb8n507 g.classGroup line{stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 .classLabel .box{stroke:none;stroke-width:0;fill:#ECECFF;opacity:0.5}#mermaid-svg-GfjTDdNj1pb8n507 .classLabel .label{fill:#9370db;font-size:10px}#mermaid-svg-GfjTDdNj1pb8n507 .relation{stroke:#9370db;stroke-width:1;fill:none}#mermaid-svg-GfjTDdNj1pb8n507 .dashed-line{stroke-dasharray:3}#mermaid-svg-GfjTDdNj1pb8n507 #compositionStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #compositionEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #aggregationStart{fill:#ECECFF;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #aggregationEnd{fill:#ECECFF;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #dependencyStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #dependencyEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #extensionStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 #extensionEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 .commit-id,#mermaid-svg-GfjTDdNj1pb8n507 .commit-msg,#mermaid-svg-GfjTDdNj1pb8n507 .branch-label{fill:lightgrey;color:lightgrey;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .pieTitleText{text-anchor:middle;font-size:25px;fill:#000;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .slice{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 g.stateGroup text{fill:#9370db;stroke:none;font-size:10px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 g.stateGroup text{fill:#9370db;fill:#333;stroke:none;font-size:10px}#mermaid-svg-GfjTDdNj1pb8n507 g.statediagram-cluster .cluster-label text{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 g.stateGroup .state-title{font-weight:bolder;fill:#000}#mermaid-svg-GfjTDdNj1pb8n507 g.stateGroup rect{fill:#ECECFF;stroke:#9370db}#mermaid-svg-GfjTDdNj1pb8n507 g.stateGroup line{stroke:#9370db;stroke-width:1}#mermaid-svg-GfjTDdNj1pb8n507 .transition{stroke:#9370db;stroke-width:1;fill:none}#mermaid-svg-GfjTDdNj1pb8n507 .stateGroup .composit{fill:white;border-bottom:1px}#mermaid-svg-GfjTDdNj1pb8n507 .stateGroup .alt-composit{fill:#e0e0e0;border-bottom:1px}#mermaid-svg-GfjTDdNj1pb8n507 .state-note{stroke:#aa3;fill:#fff5ad}#mermaid-svg-GfjTDdNj1pb8n507 .state-note text{fill:black;stroke:none;font-size:10px}#mermaid-svg-GfjTDdNj1pb8n507 .stateLabel .box{stroke:none;stroke-width:0;fill:#ECECFF;opacity:0.7}#mermaid-svg-GfjTDdNj1pb8n507 .edgeLabel text{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .stateLabel text{fill:#000;font-size:10px;font-weight:bold;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-GfjTDdNj1pb8n507 .node circle.state-start{fill:black;stroke:black}#mermaid-svg-GfjTDdNj1pb8n507 .node circle.state-end{fill:black;stroke:white;stroke-width:1.5}#mermaid-svg-GfjTDdNj1pb8n507 #statediagram-barbEnd{fill:#9370db}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-cluster rect{fill:#ECECFF;stroke:#9370db;stroke-width:1px}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-cluster rect.outer{rx:5px;ry:5px}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-state .divider{stroke:#9370db}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-state .title-state{rx:5px;ry:5px}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-cluster.statediagram-cluster .inner{fill:white}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-cluster.statediagram-cluster-alt .inner{fill:#e0e0e0}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-cluster .inner{rx:0;ry:0}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-state rect.basic{rx:5px;ry:5px}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-state rect.divider{stroke-dasharray:10,10;fill:#efefef}#mermaid-svg-GfjTDdNj1pb8n507 .note-edge{stroke-dasharray:5}#mermaid-svg-GfjTDdNj1pb8n507 .statediagram-note rect{fill:#fff5ad;stroke:#aa3;stroke-width:1px;rx:0;ry:0}:root{--mermaid-font-family: '"trebuchet ms", verdana, arial';--mermaid-font-family: "Comic Sans MS", "Comic Sans", cursive}#mermaid-svg-GfjTDdNj1pb8n507 .error-icon{fill:#522}#mermaid-svg-GfjTDdNj1pb8n507 .error-text{fill:#522;stroke:#522}#mermaid-svg-GfjTDdNj1pb8n507 .edge-thickness-normal{stroke-width:2px}#mermaid-svg-GfjTDdNj1pb8n507 .edge-thickness-thick{stroke-width:3.5px}#mermaid-svg-GfjTDdNj1pb8n507 .edge-pattern-solid{stroke-dasharray:0}#mermaid-svg-GfjTDdNj1pb8n507 .edge-pattern-dashed{stroke-dasharray:3}#mermaid-svg-GfjTDdNj1pb8n507 .edge-pattern-dotted{stroke-dasharray:2}#mermaid-svg-GfjTDdNj1pb8n507 .marker{fill:#333}#mermaid-svg-GfjTDdNj1pb8n507 .marker.cross{stroke:#333}:root { --mermaid-font-family: "trebuchet ms", verdana, arial;}#mermaid-svg-GfjTDdNj1pb8n507 {color: rgba(0, 0, 0, 0.75);font: ;}样本信息snakemake将要进行分析的样本加准备好进入snakemake样本信息snakemake
#mermaid-svg-CtmJLSOjPch8RHVP .label{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);fill:#333;color:#333}#mermaid-svg-CtmJLSOjPch8RHVP .label text{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .node rect,#mermaid-svg-CtmJLSOjPch8RHVP .node circle,#mermaid-svg-CtmJLSOjPch8RHVP .node ellipse,#mermaid-svg-CtmJLSOjPch8RHVP .node polygon,#mermaid-svg-CtmJLSOjPch8RHVP .node path{fill:#ECECFF;stroke:#9370db;stroke-width:1px}#mermaid-svg-CtmJLSOjPch8RHVP .node .label{text-align:center;fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .node.clickable{cursor:pointer}#mermaid-svg-CtmJLSOjPch8RHVP .arrowheadPath{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .edgePath .path{stroke:#333;stroke-width:1.5px}#mermaid-svg-CtmJLSOjPch8RHVP .flowchart-link{stroke:#333;fill:none}#mermaid-svg-CtmJLSOjPch8RHVP .edgeLabel{background-color:#e8e8e8;text-align:center}#mermaid-svg-CtmJLSOjPch8RHVP .edgeLabel rect{opacity:0.9}#mermaid-svg-CtmJLSOjPch8RHVP .edgeLabel span{color:#333}#mermaid-svg-CtmJLSOjPch8RHVP .cluster rect{fill:#ffffde;stroke:#aa3;stroke-width:1px}#mermaid-svg-CtmJLSOjPch8RHVP .cluster text{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP div.mermaidTooltip{position:absolute;text-align:center;max-width:200px;padding:2px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);font-size:12px;background:#ffffde;border:1px solid #aa3;border-radius:2px;pointer-events:none;z-index:100}#mermaid-svg-CtmJLSOjPch8RHVP .actor{stroke:#ccf;fill:#ECECFF}#mermaid-svg-CtmJLSOjPch8RHVP text.actor>tspan{fill:#000;stroke:none}#mermaid-svg-CtmJLSOjPch8RHVP .actor-line{stroke:grey}#mermaid-svg-CtmJLSOjPch8RHVP .messageLine0{stroke-width:1.5;stroke-dasharray:none;stroke:#333}#mermaid-svg-CtmJLSOjPch8RHVP .messageLine1{stroke-width:1.5;stroke-dasharray:2, 2;stroke:#333}#mermaid-svg-CtmJLSOjPch8RHVP #arrowhead path{fill:#333;stroke:#333}#mermaid-svg-CtmJLSOjPch8RHVP .sequenceNumber{fill:#fff}#mermaid-svg-CtmJLSOjPch8RHVP #sequencenumber{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP #crosshead path{fill:#333;stroke:#333}#mermaid-svg-CtmJLSOjPch8RHVP .messageText{fill:#333;stroke:#333}#mermaid-svg-CtmJLSOjPch8RHVP .labelBox{stroke:#ccf;fill:#ECECFF}#mermaid-svg-CtmJLSOjPch8RHVP .labelText,#mermaid-svg-CtmJLSOjPch8RHVP .labelText>tspan{fill:#000;stroke:none}#mermaid-svg-CtmJLSOjPch8RHVP .loopText,#mermaid-svg-CtmJLSOjPch8RHVP .loopText>tspan{fill:#000;stroke:none}#mermaid-svg-CtmJLSOjPch8RHVP .loopLine{stroke-width:2px;stroke-dasharray:2, 2;stroke:#ccf;fill:#ccf}#mermaid-svg-CtmJLSOjPch8RHVP .note{stroke:#aa3;fill:#fff5ad}#mermaid-svg-CtmJLSOjPch8RHVP .noteText,#mermaid-svg-CtmJLSOjPch8RHVP .noteText>tspan{fill:#000;stroke:none}#mermaid-svg-CtmJLSOjPch8RHVP .activation0{fill:#f4f4f4;stroke:#666}#mermaid-svg-CtmJLSOjPch8RHVP .activation1{fill:#f4f4f4;stroke:#666}#mermaid-svg-CtmJLSOjPch8RHVP .activation2{fill:#f4f4f4;stroke:#666}#mermaid-svg-CtmJLSOjPch8RHVP .mermaid-main-font{font-family:"trebuchet ms", verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .section{stroke:none;opacity:0.2}#mermaid-svg-CtmJLSOjPch8RHVP .section0{fill:rgba(102,102,255,0.49)}#mermaid-svg-CtmJLSOjPch8RHVP .section2{fill:#fff400}#mermaid-svg-CtmJLSOjPch8RHVP .section1,#mermaid-svg-CtmJLSOjPch8RHVP .section3{fill:#fff;opacity:0.2}#mermaid-svg-CtmJLSOjPch8RHVP .sectionTitle0{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .sectionTitle1{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .sectionTitle2{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .sectionTitle3{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .sectionTitle{text-anchor:start;font-size:11px;text-height:14px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .grid .tick{stroke:#d3d3d3;opacity:0.8;shape-rendering:crispEdges}#mermaid-svg-CtmJLSOjPch8RHVP .grid .tick text{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .grid path{stroke-width:0}#mermaid-svg-CtmJLSOjPch8RHVP .today{fill:none;stroke:red;stroke-width:2px}#mermaid-svg-CtmJLSOjPch8RHVP .task{stroke-width:2}#mermaid-svg-CtmJLSOjPch8RHVP .taskText{text-anchor:middle;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .taskText:not([font-size]){font-size:11px}#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutsideRight{fill:#000;text-anchor:start;font-size:11px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutsideLeft{fill:#000;text-anchor:end;font-size:11px}#mermaid-svg-CtmJLSOjPch8RHVP .task.clickable{cursor:pointer}#mermaid-svg-CtmJLSOjPch8RHVP .taskText.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutsideLeft.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutsideRight.clickable{cursor:pointer;fill:#003163 !important;font-weight:bold}#mermaid-svg-CtmJLSOjPch8RHVP .taskText0,#mermaid-svg-CtmJLSOjPch8RHVP .taskText1,#mermaid-svg-CtmJLSOjPch8RHVP .taskText2,#mermaid-svg-CtmJLSOjPch8RHVP .taskText3{fill:#fff}#mermaid-svg-CtmJLSOjPch8RHVP .task0,#mermaid-svg-CtmJLSOjPch8RHVP .task1,#mermaid-svg-CtmJLSOjPch8RHVP .task2,#mermaid-svg-CtmJLSOjPch8RHVP .task3{fill:#8a90dd;stroke:#534fbc}#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutside0,#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutside2{fill:#000}#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutside1,#mermaid-svg-CtmJLSOjPch8RHVP .taskTextOutside3{fill:#000}#mermaid-svg-CtmJLSOjPch8RHVP .active0,#mermaid-svg-CtmJLSOjPch8RHVP .active1,#mermaid-svg-CtmJLSOjPch8RHVP .active2,#mermaid-svg-CtmJLSOjPch8RHVP .active3{fill:#bfc7ff;stroke:#534fbc}#mermaid-svg-CtmJLSOjPch8RHVP .activeText0,#mermaid-svg-CtmJLSOjPch8RHVP .activeText1,#mermaid-svg-CtmJLSOjPch8RHVP .activeText2,#mermaid-svg-CtmJLSOjPch8RHVP .activeText3{fill:#000 !important}#mermaid-svg-CtmJLSOjPch8RHVP .done0,#mermaid-svg-CtmJLSOjPch8RHVP .done1,#mermaid-svg-CtmJLSOjPch8RHVP .done2,#mermaid-svg-CtmJLSOjPch8RHVP .done3{stroke:grey;fill:#d3d3d3;stroke-width:2}#mermaid-svg-CtmJLSOjPch8RHVP .doneText0,#mermaid-svg-CtmJLSOjPch8RHVP .doneText1,#mermaid-svg-CtmJLSOjPch8RHVP .doneText2,#mermaid-svg-CtmJLSOjPch8RHVP .doneText3{fill:#000 !important}#mermaid-svg-CtmJLSOjPch8RHVP .crit0,#mermaid-svg-CtmJLSOjPch8RHVP .crit1,#mermaid-svg-CtmJLSOjPch8RHVP .crit2,#mermaid-svg-CtmJLSOjPch8RHVP .crit3{stroke:#f88;fill:red;stroke-width:2}#mermaid-svg-CtmJLSOjPch8RHVP .activeCrit0,#mermaid-svg-CtmJLSOjPch8RHVP .activeCrit1,#mermaid-svg-CtmJLSOjPch8RHVP .activeCrit2,#mermaid-svg-CtmJLSOjPch8RHVP .activeCrit3{stroke:#f88;fill:#bfc7ff;stroke-width:2}#mermaid-svg-CtmJLSOjPch8RHVP .doneCrit0,#mermaid-svg-CtmJLSOjPch8RHVP .doneCrit1,#mermaid-svg-CtmJLSOjPch8RHVP .doneCrit2,#mermaid-svg-CtmJLSOjPch8RHVP .doneCrit3{stroke:#f88;fill:#d3d3d3;stroke-width:2;cursor:pointer;shape-rendering:crispEdges}#mermaid-svg-CtmJLSOjPch8RHVP .milestone{transform:rotate(45deg) scale(0.8, 0.8)}#mermaid-svg-CtmJLSOjPch8RHVP .milestoneText{font-style:italic}#mermaid-svg-CtmJLSOjPch8RHVP .doneCritText0,#mermaid-svg-CtmJLSOjPch8RHVP .doneCritText1,#mermaid-svg-CtmJLSOjPch8RHVP .doneCritText2,#mermaid-svg-CtmJLSOjPch8RHVP .doneCritText3{fill:#000 !important}#mermaid-svg-CtmJLSOjPch8RHVP .activeCritText0,#mermaid-svg-CtmJLSOjPch8RHVP .activeCritText1,#mermaid-svg-CtmJLSOjPch8RHVP .activeCritText2,#mermaid-svg-CtmJLSOjPch8RHVP .activeCritText3{fill:#000 !important}#mermaid-svg-CtmJLSOjPch8RHVP .titleText{text-anchor:middle;font-size:18px;fill:#000;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP g.classGroup text{fill:#9370db;stroke:none;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family);font-size:10px}#mermaid-svg-CtmJLSOjPch8RHVP g.classGroup text .title{font-weight:bolder}#mermaid-svg-CtmJLSOjPch8RHVP g.clickable{cursor:pointer}#mermaid-svg-CtmJLSOjPch8RHVP g.classGroup rect{fill:#ECECFF;stroke:#9370db}#mermaid-svg-CtmJLSOjPch8RHVP g.classGroup line{stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP .classLabel .box{stroke:none;stroke-width:0;fill:#ECECFF;opacity:0.5}#mermaid-svg-CtmJLSOjPch8RHVP .classLabel .label{fill:#9370db;font-size:10px}#mermaid-svg-CtmJLSOjPch8RHVP .relation{stroke:#9370db;stroke-width:1;fill:none}#mermaid-svg-CtmJLSOjPch8RHVP .dashed-line{stroke-dasharray:3}#mermaid-svg-CtmJLSOjPch8RHVP #compositionStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #compositionEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #aggregationStart{fill:#ECECFF;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #aggregationEnd{fill:#ECECFF;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #dependencyStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #dependencyEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #extensionStart{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP #extensionEnd{fill:#9370db;stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP .commit-id,#mermaid-svg-CtmJLSOjPch8RHVP .commit-msg,#mermaid-svg-CtmJLSOjPch8RHVP .branch-label{fill:lightgrey;color:lightgrey;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .pieTitleText{text-anchor:middle;font-size:25px;fill:#000;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .slice{font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP g.stateGroup text{fill:#9370db;stroke:none;font-size:10px;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP g.stateGroup text{fill:#9370db;fill:#333;stroke:none;font-size:10px}#mermaid-svg-CtmJLSOjPch8RHVP g.statediagram-cluster .cluster-label text{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP g.stateGroup .state-title{font-weight:bolder;fill:#000}#mermaid-svg-CtmJLSOjPch8RHVP g.stateGroup rect{fill:#ECECFF;stroke:#9370db}#mermaid-svg-CtmJLSOjPch8RHVP g.stateGroup line{stroke:#9370db;stroke-width:1}#mermaid-svg-CtmJLSOjPch8RHVP .transition{stroke:#9370db;stroke-width:1;fill:none}#mermaid-svg-CtmJLSOjPch8RHVP .stateGroup .composit{fill:white;border-bottom:1px}#mermaid-svg-CtmJLSOjPch8RHVP .stateGroup .alt-composit{fill:#e0e0e0;border-bottom:1px}#mermaid-svg-CtmJLSOjPch8RHVP .state-note{stroke:#aa3;fill:#fff5ad}#mermaid-svg-CtmJLSOjPch8RHVP .state-note text{fill:black;stroke:none;font-size:10px}#mermaid-svg-CtmJLSOjPch8RHVP .stateLabel .box{stroke:none;stroke-width:0;fill:#ECECFF;opacity:0.7}#mermaid-svg-CtmJLSOjPch8RHVP .edgeLabel text{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .stateLabel text{fill:#000;font-size:10px;font-weight:bold;font-family:'trebuchet ms', verdana, arial;font-family:var(--mermaid-font-family)}#mermaid-svg-CtmJLSOjPch8RHVP .node circle.state-start{fill:black;stroke:black}#mermaid-svg-CtmJLSOjPch8RHVP .node circle.state-end{fill:black;stroke:white;stroke-width:1.5}#mermaid-svg-CtmJLSOjPch8RHVP #statediagram-barbEnd{fill:#9370db}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-cluster rect{fill:#ECECFF;stroke:#9370db;stroke-width:1px}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-cluster rect.outer{rx:5px;ry:5px}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-state .divider{stroke:#9370db}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-state .title-state{rx:5px;ry:5px}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-cluster.statediagram-cluster .inner{fill:white}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-cluster.statediagram-cluster-alt .inner{fill:#e0e0e0}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-cluster .inner{rx:0;ry:0}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-state rect.basic{rx:5px;ry:5px}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-state rect.divider{stroke-dasharray:10,10;fill:#efefef}#mermaid-svg-CtmJLSOjPch8RHVP .note-edge{stroke-dasharray:5}#mermaid-svg-CtmJLSOjPch8RHVP .statediagram-note rect{fill:#fff5ad;stroke:#aa3;stroke-width:1px;rx:0;ry:0}:root{--mermaid-font-family: '"trebuchet ms", verdana, arial';--mermaid-font-family: "Comic Sans MS", "Comic Sans", cursive}#mermaid-svg-CtmJLSOjPch8RHVP .error-icon{fill:#522}#mermaid-svg-CtmJLSOjPch8RHVP .error-text{fill:#522;stroke:#522}#mermaid-svg-CtmJLSOjPch8RHVP .edge-thickness-normal{stroke-width:2px}#mermaid-svg-CtmJLSOjPch8RHVP .edge-thickness-thick{stroke-width:3.5px}#mermaid-svg-CtmJLSOjPch8RHVP .edge-pattern-solid{stroke-dasharray:0}#mermaid-svg-CtmJLSOjPch8RHVP .edge-pattern-dashed{stroke-dasharray:3}#mermaid-svg-CtmJLSOjPch8RHVP .edge-pattern-dotted{stroke-dasharray:2}#mermaid-svg-CtmJLSOjPch8RHVP .marker{fill:#333}#mermaid-svg-CtmJLSOjPch8RHVP .marker.cross{stroke:#333}:root { --mermaid-font-family: "trebuchet ms", verdana, arial;}#mermaid-svg-CtmJLSOjPch8RHVP {color: rgba(0, 0, 0, 0.75);font: ;}

样本一
fastp
bowtie2
metaphlan3
样本二
fastp
bowtie2
metaphlan3
样本三
fastp
bowtie2
metaphlan3
合并
prifiles

构建流程----非snakemake 的并运行,模仿snakemake的运行模式

细菌宏基因组分析中菌群中种群分布情况的获取

# 进行第一步质滤
mkdir -p 1.assay/fastp/ 1.assay/fastp/json/ 1.assay/fastp/html/ 1.assay/logs/fastp/ 2.runid/
cat sample.txt|while read line;doecho $line>lineinput_1=$(awk '{print$2}' line) input_2=$(awk '{print$3}' line) id=$(awk '{print$1}' line) export id=$idexport input_1=$input_1export input_2=$input_2echo /ldfssz1/ST_META/share/User/zhujie/.conda/envs/bioenv/bin/fastp -i $input_1 -I $input_2 -o 1.assay/fastp/$id.trimmed.1.fq.gz -O 1.assay/fastp/$id.trimmed.2.fq.gz -w 8 --length_required 30 --disable_adapter_trimming -j 1.assay/fastp/json/$id.json -h 1.assay/fastp/html/$id.html 2\>1.assay/logs/fastp/$id.fastp.log >fastp.shqsub -clear -cwd -l vf=8g,num_proc=4 -P P18Z10200N0127 -binding linear:4  -q st.q -o o.txt -e e.txt fastp.sh # sh  fastp.sh $id $input_1 $input_2#qsub -clear -cwd -l vf=4g,num_proc=4 -P P18Z10200N0127 -binding linear:4  -q st.q dofastp.sh#bowtie2去除人类的基因,进行比对,将人类的基因去除,由于我们要分析的数据为细菌基因,而人类基因常常占大多数所以去除改部分基因可以提高细菌中物种显著水平
bowtie2_index="/ldfssz1/ST_META/P18Z10200N0127_VC/tianliu/todo/03.stlfr/metastlfr/database/hg38/hg38"
#由于任务集群上运行的而该步骤要等前一步骤的样本完成才能进行所以要增加一个守门关#中间增加一个控制程序检测上一个程序运行了没有运行了就可以进行下一步操作
make -p 1.assay/logs/bowtie2/ 2.runid/ 1.assay/bowtie2/
cut -f 1 sample.txt|while read id;
doecho   /ldfssz1/ST_META/share/User/zhujie/.conda/envs/bioenv/bin/bowtie2 --very-sensitive -p 8 -x $bowtie2_index -1 1.assay/fastp/$id.trimmed.1.fq.gz -2 1.assay/fastq/$id.trimmed.2.fq.gz 2\> 1.assay/logs/bowtie2/$id.bowtie2_log \| samtools fastq -N -c 5 -f 12 -F 256 -1 1.assay/bowtie2/$id.rmhost.1.fq.gz -2 1.assay/bowtie2/$id.rmhost.2.fq.gz >sam.shqsub -clear -cwd -l vf=4g,num_proc=4 -P P18Z10200N0127 -binding linear:4  -q st.q -o o.txt -e e.txt sam.sh |awk '{print $3}' 1>>2.runid/sam.qsub.iddone

通过之前的操作获得的数据是样本中物种的·丰富度的信息,接下来就是对数据进行分析试图寻找一些mark或者是在群落相关水平的影响的分析

通常来说会使用多种方法对数据进行处理,得到不同的结果,通过对比几个方案的不同数据,来获取最佳的数据进行后续处理同时也是一个证明数据,为你自己所选择的方法进行作证。

以下是两种方法的数据质量的对比,一种是上面的方法,一种是不知是什么方法,但是样本是一样的作图一般才学加上考研才回来太久没练手了。



T5为上面的方法,发现虽然T5方法unkonw更多,这其实只是数据库的问题,不影响T5方法的准确性,并且T5获取的菌种更多。使用T1方法还有很多未在库中的unkonw细菌没有识别出来。

总结

关于宏基因组分析把握不到位还有很多地方要进步,后期对种群数据分析绘图能力还要加强。争取除夕前学会使用ggplot。以及基础的分析方法。

第一次基因数据处理从集群到数据处理结果---构建菌群物种丰度的图谱相关推荐

  1. 一周碎碎念,2021.11.7,两个MGR集群间还可以构建传统的主从复制通道吗

    叨叨最近遇到的一些事以及见闻.思考. 1. GreatSQL编译环境Dockerfile更新了 用于构建GreatSQL编译环境的Dockerfile发现几个小瑕疵,于是更新了下. 利用Docker环 ...

  2. spark 广播变量大数据_大数据处理 | Spark集群搭建及基本使用

    点击蓝字关注我 前面用了一篇文章详细的介绍了集群HDFS文件系统的搭建,HDFS文件系统只是一个用于存储数据的系统,它主要是用来服务于大数据计算框架,例如MapReduce.Spark,本文就接着上一 ...

  3. Java集群优化——dubbo+zookeeper构建高可用分布式集群

    为什么80%的码农都做不了架构师?>>>    不久前,我们讨论过Nginx+tomcat组成的集群,这已经是非常灵活的集群技术,但是当我们的系统遇到更大的瓶颈,全部应用的单点服务器 ...

  4. linux 计算集群搭建,使用centos构建服务器计算集群

    yum源的搭建 排除列表 #新建文件 [root@centos-test ~]# vim /usr/local/local_mirror/exclude.list #填入内容 SRPMS aarch6 ...

  5. Redis集群进阶之路

    Redis集群规范 本文档基于Redis 3.X或更高版本,讲解Redis集群算法以及设计原理.此官方文档长期更新且随着Redis新版本特性的变化变动,详细请留意官网. 官网地址:https://re ...

  6. 集群概述及原理笔记(1)

    it你好linux学习文档之集群概述及原理笔记(1) 一 前言 目前,越来越多的网站采用Linux操作系统,提供邮件.Web.文件存储.数据库等服务.也有非常多的公司在企业内部网中利用Linux服务器 ...

  7. Redis面试连环问:集群、复制以及与其他NOSQL数据库的区别?

    欢迎关注方志朋的博客,回复"666"获面试宝典 前言 面试官问:您熟悉redis?搭建过redis集群?rediscluster底层的主从复制rdb如何工作的?-- 接下来,就是一 ...

  8. Zookeeper源码分析:集群模式启动概述

    参考资料 <<从PAXOS到ZOOKEEPER分布式一致性原理与实践>> zookeeper-3.0.0 Zookeeper概述 Zookeeper是一个分布式的,开放源码的分 ...

  9. Hadoop集群部署模型纵览1

    vSphere Big Data Extensions(简称BDE)支持多种部署方式来构建Hadoop集群.按: 存储/计算绑定模型:将存储节点(Data Node)和计算节点(Task Tracke ...

  10. 基于Infiniband高性能集群硬件配置方案

       摩尔定律的一再验证残酷的揭示了一个现实:速度是技术发展的终极目标.高性能计算领域也是一样,如何使高性能计算平台运行的更快.更高效一直是服务器厂商研究的方向,曙光作为中国高性能计算的领头羊,作为高 ...

最新文章

  1. 在 Linux 上给用户赋予指定目录的读写权限
  2. RESTful再理解
  3. 【译】.NET Core 是 .NET 的未来
  4. java设计模式---备忘录模式
  5. 成对的交换链表的节点 Swap Nodes in Pairs
  6. 1365. How Many Numbers Are Smaller Than the Current Number 有多少小于当前数字的数字
  7. python批量修改word特定位置的内容_利用python批量修改word文件名的方法示例
  8. .NET Core 如何生成信用卡卡号
  9. $(this).attr(checked, true); 设置不了
  10. WIn server 2012 Install
  11. 函数名、闭包及迭代器
  12. python生成动态链接二维码
  13. 台式计算机能不能安装蓝牙驱动,win7电脑蓝牙驱动怎么安装,详细教您怎么安装...
  14. Python爬虫从入门到精通——爬虫基础:HTTP基本原理
  15. 安卓修改软件安装时间_LaTeX软件安装教程
  16. 项目管理之团队与团队精神
  17. uni-app 手机端定位问题
  18. 浅谈XPS文件格式。
  19. Android 开源图片裁剪工具、图片显示工具分享
  20. 成都拓嘉辰丰:拼多多关键词搜不到自家商品怎么回事?

热门文章

  1. Simulink提速方式
  2. 论坛介绍 | COSCon'22  大数据(D)
  3. 【BDTC 2016】网络与通讯大数据论坛:大数据价值及合规性探索
  4. Docker进阶实战
  5. 【备忘】修复过的SCCM报错
  6. 一款基于java开发的开源监控平台
  7. 初步搭建 prometheus+ Grafana服务器性能监控平台
  8. 将越狱进行到底 Pod2g邀约众大神组建evad3rs
  9. SMETA验厂咨询,Sedex验厂工厂的自检流程有哪些
  10. 如何提高人际交往的沟通技巧