STAR aligner samtools mpileup C50

I use STAR align different RNA-SEQ sample to ref genome，then use samtools mpileup call snp enter image description here

samtools mpileup -u -q 20 -Q 25 -C50 -b /public/home/zychen/staralignments/rmduplist --VCF --output-tags INFO/AD,AD,DP,SP -f /public/home/zychen/rna-seq/filterfq2085/rice_MC1314_genomes_v7.fasta |bcftools call --keep-alts  --multiallelic-caller  --variants-only -O v - > all_test.vcf

when I did'nt add -C50, MQ will become 20 at every site if I add -C50,MQ will different

rna-seq rna-seq

• 2,394 views • link

updated 9.6 years ago by Biostar

• written 9.7 years ago by zeyuchen2016 •

0 answers

No answers yet.

More posts like this

Issues with GBS pipeline, samtools and bwa specifically
written 1.6 years ago by dimeresia •

I have been working through some GBS data trying to use GB-eaSy pipeline (https://github.com/dpwickland/GB-eaSy). But I am having issues after excluding certain samples from the …
Blog: NGS one-liner to call variants
written 2.8 years ago by barslmn •

Edit 24 Oct 2023: I use nextflow in projects. This is just a curiosity and should be treated as just experimentation. Take a look at …
Bcftools how to add SNP type to the INFO field (DEL,INS,INV,DUP..)
written 4.4 years ago by Michal Nevo •

Hi, I have a VCF file [Which I created with bcftools mpileup| bcftools call ] and I want to make sure that the Bcftools stats …
Filtering for homozygous reference alleles
written 6.9 years ago by AGE •

I would like to obtain a gvcf file for non-variant sites with a minimum coverage of 5x (homozygous reference alleles that I can at least …
Estimating cross contamination in a set of BAMS
written 8.1 years ago by Pierre Lindenbaum

Hi all, I've received a set of BAM files , the variant were called with bcftools ${bcftools_exe} mpileup -Ou -f "${REF}" \ --bam-list "${bam_list}" \ …
strand bias filtering with bcftools
written 9.4 years ago by user31888

I called SNPs with samtools/bcftools from RNASeq data as follows, and then try to keep variants showing a **strand bias p-value < 0.01** (Fischer's exact …