This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ambiguous position in SNPdb file and reference genome

Hi,

I have problem with ambiguous base in SNPdb (BQSR) file and reference genome. The first 2 line of SNPdb (BQSR) looks like:

CHROM   POS ID  REF ALT QUAL    FILTER  INFO

1   162362  rs477108567 G   A   .   .   .

so that seems in position 162362 in 1 chromosome in my reference genome I should expect G, but there is A. What is a problem. Those 2 files represent the same reference. Any ideas?

Thank you in advance.

snp reference ngs
• 1,747 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by misterie misterie

SNPdb file is ARS1.2PlusY_BQSR.vcf.gz, reference genome is ARS-UCD1.2_Btau5.0.1Y.fa

• 1 views • link
written by misterie misterie

1 answer

REF Base chr1:162362 is a 'G'

$ curl -s -L "https://sites.ualberta.ca/~stothard/1000_bull_genomes/ARS-UCD1.2_Btau5.0.1Y.fa.gz"  |\
  gunzip  -c | grep -v '>' | fold -w 1 | nl | grep -m 1 -w 162362

 162362 G

REF base in VCF is G and alternate base is A.

I don't see any problem.

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this