#Download using command line "NC_005956.1" / "CP003784.1"
esearch -db nucleotide -query "NC_005956.1" | efetch -format fasta > seq1.fa
esearch -db nucleotide -query "CP003784.1" | efetch -format fasta > seq2.fa
#fastANI for alignment free comparision
#conda install -c bioconda fastani
fastANI -q seq1.fa -r seq2.fa --visualize -o fastani.out
#Plot using sequences
Rscript plotSeq.R seq1.fa seq2.fa fastani.out.visual