|
vor 5 Jahren | |
---|---|---|
tasks | vor 6 Jahren | |
inputs | vor 6 Jahren | |
readme.md | vor 5 Jahren | |
workflow.wdl | vor 6 Jahren |
Author : Zhihui Li
E-mail:18210700119@fudan.edu.cn
Git: http://choppy.3steps.cn/lizhihui/fastqscreen.git
Last Updates: 28/08/2019
FastQ Screen是通过定量比对到一组参考基因组上的reads比例,验证DNA样品来源的工具。它把把常见的基因组建个索引,然后逐个比对,哪个比对最高,就是哪个污染。当需要确认来源基因组时,需要针对多个基因组进行比对。它通过定量比对到一组参考基因组上的reads比例,验证DNA样品来源的工具。FastQ Screen旨在常规地用作一种质量控制措施和分析DNA来源不确定或具有多种来源的样品。针对多个基因组进行比对对检测污染或鉴定样品互换也是有用的,若不检测,可能导致错误的实验结果。fastqscreen
是用于 Choppy-pipe 系统使用的 APP。
Tips:2020.2.19日更新,目前fastqscreen不再比对mapping上老鼠上的基因比例,无需重装APP。
在终端中输入以下命令即可快速安装本APP。
1.安装
$ source activate choppy-py3
$ choppy install lizhihui/fastqscreen
$ choppy apps
2.使用
$ choppy samples fastqscreen-latest --out Projectname_fastqscreen_date_people.csv
$ choppy batch fastqscreen-latest Projectname_fastqscreen_date_people.csv --project-name Projectname_fastqscreen_date_people
按照上述步骤安装成功之后,可以通过下面简单的命令即可使用APP:
# Generate samples file
$ choppy samples fastqscreen-latest --out Projectname_fastqscreen_date_people.csv
Projectname_fastqscreen_date_people.csv
包含以下几个需要填写的参数:
read,sample_id
# read 双端测序数据的R1、R2端在阿里云上的路径信息
# sample_id 每个样本任务的识别码。注意:同一个samples文件中,不同样本的ID应该不同
在配置好Projectname_fastqscreen_date_people.csv
文件后,使用以下命令可以提交计算任务:
$ choppy batch fastqscreen-latest Projectname_fastqscreen_date_people.csv --project-name Projectname_fastqscreen_date_people
提交成功后,即可在工作目录下找到生成的目录名为Projectname_fastqscreen_date_people,里面包含了本次提交任务的所有样本信息。
任务成功结束后,便可以在阿里云相应的OSS端生成相应的结果文件。包括测序原始数据对多个基因组进行比对产生的包含来源基因组信息的html
、png
、txt
文件。
FastQ Screen把常见的基因组建个索引,然后逐个比对,哪个比对最高,就是哪个污染。是高通量测序分析中常见的工具,用于评估样品是否被其他物种污染,常见的有支原体污染,接头污染等,一般会整合到分析流程中,作为质控的一部分。
运行APP后,
每个sample对应一个文件夹,内部结构如下:
call-fastq_screen
_screen.html
_screen.png
_screen.txt
fastqscreen: 0.12.0