GATK4--AnalyzeCovariates报错
2019-11-26 本文已影响0人
Ninininiba
以下命令解决了GATK4按照目前GATK3来做AnalyzeCovariates的报错：
“A USER ERROR has occurred: b is not a recognized option”
My commands are as below:
# BQSR - first pass
gatk BaseRecalibrator \
-I=sample_dedup.bam \
-R=/Reference/genome_ref.fasta \
--known-sites mysnp.vcf \
-O=**sample_recal_pass1.tabl**
gatk ApplyBQSR \
-I sample_dedup.bam \
-R /Reference/genome_ref.fasta \
--bqsr-recal-file **sample_recal_pass1.tabl** \
-O **sample_dedup.recal.pass1.bam**
# BQSR - second pass
gatk BaseRecalibrator \
-I **sample_dedup.recal.pass1.bam** \
-R /Reference/genome_ref.fasta \
--known-sites mysnp.vcf \
-O **sample_recal_pass2.tabl**
gatk ApplyBQSR \
-I sample_dedup.recal.pass1.bam \
-R /Reference/genome_ref.fasta \
--bqsr-recal-file **sample_recal_pass2.tabl** \
-O **sample_dedup.recal.pass2.bam**
##Then use the first pass BAM file `sample_dedup.recal.pass1.bam` for downstream processing.
# compare outputs
gatk AnalyzeCovariates \
--before-report-file **sample_recal_pass1.tabl** \
--after-report-file **sample_recal_pass2.tabl** \
--plots-report-file output.pdf
· [Share on Twitter](https://twitter.com/share?url=https%3A%2F%2Fgatkforums.broadinstitute.org%2Fgatk%2Fdiscussion%2Fcomment%2F46067%23Comment_46067) [Share on Google+](https://gatkforums.broadinstitute.org/gatk/post/googleplus/comment?id=46067)
