Problem with genome assembly de novo on STACKS
Hi! I realized my genome assembly de novo in STACKS and I used this command to called the SNPs in VCF format.
populations -P ../stacks -M ../stacks/popmap.lst --fstats --vcf --genepop --threads 80
I have a problem with the VCF, I don´t have any value of quality (column "QUAL"), I don´t know what´s wrong.
Thank you so much! and please feel free to ask me any additional info.
• 1,219 views
•
link
written
by
Pueblo errante •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to Perform Lineage Analysis Using SNPs Derived from a De Novo Assembly?
written by alex •
-
Converting custom VCF to standard VCF
written by avelarbio46 •
-
Only output SNPs in bcftools
written by NewtoBioinfo •
-
Determine level of heterozygosity
written by max_19
-
DE novo call from trio
written by cristina_sabiers
-
breseq doesn't identify true SNPs when mapped a Wild Type genome with the Mutant reads and vice ver…
written by ashishbeck96 •
-
Comparing vcf files
written by paraskevopou •
-
Convert chromosome headings in VCF file from Trinity format to chr:pos format
written by TrentGenomics •
-
How to differentiate Homozygous reference calls from No calls in VCF
written by second_exon
-
Filter VCF on Number of genotyped individuals
written by mek362 •
Please provide more information. From the description it sounds like you are doing de novo genome assembly, but I can only assume you mean you are de novo assembling your RAD-Seq type reads with the de novo pipeline in STACKS (also not sure of the version of STACKS you are using or what commands you used to get to the
populationsstep - those need to be provided too if you want any additional help). It is also very possible that there is noQUALvalue for STACKS's VCF output. You could filter on other values such asDP, etc.