Freebayes Variant calling
1
0
Entering edit mode
8.8 years ago

Hi, All I'm New to NGS analysis, i want to compare a variant's from different aligner (BWA,Bowtie,Novoalign) Freebayes variant calling program. while doing for Indel i found the command in in freebayes which used for Ignoring SNP alleles.

The command and the argument i have used in given below

./Tools/freebayes/bin/freebayes -f Reference/Reference/human_g1k_v37.fasta -I -X -u -v Output/Novoalign/Variant/Novoalign_InDel_freebayes.vcf Output/Novoalign/Alignment/Novoalign_realignment_recal.bam

Now if i run this my program runs but it stops at certain position saying the following message.

variant at 1:32170247

deletion:7:

alt is empty

and i'm getting the vcf file of 1 mb till the exact position in chromosome 1

If i run the arguments which Ignore insertion and deletion alleles. I'm getting an file size of 1gb with all the chromosome.

I dont know why it stuck at particular postion the same follows for different aligner such as BWA or Novoalign or bowtie it stuck at certain position and ends the program

Does any one knows why ?

Freebayes NGS • 3.9k views
ADD COMMENT
0
Entering edit mode
8.8 years ago
bedeabc ▴ 110

I'd consider cross posting to the Google Group https://groups.google.com/forum/#!forum/freebayes

ADD COMMENT

Login before adding your answer.

Traffic: 1823 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6