How to visually display the results of blast double sequence alignment using R language package circlize

How to use the R language package circlize to visually display the results of blast double sequence alignment? in view of this problem, this article introduces the corresponding analysis and solution in detail, hoping to help more partners who want to solve this problem to find a more simple and feasible method.

Circlize is a powerful package, and it is very convenient to draw circles in R language.

Today's article records the code that uses circlize as a package to draw a circle diagram to show the result of blast double sequence alignment.

Articles on plant mitochondrial genomes usually analyze gene transfer between organelle genomes, and the basic analysis method is blast comparison. Visual display can be done with this circle diagram.

The first is to use blast to build a library for comparison.

Makeblastdb-in mt.fasta-dbtype nucl-out mt

Blastn-query cp.fasta-db mt-outfmt 6 > output.txt

Then there is the part of preparing data to draw the outermost circle to show the two sequences.

Df df

Chr x y

1 chloroplast 1 0

2 chloroplast 131478 1

3 mitochondrial 1 0

4 mitochondrial 444567 1

Then read the output of blast.

Df1