Question: SNP calling through bcftools
0
gravatar for palak.siddhu
2.0 years ago by
palak.siddhu0 wrote:

Hello, I am doing SNP calling, using the pipeline as follows: 1. index the reference fasta file using bwa. 2. align the read file using bwa aln command. 3. converted the .sai file to .sam file. 4. converted .sam to the .bam file using samtools. 5.constructed mpileup file using samtools. 6.called SNP in the form of .vcf file using bcftools.

but as a result, i m getting empty vcf file.and during the process, there is no error shown by the tools I have used. what may be the reason for this I m not getting. please help.

snp rna-seq software error • 1.5k views
ADD COMMENTlink modified 11 months ago by Lakshman Teja30 • written 2.0 years ago by palak.siddhu0

Steps what you are following looks fine. But if mention what commands you used then it would easy to solve your issue.

bcftools-1.2/bcftools call -Avc input.bcf -O v -o output.vcf

One thing you can just check is the mpileup file and see the columns what it contains for each position manually against your reference (if its small size).

Follow the link and read more about file formats here.

Also, u can use VarScan and know the read count for the reference and variants.

java -jar VarScan.v2.2.11.jar readcounts input.mpileup --min-coverage 20 --min-base-qual 30 --output-file output_readcound.txt

ADD REPLYlink written 2.0 years ago by mks002160
0
gravatar for Lakshman Teja
11 months ago by
BANGALORE, INDIA
Lakshman Teja30 wrote:

Hello, You can use below commands

Generating variation BCF(binary) file for variation

~/Programs/samtools-0.1.17/samtools mpileup -ugf mtb_h37rv.fna aln_sorted.bam | ~/Programs/samtools-0.1.17/bcftools/bcftools view -bvcg - > var.raw.bcf

Converty BCF(binary) variation file to VCF(text) variations file

~/Programs/samtools-0.1.17/bcftools/bcftools view var.raw.bcf > var.raw.vcf

ADD COMMENTlink written 11 months ago by Lakshman Teja30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 852 users visited in the last hour