<p>First focus on calling SNPs accurately. See the <a href="http://bio-bwa.sourceforge.net/" rel="nofollow noreferrer">BWA</a>, <a href="https://www.htslib.org/" rel="nofollow noreferrer">SAMtools</a>, <a href="https://samtools.github.io/bcftools/" rel="nofollow noreferrer">BCFtools</a> pipeline (as <a href="https://www.ebi.ac.uk/sites/ebi.ac.uk/files/content.ebi.ac.uk/materials/2014/140217_AgriOmics/dan_bolser_snp_calling.pdf" rel="nofollow noreferrer">described here</a>) or <a href="https://gatk.broadinstitute.org/hc/en-us" rel="nofollow noreferrer">GATK4</a> (as <a href="https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/" rel="nofollow noreferrer">described here</a>).</p>
<p>Once you have the [BV]CFs, then just code up the analysis to see changes over generations.</p>

CC BY-SA 4.0
2021-11-08T12:32:28.993Z
2021-11-08T12:32:28.993Z
Dan Bolser
<p>I have a population in the hundreds and did an Illumina pool-seq (ended up with ~100bp reads for the whole population).</p>
<p>I have the reference genome for the species and now need to align all the population reads to the reference genome, figure out which are SNPs and finally determine which ones are changing from one generation to the next (about 20generations difference).</p>
<p>Right now I'm trying to use popoolation2 <a href="https://sourceforge.net/p/popoolation2/wiki/Main/" rel="nofollow noreferrer">https://sourceforge.net/p/popoolation2/wiki/Main/</a>
but it seems to be a little old as it has not been updated in ~10 years.</p>
<p>In Galaxy, I couldn't find any pipelines dedicated to determine population SNP variation through generation.</p>
<p>So two questions:</p>
<ul>
<li>what software to use to determine significant variation in SNPs in a population (in poolSeq)?</li>
<li>can you suggest me a pipeline for this? if one exists.</li>
</ul>

<snp><phylogenetics>
What to use to determine SNPs in a population? How to find significant SNP variation between generations?
# What to use to determine SNPs in a population? How to find significant SNP variation between generations?

(Asked by: username_0 on 2021-11-02)

I have a population in the hundreds and did an Illumina pool-seq (ended up with ~100bp reads for the whole population).

I have the reference genome for the species and now need to align all the population reads to the reference genome, figure out which are SNPs and finally determine which ones are changing from one generation to the next (about 20generations difference).

Right now I'm trying to use popoolation2 <https://sourceforge.net/p/popoolation2/wiki/Main/>
but it seems to be a little old as it has not been updated in ~10 years.

In Galaxy, I couldn't find any pipelines dedicated to determine population SNP variation through generation.

So two questions:

* what software to use to determine significant variation in SNPs in a population (in poolSeq)?
* can you suggest me a pipeline for this? if one exists.

**Tags:** snp, phylogenetics

---

## 1 answers

### Answer by username_1 on 2021-11-08. Score: 1

First focus on calling SNPs accurately. See the [BWA](http://bio-bwa.sourceforge.net/), [SAMtools](https://www.htslib.org/), [BCFtools](https://samtools.github.io/bcftools/) pipeline (as [described here](https://www.ebi.ac.uk/sites/ebi.ac.uk/files/content.ebi.ac.uk/materials/2014/140217_AgriOmics/dan_bolser_snp_calling.pdf)) or [GATK4](https://gatk.broadinstitute.org/hc/en-us) (as [described here](https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/)).

Once you have the [BV]CFs, then just code up the analysis to see changes over generations.


---

