Hi Biostars,
My advisor wants me to run PSMC on a data set we've generated (because I'm the only one in the lab who can sort of script and use Linux). I've read the manual on the github page, the original paper (and googled for answers), but I fundamentally don't understand the -p parameter or how to estimate the proportion of missing heterozygotes due to low coverage.
If someone has had to run PSMC any tips would be super helpful. I've already got my fastq consensus and can get psmc to run with the example parameters on the github page (which are probably not suitable for my data). Thanks in advance.
•
written
by
memory_donk
1 answer
The first thing I would do is try several rates of -p. How much does it affect your Ne estimates?
The second thing you could do it use the lander waterman to take an education guess at your variant calling dropout:
Last thing to play with is the depth cutoffs.
I've been where you are. Felt your pain.
Log in to answer this question.
-
Estimating per-variant heritability from summary stats
written by Daniel •
-
Create consensus sequence with ambiguity codes using HIFIASM allignment
written by david.lerner •
-
Error in path.expand(file) : invalid 'path' argument
written by Camila Martínez •
-
Extract unmapped sequenced from a sam file output from minimap2
written by ja569116 •
-
Salmon only processes human reads from PDX data
written by jrleary
-
PSMC: fq2psmc conversion problems
written by dthorbur
-
Gene Ontology term enrichment tools that count duplicate/redundant genes
written by memory_donk
-
makeblastdb and blast very slow for 1 Ensembl genome, not for others
written by memory_donk
-
Method for generating a reference-based genome assembly in fasta with no reference genome bases
written by memory_donk
-
Single end or paired end for variant calling
written by memory_donk
