This is a test version of Biostars. For the public version, visit https://www.biostars.org.
can variants be called between two genomes?

Usually, SNP calling is between a reference genome sequence (fasta) and reads (fastq), is it possible to use the commonly used pipeline (BWA alignment and GATK) to call variants between two genomes (fasta files)? For example, I need to find out the SNPs between bacterial strains, which genome sequences have been fully assembled and available. Thank you for your suggestions.

snp alignment next-gen sequence
• 3,136 views • link
written by Charles Yin Charles Yin

2 answers

You could simulate reads from one of the genomes using wgsim and then use a normal pipeline to call the variants. As an alternative, there is a tool called LAST that was built specifically for a genome-genome comparison (I have never used this, though).

• 1 views • link
written by Lisa Ha Lisa Ha

Thank you for your solution!

• 1 views • link
written by Charles Yin Charles Yin

Log in to answer this question.

More posts like this