Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: bcftools - creating indels
0
gravatar for rayzghimire
4 months ago by
rayzghimire • 0
rayzghimire • 0 wrote:

Hi, So I am converting fastq files into vcf. I wanted to know while calling variants how to include both SNPs and indels using samtools or bcftools.

Thank you

sequencing snp indels • 161 views
ADD COMMENT • link •
modified 4 months ago by Istvan Albert ♦♦ 86k • written 4 months ago by rayzghimire • 0
0
gravatar for Istvan Albert
4 months ago by
Istvan Albert ♦♦ 86k
University Park, USA
Istvan Albert ♦♦ 86k wrote:

Your verbiage seems a bit all over.

You can't convert FASTQ to VCF - it is not a conversion problem.

Perhaps you mean to align FASTQ to a reference, create a BAM alignment file, then call SNPs from the BAM file with bcftools.

If so you don't need to specifically ask for indels, they would be produced along the way.
ADD COMMENT • link written 4 months ago by Istvan Albert ♦♦ 86k

Hi Albert, You are exactly right, that is what I wanted to do. Sorry should have been more clear. A little follow up on that- $ bcftools call --ploidy 1 -m -v -o results/bcf/SRR2584866_variants.vcf results/bcf/SRR2584866_raw.bcf I am using this command to call SNP is there another command in bcftools to call indels?

ADD REPLY • link written 4 months ago by rayzghimire • 0

when we call a tool a SNP caller it almost always also calls insertions and deletions as well.

SNP calling is technically a misnomer. Most variation callers are SNV callers (bcftools included) and are able to produce all small nucleotide variations including insertions and deletions.

ADD REPLY • link written 4 months ago by Istvan Albert ♦♦ 86k
Please log in to add an answer.

Similar posts • Search »

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