I have genomes from clinical bacterial infections. I'm trying to find out whether any of the isolates have infected multiple people.
As such I want to find out how many SNPs difference each genome has across the dataset.
There aren't many good quality reference genomes available so individual mapping is hard, relying on the draft assemblies produced within.
Any good tools/pipelines to achieve this? I've used MEGA's pairwise matrix, but I'm dubious about the results. I'm not sure if the quality of the assemblies will have a large impact on the % of distance detected (i.e. if genome A's assemblies is 2.541mb and genome B's assemlies is 2.543, will it score that extra .002 as distance even though it's likely an error)?
•
written
by
f2369583 •
1 answer
I think the fastest way is to use a tool which don't need a reference genome for finding your Snps (like DiscoSnp ++ http://colibread.inria.fr/software/discosnp/) and then compares them to each other for each of your datasets.
Log in to answer this question.
-
Illumina and ONT reads comparison
written by magatam •
-
How evaluate assembly when multiple reference strains are available?
written by MetaPaul •
-
How many SNPs are too many?
written by SushiRoll •
-
Searching for tool to condense ~900 draft genomes from different isolates into one mega sequence
written by Dr. Jason •
-
download all E coli RefSeq genome fasta nucleotide files (NOT plasmids)
written by Siddharth •
-
DE for a de novo species
written by rnaseq2018 •
-
Assembly method validation - am I doing it right?
written by f2369583 •
-
How to use de novo assembly to build onto an existing reference-guided assembly.
written by Pryce Michener •
-
De-Novo Genome Assemblies Comparison
written by Coryza
-
What's a good Variant effect predictor for draft genomes?
written by madkitty
Do you have a VCF file? Look into vcftools.