I am a high schooler new to R trying to do a bioinformatics project using FASTA files. I want to find a quantitive value for the difference between my test and reference sequence so I can run a chi square test and determine if the difference is statistically significant.

When I used the NCBI BLAST for global align, my sequences were too large. Are there any R packages that would be recommended to compare FASTA/whole genome files?