can variants be called between two genomes?
2
1
Entering edit mode
2.8 years ago
genebow ▴ 160

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 • 1.4k views
ADD COMMENT
0
Entering edit mode
2.8 years ago
Lisa Ha ▴ 90

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).

ADD COMMENT
0
Entering edit mode
2.8 years ago
genebow ▴ 160

Thank you for your solution!

ADD COMMENT

Login before adding your answer.

Traffic: 2404 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6