本篇內(nèi)容主要講解“GATK BQSR的作用是什么”，感興趣的朋友不妨來(lái)看看。本文介紹的方法操作簡(jiǎn)單快捷，實(shí)用性強(qiáng)。下面就讓小編來(lái)帶大家學(xué)習(xí)“GATK BQSR的作用是什么”吧!

超過(guò)10年行業(yè)經(jīng)驗(yàn)，技術(shù)領(lǐng)先，服務(wù)至上的經(jīng)營(yíng)模式，全靠網(wǎng)絡(luò)和口碑獲得客戶，為自己降低成本，也就是為客戶降低成本。到目前業(yè)務(wù)范圍包括了：成都網(wǎng)站制作、成都網(wǎng)站設(shè)計(jì)，成都網(wǎng)站推廣，成都網(wǎng)站優(yōu)化，整體網(wǎng)絡(luò)托管，小程序定制開(kāi)發(fā)，微信開(kāi)發(fā)，APP應(yīng)用開(kāi)發(fā)，同時(shí)也可以讓客戶的網(wǎng)站和網(wǎng)絡(luò)營(yíng)銷和我們一樣獲得訂單和生意！

BQSR 全稱叫做 Base Quality Score Recalibration， 可以理解為堿基質(zhì)量校正。對(duì)于變異位點(diǎn)的鑒定，堿基質(zhì)量是非常重要的。比如測(cè)序識(shí)別到的一個(gè)位點(diǎn)，其堿基和參考基因組上的堿基不同，但是其質(zhì)量值特別低，此時(shí)可以認(rèn)為是一個(gè)測(cè)序錯(cuò)誤，而不是一個(gè)SNP位點(diǎn)。

在測(cè)序的原始數(shù)據(jù)中，本身就提供了每個(gè)堿基對(duì)應(yīng)的質(zhì)量值，但是GATK官方認(rèn)為測(cè)序儀提供的堿基質(zhì)量值，是不準(zhǔn)確的，存在誤差的。

某個(gè)位點(diǎn)前后的堿基的種類，稱之為上下文環(huán)境，會(huì)對(duì)這個(gè)堿基的質(zhì)量值產(chǎn)生影響。對(duì)于A,T,C,G 4種堿基，共有4 x 4 =16 種上下文環(huán)境，左側(cè)的圖是利用fastq文件中測(cè)序儀給出的堿基質(zhì)量值做的圖，可以看到，對(duì)于不同的上下文環(huán)境，堿基質(zhì)量值分布不同；右圖為經(jīng)過(guò)BQSR校正之后，不同上下文環(huán)境中堿基質(zhì)量的分布?？梢钥吹剑Ｕ?，不同的上下文環(huán)境的堿基質(zhì)量分布基本相同。也就是說(shuō)，BQSR消除了上下文環(huán)境對(duì)堿基質(zhì)量的影響。

在堿基質(zhì)量校正時(shí)，主要考慮下列3個(gè)因素：

1. 堿基在reads中的位置

2. 堿基的上下文環(huán)境

3. 堿基原始的質(zhì)量值

根據(jù)這3這個(gè)因素，首先計(jì)算出原始?jí)A基質(zhì)量中錯(cuò)誤的分布模型，然后利用這個(gè)模型對(duì)堿基質(zhì)量校正，生成新的堿基質(zhì)量值。

執(zhí)行BQSR分析包含以下三步

1. 根據(jù)原始bam文件中的堿基質(zhì)量值計(jì)算出系統(tǒng)誤差的分布

命令如下

gatk BaseRecalibrator \
    -R ${ref_fasta} \
    -I ${input_bam} \
    --use-original-qualities \
    -O ${recalibration_report_filename} \
    --known-sites ${dbSNP_vcf} \
    --known-sites ${sep=" --known-sites " known_indels_sites_VCFs}

在計(jì)算的過(guò)程中， 不考慮已知的變異位點(diǎn)的堿基質(zhì)量，--known-sites指定已知變異位點(diǎn)對(duì)應(yīng)的vcf文件。這一步對(duì)單個(gè)樣本進(jìn)行操作，每個(gè)樣本生成一個(gè)錯(cuò)誤模型文件。

2.  綜合多個(gè)樣本的模型，生成一個(gè)總的模型

命令如下

gatk GatherBQSRReports \
    -I ${sep=' -I ' input_bqsr_reports} \
    -O ${output_report_filename}

3. 根據(jù)之前計(jì)算的模型對(duì)堿基質(zhì)量進(jìn)行校正

命令如下：

gatk ApplyBQSR \
    -R ${ref_fasta} \
    -I ${input_bam} \
    -O ${output_bam_basename}.bam \
    -bqsr ${recalibration_report} \
    --static-quantized-quals 10 --static-quantized-quals 20 --static-quantized-quals 30 \
    --add-output-sam-program-record \
    --create-output-bam-md5 \
    --use-original-qualities

BQSR會(huì)對(duì)輸入的bam文件中的堿基質(zhì)量值進(jìn)行替換，替換為校正之后的質(zhì)量值，而原先的質(zhì)量值保存在OQtag 中，示意圖如下

到此，相信大家對(duì)“GATK BQSR的作用是什么”有了更深的了解，不妨來(lái)實(shí)際操作一番吧！這里是創(chuàng)新互聯(lián)網(wǎng)站，更多相關(guān)內(nèi)容可以進(jìn)入相關(guān)頻道進(jìn)行查詢，關(guān)注我們，繼續(xù)學(xué)習(xí)！

本文名稱：GATKBQSR的作用是什么
本文路徑：http://jinyejixie.com/article20/gdjpco.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供小程序開(kāi)發(fā)、自適應(yīng)網(wǎng)站、動(dòng)態(tài)網(wǎng)站、網(wǎng)站設(shè)計(jì)公司、企業(yè)網(wǎng)站制作、搜索引擎優(yōu)化

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請(qǐng)盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng)，如需處理請(qǐng)聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來(lái)源： 創(chuàng)新互聯(lián)