Mappability calculation based on 150 bp reads after mapping with bwa · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Mappability calculation based on 150 bp reads after mapping with bwa

Hi,

I am trying to apply some filters on whole exome sequencing data. Firstly I did the mapping using bwa and then I followed the proposed pipeline from GATK for Calling variants on cohorts of samples using the HaplotypeCaller in GVCF mode. As next step I am trying to apply the mappability filter <1 (based on 150bp reads) but I have confused with the term "mappability" . Can anyone tell me how can I calculate it ?

Thanks a lot, Maria

mappability we calculate cab score how

• 936 views • link

written 4.8 years ago by astrinaki_maria •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

GATK FilterMutectCalls
written 7 months ago by sk_24 •

I am calling somatic variants using Mutect2 in paired mode. After variant calling, I applied FilterMutectCalls, followed by annotation using Funcotator. I wanted to ask …
Seeking Advice on Multisample Variant Calling Pipeline with BWA-MEME and GATK DRAGEN Mode
written 2.0 years ago by George •

Hello everyone, I am in the process of building a pipeline for multisample variant calling and have some questions regarding tool usage and their safety. …
The GATK "The given bam input has no sample names." error
written 3.0 years ago by Begonia_pavonina •

for f in MINIMAP BWA ; do ~/gatk-4.2.0.0/gatk HaplotypeCaller --reference /home/tmichel/projects/rbge/HybSeq_thibauld/reference_genomes/Begonia_loranthoides_scaffold.fasta --input Hillebrandia_sorted.$f.bam --output Hillebrandia.$f.g.vcf.gz --emit-ref-confidence GVCF ; done I have used GATK to call …
GATK haplotypeCaller output MQ in vcf for non variant sites
written 3.6 years ago by octpus616 •

Hi I am using following command to run `gatk haplotypeCaller` ``` {params.gatk_path} {params.gatk_params} HaplotypeCaller \ -R ${{ref}} \ -I {input.bam} \ -O {output.gvcf} \ -ERC …
gVCF size human genome
written 5.1 years ago by quentin54520 •

Hello all, I got some gvcf using haplotypecaller in gvcf mode (-ERC GVCF). In input I am using human genomes, at around 35 / 40x. …
How to convert gvcf file to normal vcf file ?
written 5.4 years ago by DareDevil

Hi, I followed following steps to arrive at gvcf file: #Mark duplicates in the aligned file gatk MarkDuplicates -I s1.bam --METRICS_FILE s1_markdups_gatk_metrics.txt -O s1_sorted_gatkMarkDups.bam --VALIDATION_STRINGENCY …
monozygotic twins ped file specification
written 5.7 years ago by cocchi.e89

Hi, I got some monozygotic twins samples (and a technical replicate of the same sample sequences twice) and I am running a `GATK GenotypeGVCFs` on …
How to force GATK haplotypecaller to output all genomic position (even non variant positions)
written 6.3 years ago by Nicolas Rosewick

Hi, Is there a way to force GATK haplotypecaller to output all genomic positions in the vcf; not only variant ones ? Example here : …
GATK: HaplotypeCaller gVCF and multisample
written 11.6 years ago by iraun

Hi all, Anyone can explain me what is the main difference between using GATK HC in gVCF mode instead of in multi-sample mode? I know …
Gatk - Haplotypecaller Is So Slow, What Is Faster And As Good?
written 13.0 years ago by newDNASeqer

<p>I have 15 exome-seq samples, and have been using <a href='http://bio-bwa.sourceforge.net/'>BWA</a>-PiCard-GATK pipeline to do the variant calling. I did not realize GATK is so slow …