Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

RNAseq Library to SNP Call using GATK

$
0
0

Hi ,

I am working on SNP detection on multiple pair end RNAseq library:

after indel realignment, I want perform base quality score recalibration by using these two steps, I do not have any information about known SNP,INDEL.VCF. If I do not provide this option is it ok?

java -jar GenomeAnalysisTK.jar \
-T BaseRecalibrator \
-R refgenome.fasta\
-knownSites known_snps_indels.vcf \
-I sample1.sorted.dedup.realigned.fixmate.bam \
-o sample1.sorted.dedup.realigned.fixmate.recal_data.table \
-cov ReadGroupCovariate \
-cov QualityScoreCovariate \
-cov CycleCovariate

java -jar GenomeAnalysisTK.jar \
-T PrintReads \
-R refgenome.fasta \
-BQSR sample1.sorted.dedup.realigned.fixmate.recal_data.table \
-I sample1.sorted.dedup.realigned.fixmate.bam \
-o sample1.sorted.dedup.realigned.fixmate.recal.bam


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>