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
• 67 views
•
link
0 answers
No answers yet.
Log in to answer this question.
show us an example of mnv that should have been atomized please.