Question: Count covariates ERROR
0
gravatar for cristina_sabiers
16 months ago by
Spain
cristina_sabiers40 wrote:

WELL..IM GOING TO DIE XD

This pipeline should really need to be updated.

http://seqanswers.com/wiki/How-to/exome_analysis

Im now at step, Quality score recalibration , This step didnt work for me

1) Count covariates:

java -Xmx4g -jar /home/cri/Desktop/GATK/GenomeAnalysisTK.jar \
-l INFO \
-R hg19.fa \
--DBSNP dbsnp132.txt \
-I input.header.realigned.bam \
-T CountCovariates \
-cov ReadGroupCovariate \
-cov QualityScoreCovariate \
-cov CycleCovariate \
-cov DinucCovariate \
-recalFile input.recal_data.csv

so instead I used this one

 java -jar /home/cri/Desktop/GATK/GenomeAnalysisTK.jar \
   -T BaseRecalibrator \
   -R hg19.fa \
   -I input.header.realigned.bam \
   -knownSites /media/cri/CRIS_DATA/All_20160407.vcf \
   -o recal_data.table

wget -c ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/GATK/All_20160407.vcf.gz

sadly got a new ERROR...

ERROR MESSAGE: The platform (Ion Torrent) associated with read group GATKSAMReadGroupRecord @RG:1 is not a recognized platform. Allowable options are ILLUMINA,SLX,SOLEXA,SOLID,454,LS454,COMPLETE,PACBIO,IONTORRENT,CAPILLARY,HELICOS,UNKNOWN

samtools view -H input.header.realigned.bam | grep '^@RG'
@RG ID:1    PU:RUN  LB:LIBRARY  SM:SAMPLE   CN:BCM  PL:Ion Torrent

Is because PL is Ion Torrent instead IONTORRENT?

Any suggestion how I can fix this issue?

Thanks

genome • 503 views
ADD COMMENTlink modified 9 days ago by Biostar ♦♦ 20 • written 16 months ago by cristina_sabiers40
2
gravatar for russhh
16 months ago by
russhh2.7k
UK, U. Glasgow
russhh2.7k wrote:

The given error message would go away if you use IONTORRENT. Whether any further errors will crop up, I'm unsure...

ADD COMMENTlink written 16 months ago by russhh2.7k

Thanks Russhh, I changed PL for IONTORRENT...Is working now :)

ADD REPLYlink written 16 months ago by cristina_sabiers40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1468 users visited in the last hour