Get RMS MQ for specific POS

Hi,

I have a variat calling vcf files with a lots of snp and its bam files, can I calculate MQ (RMS MQ) for each position in vcf files or/and each position in bam file?

Regrads

Zhang

vcf bam mq

• 597 views • link

written 3.4 years ago by octpus616 •

0 answers

No answers yet.

More posts like this

how simply mutation annotation in R?
written 2.9 years ago by octpus616 •

Hi, I am working with a pre-exists mutation list in `Granges`, this data have some fields `chr`, `position`,`ref`, `alt`, `strand`, because our positions are limited …
Why is NA18498 sample missing from 1000 genome GRCh38 vcf file?
written 3.3 years ago by octpus616 •

Hi I am trying to compare the 1000 genome phase3 GRCh19 vcf file and the GRCh38 （low coverage） vcf file. I noticed that one sample, …
rtracklayer import.chain not accpect chain file
written 3.3 years ago by octpus616 •

Hi, I am trying to import chain with `import.chain()` in `rtracklayer` ``` library(rtracklayer) chain <- '/path/to/UCSCliftOver/liftOverChain/hg19-chm13v2.over.chain' ch <- import.chain(chain) ``` ``` Error in .local(con, format, …
what is IGV visualization lighter color variants stand for
written 3.4 years ago by octpus616 •

Hi I am recently trying to use `IGV` to visualize variations in vcf file, I noted some variants mark as a lighter color, what is …
bcftools merge without multiallelic sites
written 3.4 years ago by octpus616 •

Hi I have tried to use `bcftools merge` to merge two `vcf` files ``` # vcf a chr1 10000 A T ..... 0/0 chr1 10003 …
Annotate a series of genomic intervals using R as detailed as possible
written 3.5 years ago by octpus616 •

HI, I plan to annotate a series of genomic intervals using R. I am interested in their location on the genome and the corresponding features, …
Calculation method of RMSMappingQuality
written 3.6 years ago by octpus616 •

HI, I am trying calculate RMS (root-mean-square) MQ (MappingQuality) from a bam file for ever sites (even no variation), I visited GATK's blog: [https://gatk.broadinstitute.org/hc/en-us/articles/360037591751-RMSMappingQuality](https://gatk.broadinstitute.org/hc/en-us/articles/360037591751-RMSMappingQuality) I …
Genotype SNP with IGV
written 5.6 years ago by biom.andressa •

Hi I have bam files of WXS and I'm visualizing it with IGV. Can I genotype a SNP only by visualizing its reads and bases …
SNP calling with SAMtools on multiple samples
written 7.4 years ago by madeira.ag •

Hello. I found similar questions here, but couldnt find exactly what I am searching for. I have several (+180) .bam files, each from a different …
IGV and Variant Call
written 11.2 years ago by a_biologist •

Hello guys, I have two bam files that I load in IGV. Each of those bams I did Variant Calls on them. I did take …