know sites of rice for BaseRecalibrator

how can i get the know sites of the rice for BaseRecalibrator? samtools（mpileup）+bcftools? or gatk Haplotypecaller? finally i want to get the credible snps and indels from my resequence data of rice

dna-seq

• 932 views • link

written 6.3 years ago by 543541656 •

0 answers

No answers yet.

More posts like this

Read group addition order in RNA-seq short variant discovery (SNPs + Indels)
written 9 months ago by iamsmor •

Hello, I’m following the GATK best practices for RNA-seq short variant discovery (SNPs + Indels) and wondering about the correct point to add Read Groups …
Fake indels and artificially low coverage in bcftools mpileup when using close references
written 4.3 years ago by alnam •

Hi, I am working on targetted ONT long reads, mapping them with minimap2 using the following options ~/apps/minimap2/minimap2 -ax map-ont --secondary=no And then doing SNP …
Indel realignment: GATK vs. BCFtools
written 5.8 years ago by martytriggers •

Hello everyone, I'm currently calling SNPs for multiple samples. I have done the pre-processing e.g. trimmed, mapped to reference (BWA mem), added readgroups, marked duplicates, …
is the output.table correct after BQSR?
written 6.3 years ago by 543541656 •

Hi， after I ran the BQSR command gatk --java-options "-Xmx4G" BaseRecalibrator -R /home/riceUsers/fyx/bm4/reference/1_genomic.fna -I 1_dedup.bam --known-sites 1_round1_passed_snp.vcf --known-sites 1_round1_passed_indel.vcf -O 1_round1.table my cotents of the …
How to FLAG low quality SNPs in samtools VCF
written 6.6 years ago by hellbio

Dear all, I have been using the below command to remove low quality SNPs from VCF generated through samtools mpileup and bcftools. > vcfutils.pl varFilter …
get consensus sequence from bam file with insertions
written 7.2 years ago by Vvedenskiia •

Hello everyone! I want to get consensus sequence from my NGS data. I created sorted bam, and visualized it in IGV. I know that alligment …
bcftools 1.4.1 call parameters
written 7.6 years ago by cetin.m

I first do samtools mpileup: samtools mpileup -uf /mnt/NAS/share/Ref/human/hs37d5.fa -l 1000G.bed -q 30 -Q 30 in.bam -o out.bcf Then: I can call snps from the …
DNAseq pipeline steps
written 8.1 years ago by shuksi1984 •

Following are my DNAseq pipeline steps keeping GATKBestPrac in mind: Step1-Quality check of raw data (fastqc) Step2-Alignment with reference genome(bwa) Step3-SAM to BAM conversion and …
direct call of indels/snps with bcftools?
written 10.0 years ago by cristina_sabiers

I am wondering... With vcftools I can call directly my snps and indels from my vcf file, but with bcftools I only found this way …
SAMtools mpileup and bcftools doesn't call InDels...
written 11.1 years ago by Kirill Tsyganov

Hi Guys, I have been researching fair bit over the last few days about SNPs/InDels calling, My project pretty much exactly summarised [here][1]. My NGS …