Is it possible to have BWA align contigs to reference and have output visualised via IGV/Tablet with Gaps and SNPs highlighted?
I'm trying this strategy because I need to visualise alignment of 10-20 regions of the reference genome but don't have this regions in separated fasta files for other strains we want to identify SNPs/indels for some reference coordinates. All bacterial data is already assembled.
You may be able to use Mauve for this without doing alignments with bwa.