Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • mcetin
    started a topic Bcftools Mpileup bias for alternative alleles?

    Bcftools Mpileup bias for alternative alleles?

    I have called snps from a bam file using the below commands:

    Code:
    /usr/local/sw/bcftools-1.9/bcftools mpileup -f /mnt/NEOGENE2/share/ref/genomes/hsa/hs37d5.fa -I /mnt/NEOGENE2/share/dna/hsa/comparative_seqs/lazaridis2014/Loschbour.hg19_1000g.dq.bam.hs37d5.fa.cons.90perc.bam -R /mnt/NAS/projects/2018_MCetin_Selection/imputation/Genomic_Hist_SE_Europe/new225_fixed_chr21.bed -q 30 -Q 30 -o Loschbour_ori_chr21_bcftools.bcf
    
    /usr/local/sw/bcftools-1.9/bcftools call -mV indels Loschbour_ori_chr21_bcftools.bcf > Loschbour_ori_chr21_bcftools.vcf

    Then I calculated the percentage of reference and alternative alleles that are in the vcf file. I got

    0/0 10019

    0/1 3412

    1/1 3281

    % ref= 70.1591670656

    % alt= 29.8408329344

    The ratio of alt alleles seem too high to me, since when I do the same thing with GATK Unified genotyper or haplotype caller, I get 5% alternative, which is I think more in line with the expectations from population genetics (Most positions should be reference)

    Has Anyone seen this before?

Latest Articles

Collapse

  • seqadmin
    Genetic Variation in Immunogenetics and Antibody Diversity
    by seqadmin



    The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
    11-06-2024, 07:24 PM
  • seqadmin
    Choosing Between NGS and qPCR
    by seqadmin



    Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
    10-18-2024, 07:11 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 11:09 AM
0 responses
24 views
0 likes
Last Post seqadmin  
Started by seqadmin, Today, 06:13 AM
0 responses
20 views
0 likes
Last Post seqadmin  
Started by seqadmin, 11-01-2024, 06:09 AM
0 responses
30 views
0 likes
Last Post seqadmin  
Started by seqadmin, 10-30-2024, 05:31 AM
0 responses
21 views
0 likes
Last Post seqadmin  
Working...
X