This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Why were no loci realigned after bcftools normalize?

Why didn't any loci get realigned after using bcftools normalize

Input:

bcftools norm -a --atom-overlaps '*' -f $REF selectionvariants_Qfilt_noStJohn.recode.vcf -o selectionvariants_prim_norm_noStJohn.vcf

part of the script output:

bcftools 1.16
Using htslib 1.21

VCFtools - 0.1.14
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
    --vcf /carc/scratch/projects/mosborne/mosborne2016445/spinedace_genomics/FinalTotalRawVariants.vcf
    --keep keep_no_StJohn.txt
    --recode-INFO-all
    --max-missing 0.5
    --out variants50_noStJohn
    --recode

Keeping individuals in 'keep' list
After filtering, kept 446 out of 475 Individuals
Outputting VCF file...
After filtering, kept 5356785 out of a possible 6799531 Sites
Run Time = 4849.00 seconds

VCFtools - 0.1.14
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
    --vcf variants50_noStJohn.recode.vcf
    --recode-INFO-all
    --mac 3
    --maf 0.02
    --minDP 5
    --min-meanDP 12
    --out selectionvariants_Qfilt_noStJohn
    --recode

**After filtering, kept 446 out of 446 Individuals
Outputting VCF file...
After filtering, kept 16856 out of a possible 5356785 Sites
Run Time = 387.00 seconds
Lines   total/split/realigned/skipped:  24440/0/0/0**

VCFtools - 0.1.14
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
    --vcf selectionvariants_prim_norm_noStJohn.vcf
    --recode-INFO-all
    --out selectionSNPs_noStJohn
    --recode
    --remove-indels
data genetic bcftools snpfiltering normalize

show us an example of mnv that should have been atomized please.

0 answers

No answers yet.

Log in to answer this question.