Skip to contents

Compare refseq slots between two phyloseq objects

Source: R/compare_refseq.R
compare_refseq.Rd

lifecycle-experimental

Performs a detailed comparison of the reference sequences (refseq slot) between two phyloseq objects. Identifies shared and unique ASVs/OTUs both by taxa name and by actual DNA sequence content.

This is useful to detect subtle differences between two phyloseq objects that may share the same samples but differ in their ASV/OTU composition, e.g. after different bioinformatics pipelines or filtering steps.

When unique (non-shared) sequences exist, the mean nearest-neighbor k-mer distance from each unique sequence to the other set is computed. This gives a sense of how different the unmatched sequences are from their closest counterpart in the other object.

Usage

compare_refseq(
  physeq1,
  physeq2 = NULL,
  name1 = NULL,
  name2 = NULL,
  k = 5,
  max_seqs = 500,
  seed = NULL,
  verbose = TRUE
)

Arguments

physeq1

(phyloseq or list_phyloseq, required) First phyloseq object, which must have a refseq slot. Alternatively, a list_phyloseq object; in that case the first two phyloseq objects are used (with their names) and physeq2 is ignored.

physeq2

(phyloseq, default NULL) Second phyloseq object. Must have a refseq slot. Ignored when physeq1 is a list_phyloseq.

name1

(character, default NULL) Label for the first phyloseq object in the output. If NULL, inferred from the list_phyloseq names or defaults to "physeq1".

name2

(character, default NULL) Label for the second phyloseq object in the output. If NULL, inferred from the list_phyloseq names or defaults to "physeq2".

k

(integer, default 5) k-mer size for nearest-neighbor distance computation on unique sequences.

max_seqs

(integer, default 500) Maximum number of unique sequences per object to use for distance computation. If there are more, a random sample is drawn. Set to Inf to use all (may be slow).

seed

(integer, default NULL) Random seed for reproducible sampling when max_seqs is exceeded.

verbose

(logical, default TRUE) If TRUE, print a summary of the comparison.

Value

A list of class "compare_refseq" with components:

n_taxa

Named integer vector with the number of taxa in each object.

shared_names

Character vector of taxa names present in both objects.

unique_names_1

Taxa names only in physeq1.

unique_names_2

Taxa names only in physeq2.

shared_seqs

Character vector of unique DNA sequences found in both objects (regardless of taxa name).

unique_seqs_1

DNA sequences only in physeq1.

unique_seqs_2

DNA sequences only in physeq2.

same_name_diff_seq

A data frame of taxa with the same name but different sequences across objects. Columns: taxa_name, seq_physeq1, seq_physeq2.

same_seq_diff_name

A data frame of sequences shared between objects under different names. Columns: sequence, name_physeq1, name_physeq2.

mean_nn_kmer_dist_1

Mean nearest-neighbor k-mer distance from sequences unique to physeq1 to the closest sequence in physeq2. NA if no unique sequences.

mean_nn_kmer_dist_2

Mean nearest-neighbor k-mer distance from sequences unique to physeq2 to the closest sequence in physeq1. NA if no unique sequences.

Author

Adrien Taudiere

Examples

# Compare a phyloseq object with itself (all shared)
res <- compare_refseq(data_fungi_mini, data_fungi_mini)
#> == Reference Sequence Comparison ==
#> physeq1 : 45 taxa
#> physeq2 : 45 taxa
#> 
#> -- By taxa name --
#>   Shared : 45 
#>   Only in physeq1 : 0 
#>   Only in physeq2 : 0 
#> 
#> -- By DNA sequence --
#>   Shared : 45 
#>   Only in physeq1 : 0 
#>   Only in physeq2 : 0 
#> 
res
#> $name1
#> [1] "physeq1"
#> 
#> $name2
#> [1] "physeq2"
#> 
#> $k
#> [1] 5
#> 
#> $n_taxa
#> physeq1 physeq2 
#>      45      45 
#> 
#> $shared_names
#>  [1] "ASV7"   "ASV8"   "ASV12"  "ASV18"  "ASV25"  "ASV26"  "ASV27"  "ASV29" 
#>  [9] "ASV32"  "ASV34"  "ASV35"  "ASV41"  "ASV42"  "ASV46"  "ASV47"  "ASV48" 
#> [17] "ASV49"  "ASV50"  "ASV53"  "ASV54"  "ASV58"  "ASV59"  "ASV61"  "ASV62" 
#> [25] "ASV63"  "ASV64"  "ASV67"  "ASV68"  "ASV71"  "ASV72"  "ASV75"  "ASV77" 
#> [33] "ASV82"  "ASV83"  "ASV85"  "ASV91"  "ASV93"  "ASV94"  "ASV99"  "ASV100"
#> [41] "ASV101" "ASV104" "ASV105" "ASV107" "ASV108"
#> 
#> $unique_names_1
#> character(0)
#> 
#> $unique_names_2
#> character(0)
#> 
#> $shared_seqs
#>  [1] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCGCTTTCTTGTGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCAACTTTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTTGCTTCTGTGGAGACACGCTTTCAAACCGTCCGAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                      
#>  [2] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTTGCTTTTTTGGCGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCTACTTCTGGTTGCTCGGCTCCTCTCAAATTCATCAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTCGCTTCTGTTGGAGACACGCTTTCTAACCGTCCCGTAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                              
#>  [3] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCCTTTCTTTATTGAAGGGTGAGCTTGGACTTGGGGATTGCTGGTGCAAATTGGCTTCCCTTGAATGTATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACTGTGTAGCTCGGATTGTCTGGCTTCTAATCGTCCTTCACGGGACAATTTGATCTTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                         
#>  [4] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTTGCTTTTTTGGCGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCTACATTCTGGTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTCGCTTCTGTTGGAGACACGCTTTCTAACCGTCCCGTAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                             
#>  [5] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATGTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCTTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#>  [6] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCACTTTTGTGAACGTAGTGGATTGGACTTGGAGGCTTTGCCGGGCTTCACCGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGTGCTTGCTTCTGTGGAGACGCGCTTTCTAACCGTCCGAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                        
#>  [7] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGGAATTCTCAACTCATTAGATTTTTTTCTAATGCAGCTTGGACTTGGAGGTTGTGTCGGTTTTTAATAATCGACTCCTCTTAAATGCATTAGCGTGAATTCTTACGGATCGCCTTCAGTGTGATAATTATCTACGCTGTGGTGTGAAGTATCATTATGTTCATGCTTATAATCGTCTGTTCATCAGACAATTTATTATGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                          
#>  [8] "AAATGCGATACGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACCTCGTGGCTTACCACGGGGTACGCCCGTTCGAGTGTCACGTTCCCATCGCCTCACGGCGGATCTGGGCGCTGCCTTACGGCTCGCCTTAAATGCATTAGCGAGGACAACCGTGCTTGCGCTGGCGTGATAATTGTCTCGCCTCAGCTCGCATACATGGCCTGCTTCCAATCGCACGGCGTTCCGTGCACTATTCTACTTTTGACCTCGAATCGGGTGGGACTACCCGCTGAACTTAA"                                                                                    
#>  [9] "AATTGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCCTGGTATTCCGGGGAGCACGCCTGTTCGAGTGTCGTGAATACCTCAAGCTGGATGCTTTTTGTTAAGTGTTCGGCTTGGATTTGGACTCTGCCGTAACTGTTTTCTGAACGGGGACGGCTAGTCTTAAATGCATTAGCTGATCCATGCTGTGTGTTTGGTTCTACTCGGCGTGATAATTATCTGACCGCTGAGGACATTTGCTTTGTTGCGAATGGCCAGATGTGCCTTGGGTTGCTTCTAACCTGTCTGCGCGCAAGTGCGACAACTTTTTCAACTTTGACCTCGAATCAGGTGGGACTACCCGCTGAACTTA"    
#> [10] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCTCCTTGGTATTCCGAGGAGCATGTCTGTCTGAGTGTCATGAAATTCTCAACCCTTTGATTTATTTGGAAGGGCTTGGATCATGAGCGCTGCAGGCCCCAGCTGGGGTCAGCTCGCTCGAAATAGATCAGCGGTTAACCTTAGGTTTACTGGCCTGACTCGACGTGGTAGAGTGACTTTATTCGTCGAGGATCCCTGCTTGCGGGGAGCCACTGAACTTTGAGGCCGCTCCTAATTGCGTATGCACTCTTTTGGTTTAGACCTCAGATCAGGCAGGACTACCCGCTGAACTTA"                             
#> [11] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGCGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#> [12] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCATGCCTGTTTGAGTGTCATTAAATTATCAACCTTGCTCGCTTTTACCGGCTTGAGTTAGGCTTGGATGTGAGGGCTTGCCGGCTTCCTTCAGTGGATGGTCAGCTCCCTTGAAATGCATTAGTGGGATCTCTTGCGGACCGTCACTTGGTTGTGATAATTATCTACGCCATCTTGACTTTGAAGCAAACTTATGGGAACCTGCTTATAACCGTCCTTTTCAAAAGGACAATCTTTGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                       
#> [13] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATCTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCCTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#> [14] "AAATGCGATAAGTAATGTGAATTGCAGAACTCAGTGAATCATCGAATCTTTGAACGCATCTTGCACTCCATGGTATTCCGTGGAGTATGCCTGTCTGAGTGTCATGAAACTCTCAACCCCGTTACCTTTGTTTATCTTGGGCATTGGGGCTTGGATCATGAGTGTCTGCTGGACCTGTGCATCAGCTCACTCTAAATGTATCAGTGGTCTGCTTTCTCTGCCTCTGCAGATAGATCTGGCGTAATAGGTTTTATTCGCTGGGACTTCTGTGGGGAACAGCAGGGTGTTGCTTACCATTCACAAACAACTCTAAGTGTTTGGCCTCAGATCAGGCAGGACTACCCGCTGAACTTA"                                 
#> [15] "AAATGCGATAAGTAATGTGAATTGCAGATTTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTCTTTGGATTTTTGTTCATTGAAGACTTGGATGTGAGGGTTTGCCGGTTATTGTATCGGCTCCCTTTAAATGTGTTAGCGTTACCTTCAGCCGATCGAACGCCCCTCGGTATGATAATTATCTATGCCGTGTCGGGGATCGCACGAGAGTATGCTTCTAACAGTCCCTCTGGGGACATATTCATCATATCTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                         
#> [16] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCTCCGTGGTATTCCGCGGAGCATGCCTGTTTGAGTGTCACGTAACCATCGCCCCTCAGATCGTTAACGCGGTCTGTCGAGGTGGATTTGGACTGTGCCGTAGTTGGCTCGTCTTGAAATACATTAGCTTGCGCTTTTAGAGTGCCGGGCACCGGTGTGATAATTATCTGCGCCACAGCCTTAGGAACTCTTCAGTGGTGCTGCTTATAGTCGTCCTTACGGACAATACACTTTAAAGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                              
#> [17] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCCTGGTATTCCGGGGAGCATGCCTGTTTGAGTGTCATGGTATCCTCATCCTTCATAACTTTTTGTTATCGAAGGCATGGACTTGGAGGTCGTGCTGGTTCCTCGTTGAATCGGCTCCTCTTAAATGTATTAGCGTGAGTGTAACGGATCGCTTCGGTGTGATAATTATCTGCGCCGTGGTCGTGAAGTAACATAAGCTTGCGCTTCTAACCGTCCTTAAGCTGGACAACATAACTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                        
#> [18] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCACACCTGTTTGAGTGTCGTGAATCTCTCAACCTCCTCATCTTTTTCAAAGGTGGGCTGAGGTTGGACCTGGAGGTCTTGCAGGCGTTGAGCCTGCTCCTCTCAAATGCATTAGCGGTGCCCAGGCATCCTCGGCGTGATAAGTATCTACGTCTCTGGCATGCTGTCCGGTGGCTGTCGCTTCCAACCGTCTCCTGCGAGACAATGCTTGAAACTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                                   
#> [19] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGTGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#> [20] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCATGGAGCATGCCTGTTTGAGTGTCACGTAAACCCTCACCCTCGCGACGTAACAATCGCCCGCGGTGGACTTGGACCGTGCCGTAATCGGCTCGTCTTGAAATGCATTAGCTGGCGCTTTTAGAGTGTTGGGCGAACGGTGTGATAATTATCTGCGCCGACTAGCTCTGGGCCTCTTCAGCGGCGCCGCTTACAGCCGTCCCTCCGTGGACAACTTACATTAAAGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                             
#> [21] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGCGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#> [22] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGTTATTCTCAGCTCCTCTCAGATTTTTCTGTTGTTGGAGACTGGACTTGGAAGCTGCTGGAACGTCAGGTTCAGCTCTTCTTAAATGCATTAGTGTGAATATTTGTTGAATGTTCTCCAGTGTGATAATTGTCTACGCTGTTTGAACTTTGGCAGAAATTATATATTCATGCTTCAAACCGTCTTCGGACAACTTTTTGCATTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#> [23] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACTCTTTGGTATTCCGAAGAGTATGCCTGTTTGAGTGTCATGTATTCCTCAATCTCTTCAGTTTGTTCTGGAGTTGATTGGAGTTGGAAGCTGCTGGTTTGAGTCTGAATCAGCTCTTCTTAAATGCATTAGTGTGAATGTTACTGAATGTTCTCCAGTGTAATAATTATTTGCGCTGCTAGGACTTTGGTAGAATATGTATATATATTCATGCTTCTTAACCGTCTGTTATGGACAATACACTTTGCGTTTGACCTCAAATCAGGTAGGATTACCCGCTGAACTTA"                                    
#> [24] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACCCTTTGGTATTCCGAAGGGTACGCCCGTTTGAGTGTCATTGTACTCTCACCCTCAAGTTTTCTTGAGAGCGTGGATCTGGACGCTGCCGTGTAAACGGCTCGTCTGAAATGAATGAGTGTAACCTGTATGTCAGTGTTATTCGATGTGATAAGTATCTTCATCGGAGTAGGTTTACCTGCGCTTTCATATTGGCCGCGCTGCCAACCGTCTTCGGACAATCTTGACAATTTGACCTCAAATCGGGTGGGACTACCCGCTGAACTTA"                                                       
#> [25] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCATGCCTGTTTGAGTGTCATGGAATTCTCATCCCCTAAATTTATTTTTGGTGGGCTTGGATGTGGAGGTTGTGCCGGTTTGTTAAAAGCCGGCTCCTCTGAAATGCATTAGTGTGAATCTTACTGATCGCCTTCAGTGTGATAATTGTCTACGCTGTGGTTGTGAAGTAAATCAATGTTCATGCTTCTAATCGTCTCTTTTTGAGACAATTACTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                            
#> [26] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGGAATTCTCAACTTATAAATCCTTGTGATCTATAAGCTTGGACTTGGAGGCTTGCTGGCCCTCGTTGGTCGGCTCCTCTTGAATGCATTAGCTCGATTCCGTACGGATCGGCTCTCAGTGTGATAATTGTCTACGCTGTGACCGTGAAGTGTTTTGGCGAGCTTCTAACCGTCCATTAGGACAACTTTTTAACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                  
#> [27] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCCTTTCTTTATTGAAGGGTGAGCTTGGACTTGGGGATTGCTGGTGCAAATTGGCTTCCCTTGAATGTATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACCGTGTAGCTCGGATTGTCTGGCTTCTAATCGTCCTTCACGGGACAATTTGATCTTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                         
#> [28] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATGTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCCTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#> [29] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACCTTTTGGTATTCCGAAAGGTATGCCTGTTTGAGTGTCTTGAATCATCACTTCAAAGAGTTTTCTTTGAATGTGGATTTGGGTGTTGCCATTTGGCTCGCCTTGAATGTGTTAGTGGGTCAGCCCAAGCTAGGTGTAATAAGTTTCATCTGGTGTAAGTGTCGTCCGCTTCCAACTCGCAAGGCTTTGGCCTGCATTTTCATGACTAGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                                                
#> [30] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCATGCCTGTTTGAGTGTCATGGAATTCTCATCCCCTAAATTTATTTTTGGTGGGCTTGGATGTGGAGGTTGTGCCGGTTTGTTAAAAGCCGGCTCCTCTGAAATGCATTAGTGTGAATCTTACTGATCGCCTTCAGTGTGATAATTGTCTACGCTGTGGTTGTGAAGTAAATTAATGTTCATGCTTCTAATCGTCTCTTTTTGAGACAATTACTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                            
#> [31] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGCATTCCGAGGGGCACGCCTGTTCGAGTGTCGTGAACTCCTCCACTCCTTACCTTTTTCGAAGGGTACGGGGCTGGGATTTGAGAGTTGCGGGTCCCTGGCCGATCCGCTCTCTTTGAATGCATTAGCGAAGCCCTCGCGGCCTTGGTGTGATAATTATCTACGCCTCGGCTTAGCGCATGACCGGGCTATGCTTCGAACCGTCTCGCAAGAGACAACGACTACCAACTTGACCTCGAATCAGGCGGGACTACCCGCTGAACTTA"                                                  
#> [32] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATCTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCTTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#> [33] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCATGGTATTCCGTGGAGCATGCCTGTTTGAGTGTCACGTCAACCCTCATTCATGTGATGTTACAGTCGCATGAATGGATTTGGACCGTGTCGTAATCGACTCGTCCCGAAATGCATTAGCTTGCGCTTTTAGAGTGCTGGGCACCGGTGTGATAATTATCTGCGCCAACGCCTTAAGCCTCTTCAGCGGCGCTGCTTACAGCCGTCCTGACGGACAACTATTTTAACGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                    
#> [34] "AAATGCGATAAGTAATGTGAATTGCAGTATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCCCTTTGGTATTCCGAAGGGCATGCCTGTTTGAGTGTCATGAAAATCTCATCCCATCGGGTTTCCGATGACGATGGACTTGGGTGTTGCCGTCCCGGCTCGCCTTAAAAGTGTTAGTGGTTCCGAGAGGGAATCTGGCGTAATAAGTTTCGTCGGTGAATCTCGAAGGCTGCTCATAACCCGTCTTCTGACATCTTTTTGCTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTAA"                                                                                  
#> [35] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAACATATCAACCCCGACGGGATTTTAATCTCGACGGTGGCTTGGACTTGGGAGCCTGCTGGCCGTTGTATTTCGGTCGGCTCTCCTTGAATGCATTAGCTTGGACTTGTTTCATGGGCTGCATAGCGTGATAACTATCTACGCTTTGTACTGTCTGTGGGATACTCATGAGTCTGGCTTCTAATCGTCCTCTGGACAATTACTTTACGTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                  
#> [36] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCTCCTTGGTACTCCATGGAGCATGCCTGTTTGAGTGTCACGTAAACCCTCACCCTTGCGATGTAACAGTCGTTCGCGGTGGACTTGGACCGTGCCGTAGTTGGCTCGTCCTGAAATGCATTAGCTGGTGCTTTTAGAGTGCTGGGCGACGGTGTGATAATTATCTGCGCCAATGCCCTGGGCCTCTTCAGCGGTGCCGCTTACAGTCGTCCCCAGAGGACAACGTATCTAAAGCTTTGGCCTCAAATCAGGTAGGATTACCCGCTGAACTTA"                                                  
#> [37] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCGTCGAATCTTTGAACGCACCTTGCGCCCTTTGGTGTTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCGCTTTCTTGTGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCAACTTTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTTGCTTCTGTGGAGACACGCTTTCAAACCGTCCGAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                      
#> [38] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCGCTTTCTTGTGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCAACTTTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTCGCTTCTGTTGGAGACACGCTTTCTAACCGTCCCGTAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                   
#> [39] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGTGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#> [40] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTTTCCGGTATTCCGGAAAGCATGCCTGTTTGAGCGCCTTATCATAGTTGTGCAATTGCATTTTTTGTAGTTGCACAGTATTGGGATTTGCTGCGTGAGCAGCTCACCCTGAATGCATTAGCAATTCGACACTTGAAACTAACGTGATAAGTATCTTCGTTGGCTTGAAGTGAAGACTTGCTTCGAATCGTCTTCGGACAATACCTAATGCTAGGCATCAAATCAGGTAGGATTACCCGCTGAACTTA"                                                                          
#> [41] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCTTTGGTATTCCGAAGAGCATGCCTGTTTGAGTGTCATGAATCTCTCAACTCTTATGATCTGTAATGGATCGTATTGAGCTTGGACTTGGAGGCTTGCCGGCGTTCGTCGGCTCCTCTTAAATGCATTAGCTATGGATCACGATGTGTATGATCTACGGTGTGATATATGTCTACGCCTATGGATCTGCGCGCGAGTTACCGCCTGCAAAGGTAGCGCTCTCCGGAGCCGTGTAGCGAGTGAATACGAAGCTTCTAATCGTCTTCGGACAACTTCTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"
#> [42] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCACGCCTGTTCGAGTGTCGTGAAATTCTCAACTCAATCCTCTTGTTATGAGAGGGCTGGGCTTGGACTTGGAGGTCTTGCCGGTTGTGTTCCCTCGGGACGTCGGCTCCTCTTGAATGCATGAGTGGATCCCTTTTTGTAGGGTTTGCCCTTGGTGTGATAATTATCTACGCCGCGGGTAGCCTTGCGCGTTGGTCTGCTTCTAACCGTCCTTCAGGACATTTCATCTCAACTTGACCTCGAATCAGGCGGGACTACCCGCTGAACTTA"                                 
#> [43] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCTATTCTTTATTGAATGGTGAGCTTGGATTTGGAGATTGCTGGTGAAAATCGGCTTCTCTTGAATGAATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACTGTGTAGCTCGGATTGTCTGGCTTCTAATCGTCTTTCATGGACAATTTGATCATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                          
#> [44] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCCTTTCTTTATTGAAGGGTGAGCTTGGACTTGGGGATTGCTGGTGCAAATTGGCTTCCCTTGAATGTATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACTGTGTAGCTCGGATTGTCTGGCTTATAATCGTCCTTCACGGGACAATTTGATCTTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                         
#> [45] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCATGGTATTCCGTGGAGCATGCCTGTTTGAGTGTCACGTCAACCCTCATTCATGTGATGTTACAGTCGCATGAATGGATTTGGACCGTGTCGTAATCGACTCGTCCCGAAATGCATTAGCTTGCGCTTTTAGAGTGCTGGGCACCGGTGTGATAATTATCTGCGCCAACGCCTTAAGCCTCTTCAGTGGCGCTGCTTACAGCCGTCCTGACGGACAACTATTTTAACGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                    
#> 
#> $unique_seqs_1
#> character(0)
#> 
#> $unique_seqs_2
#> character(0)
#> 
#> $same_name_diff_seq
#> [1] taxa_name   seq_physeq1 seq_physeq2
#> <0 rows> (or 0-length row.names)
#> 
#> $same_seq_diff_name
#> [1] sequence     name_physeq1 name_physeq2
#> <0 rows> (or 0-length row.names)
#> 
#> $mean_nn_kmer_dist_1
#> [1] NA
#> 
#> $mean_nn_kmer_dist_2
#> [1] NA
#> 
#> attr(,"class")
#> [1] "compare_refseq"

# Compare with a subset
sub <- prune_taxa(taxa_names(data_fungi_mini)[1:20], data_fungi_mini)
res2 <- compare_refseq(data_fungi_mini, sub,
  name1 = "full", name2 = "subset"
)
#> == Reference Sequence Comparison ==
#> full : 45 taxa
#> subset : 20 taxa
#> 
#> -- By taxa name --
#>   Shared : 20 
#>   Only in full : 25 
#>   Only in subset : 0 
#> 
#> -- By DNA sequence --
#>   Shared : 20 
#>   Only in full : 25 
#>   Only in subset : 0 
#> 
#> 
#> -- Mean nearest-neighbor k-mer distance (k=5) --
#>   Unique in full -> subset : 0.121 
res2
#> $name1
#> [1] "full"
#> 
#> $name2
#> [1] "subset"
#> 
#> $k
#> [1] 5
#> 
#> $n_taxa
#>   full subset 
#>     45     20 
#> 
#> $shared_names
#>  [1] "ASV7"  "ASV8"  "ASV12" "ASV18" "ASV25" "ASV26" "ASV27" "ASV29" "ASV32"
#> [10] "ASV34" "ASV35" "ASV41" "ASV42" "ASV46" "ASV47" "ASV48" "ASV49" "ASV50"
#> [19] "ASV53" "ASV54"
#> 
#> $unique_names_1
#>  [1] "ASV58"  "ASV59"  "ASV61"  "ASV62"  "ASV63"  "ASV64"  "ASV67"  "ASV68" 
#>  [9] "ASV71"  "ASV72"  "ASV75"  "ASV77"  "ASV82"  "ASV83"  "ASV85"  "ASV91" 
#> [17] "ASV93"  "ASV94"  "ASV99"  "ASV100" "ASV101" "ASV104" "ASV105" "ASV107"
#> [25] "ASV108"
#> 
#> $unique_names_2
#> character(0)
#> 
#> $shared_seqs
#>  [1] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCGCTTTCTTGTGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCAACTTTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTTGCTTCTGTGGAGACACGCTTTCAAACCGTCCGAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                  
#>  [2] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTTGCTTTTTTGGCGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCTACTTCTGGTTGCTCGGCTCCTCTCAAATTCATCAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTCGCTTCTGTTGGAGACACGCTTTCTAACCGTCCCGTAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                          
#>  [3] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCCTTTCTTTATTGAAGGGTGAGCTTGGACTTGGGGATTGCTGGTGCAAATTGGCTTCCCTTGAATGTATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACTGTGTAGCTCGGATTGTCTGGCTTCTAATCGTCCTTCACGGGACAATTTGATCTTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                     
#>  [4] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTTGCTTTTTTGGCGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCTACATTCTGGTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTCGCTTCTGTTGGAGACACGCTTTCTAACCGTCCCGTAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                         
#>  [5] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATGTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCTTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                 
#>  [6] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCACTTTTGTGAACGTAGTGGATTGGACTTGGAGGCTTTGCCGGGCTTCACCGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGTGCTTGCTTCTGTGGAGACGCGCTTTCTAACCGTCCGAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                    
#>  [7] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGGAATTCTCAACTCATTAGATTTTTTTCTAATGCAGCTTGGACTTGGAGGTTGTGTCGGTTTTTAATAATCGACTCCTCTTAAATGCATTAGCGTGAATTCTTACGGATCGCCTTCAGTGTGATAATTATCTACGCTGTGGTGTGAAGTATCATTATGTTCATGCTTATAATCGTCTGTTCATCAGACAATTTATTATGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                      
#>  [8] "AAATGCGATACGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACCTCGTGGCTTACCACGGGGTACGCCCGTTCGAGTGTCACGTTCCCATCGCCTCACGGCGGATCTGGGCGCTGCCTTACGGCTCGCCTTAAATGCATTAGCGAGGACAACCGTGCTTGCGCTGGCGTGATAATTGTCTCGCCTCAGCTCGCATACATGGCCTGCTTCCAATCGCACGGCGTTCCGTGCACTATTCTACTTTTGACCTCGAATCGGGTGGGACTACCCGCTGAACTTAA"                                                                                
#>  [9] "AATTGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCCTGGTATTCCGGGGAGCACGCCTGTTCGAGTGTCGTGAATACCTCAAGCTGGATGCTTTTTGTTAAGTGTTCGGCTTGGATTTGGACTCTGCCGTAACTGTTTTCTGAACGGGGACGGCTAGTCTTAAATGCATTAGCTGATCCATGCTGTGTGTTTGGTTCTACTCGGCGTGATAATTATCTGACCGCTGAGGACATTTGCTTTGTTGCGAATGGCCAGATGTGCCTTGGGTTGCTTCTAACCTGTCTGCGCGCAAGTGCGACAACTTTTTCAACTTTGACCTCGAATCAGGTGGGACTACCCGCTGAACTTA"
#> [10] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCTCCTTGGTATTCCGAGGAGCATGTCTGTCTGAGTGTCATGAAATTCTCAACCCTTTGATTTATTTGGAAGGGCTTGGATCATGAGCGCTGCAGGCCCCAGCTGGGGTCAGCTCGCTCGAAATAGATCAGCGGTTAACCTTAGGTTTACTGGCCTGACTCGACGTGGTAGAGTGACTTTATTCGTCGAGGATCCCTGCTTGCGGGGAGCCACTGAACTTTGAGGCCGCTCCTAATTGCGTATGCACTCTTTTGGTTTAGACCTCAGATCAGGCAGGACTACCCGCTGAACTTA"                         
#> [11] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGCGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                       
#> [12] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCATGCCTGTTTGAGTGTCATTAAATTATCAACCTTGCTCGCTTTTACCGGCTTGAGTTAGGCTTGGATGTGAGGGCTTGCCGGCTTCCTTCAGTGGATGGTCAGCTCCCTTGAAATGCATTAGTGGGATCTCTTGCGGACCGTCACTTGGTTGTGATAATTATCTACGCCATCTTGACTTTGAAGCAAACTTATGGGAACCTGCTTATAACCGTCCTTTTCAAAAGGACAATCTTTGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                   
#> [13] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATCTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCCTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                 
#> [14] "AAATGCGATAAGTAATGTGAATTGCAGAACTCAGTGAATCATCGAATCTTTGAACGCATCTTGCACTCCATGGTATTCCGTGGAGTATGCCTGTCTGAGTGTCATGAAACTCTCAACCCCGTTACCTTTGTTTATCTTGGGCATTGGGGCTTGGATCATGAGTGTCTGCTGGACCTGTGCATCAGCTCACTCTAAATGTATCAGTGGTCTGCTTTCTCTGCCTCTGCAGATAGATCTGGCGTAATAGGTTTTATTCGCTGGGACTTCTGTGGGGAACAGCAGGGTGTTGCTTACCATTCACAAACAACTCTAAGTGTTTGGCCTCAGATCAGGCAGGACTACCCGCTGAACTTA"                             
#> [15] "AAATGCGATAAGTAATGTGAATTGCAGATTTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTCTTTGGATTTTTGTTCATTGAAGACTTGGATGTGAGGGTTTGCCGGTTATTGTATCGGCTCCCTTTAAATGTGTTAGCGTTACCTTCAGCCGATCGAACGCCCCTCGGTATGATAATTATCTATGCCGTGTCGGGGATCGCACGAGAGTATGCTTCTAACAGTCCCTCTGGGGACATATTCATCATATCTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#> [16] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCTCCGTGGTATTCCGCGGAGCATGCCTGTTTGAGTGTCACGTAACCATCGCCCCTCAGATCGTTAACGCGGTCTGTCGAGGTGGATTTGGACTGTGCCGTAGTTGGCTCGTCTTGAAATACATTAGCTTGCGCTTTTAGAGTGCCGGGCACCGGTGTGATAATTATCTGCGCCACAGCCTTAGGAACTCTTCAGTGGTGCTGCTTATAGTCGTCCTTACGGACAATACACTTTAAAGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                          
#> [17] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCCTGGTATTCCGGGGAGCATGCCTGTTTGAGTGTCATGGTATCCTCATCCTTCATAACTTTTTGTTATCGAAGGCATGGACTTGGAGGTCGTGCTGGTTCCTCGTTGAATCGGCTCCTCTTAAATGTATTAGCGTGAGTGTAACGGATCGCTTCGGTGTGATAATTATCTGCGCCGTGGTCGTGAAGTAACATAAGCTTGCGCTTCTAACCGTCCTTAAGCTGGACAACATAACTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                    
#> [18] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCACACCTGTTTGAGTGTCGTGAATCTCTCAACCTCCTCATCTTTTTCAAAGGTGGGCTGAGGTTGGACCTGGAGGTCTTGCAGGCGTTGAGCCTGCTCCTCTCAAATGCATTAGCGGTGCCCAGGCATCCTCGGCGTGATAAGTATCTACGTCTCTGGCATGCTGTCCGGTGGCTGTCGCTTCCAACCGTCTCCTGCGAGACAATGCTTGAAACTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                               
#> [19] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGTGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                       
#> [20] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCATGGAGCATGCCTGTTTGAGTGTCACGTAAACCCTCACCCTCGCGACGTAACAATCGCCCGCGGTGGACTTGGACCGTGCCGTAATCGGCTCGTCTTGAAATGCATTAGCTGGCGCTTTTAGAGTGTTGGGCGAACGGTGTGATAATTATCTGCGCCGACTAGCTCTGGGCCTCTTCAGCGGCGCCGCTTACAGCCGTCCCTCCGTGGACAACTTACATTAAAGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                         
#> 
#> $unique_seqs_1
#>  [1] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGCGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#>  [2] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGTTATTCTCAGCTCCTCTCAGATTTTTCTGTTGTTGGAGACTGGACTTGGAAGCTGCTGGAACGTCAGGTTCAGCTCTTCTTAAATGCATTAGTGTGAATATTTGTTGAATGTTCTCCAGTGTGATAATTGTCTACGCTGTTTGAACTTTGGCAGAAATTATATATTCATGCTTCAAACCGTCTTCGGACAACTTTTTGCATTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#>  [3] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACTCTTTGGTATTCCGAAGAGTATGCCTGTTTGAGTGTCATGTATTCCTCAATCTCTTCAGTTTGTTCTGGAGTTGATTGGAGTTGGAAGCTGCTGGTTTGAGTCTGAATCAGCTCTTCTTAAATGCATTAGTGTGAATGTTACTGAATGTTCTCCAGTGTAATAATTATTTGCGCTGCTAGGACTTTGGTAGAATATGTATATATATTCATGCTTCTTAACCGTCTGTTATGGACAATACACTTTGCGTTTGACCTCAAATCAGGTAGGATTACCCGCTGAACTTA"                                    
#>  [4] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACCCTTTGGTATTCCGAAGGGTACGCCCGTTTGAGTGTCATTGTACTCTCACCCTCAAGTTTTCTTGAGAGCGTGGATCTGGACGCTGCCGTGTAAACGGCTCGTCTGAAATGAATGAGTGTAACCTGTATGTCAGTGTTATTCGATGTGATAAGTATCTTCATCGGAGTAGGTTTACCTGCGCTTTCATATTGGCCGCGCTGCCAACCGTCTTCGGACAATCTTGACAATTTGACCTCAAATCGGGTGGGACTACCCGCTGAACTTA"                                                       
#>  [5] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCATGCCTGTTTGAGTGTCATGGAATTCTCATCCCCTAAATTTATTTTTGGTGGGCTTGGATGTGGAGGTTGTGCCGGTTTGTTAAAAGCCGGCTCCTCTGAAATGCATTAGTGTGAATCTTACTGATCGCCTTCAGTGTGATAATTGTCTACGCTGTGGTTGTGAAGTAAATCAATGTTCATGCTTCTAATCGTCTCTTTTTGAGACAATTACTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                            
#>  [6] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATGGAATTCTCAACTTATAAATCCTTGTGATCTATAAGCTTGGACTTGGAGGCTTGCTGGCCCTCGTTGGTCGGCTCCTCTTGAATGCATTAGCTCGATTCCGTACGGATCGGCTCTCAGTGTGATAATTGTCTACGCTGTGACCGTGAAGTGTTTTGGCGAGCTTCTAACCGTCCATTAGGACAACTTTTTAACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                  
#>  [7] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCCTTTCTTTATTGAAGGGTGAGCTTGGACTTGGGGATTGCTGGTGCAAATTGGCTTCCCTTGAATGTATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACCGTGTAGCTCGGATTGTCTGGCTTCTAATCGTCCTTCACGGGACAATTTGATCTTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                         
#>  [8] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATGTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCCTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#>  [9] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCACCTTTTGGTATTCCGAAAGGTATGCCTGTTTGAGTGTCTTGAATCATCACTTCAAAGAGTTTTCTTTGAATGTGGATTTGGGTGTTGCCATTTGGCTCGCCTTGAATGTGTTAGTGGGTCAGCCCAAGCTAGGTGTAATAAGTTTCATCTGGTGTAAGTGTCGTCCGCTTCCAACTCGCAAGGCTTTGGCCTGCATTTTCATGACTAGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                                                
#> [10] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCATGCCTGTTTGAGTGTCATGGAATTCTCATCCCCTAAATTTATTTTTGGTGGGCTTGGATGTGGAGGTTGTGCCGGTTTGTTAAAAGCCGGCTCCTCTGAAATGCATTAGTGTGAATCTTACTGATCGCCTTCAGTGTGATAATTGTCTACGCTGTGGTTGTGAAGTAAATTAATGTTCATGCTTCTAATCGTCTCTTTTTGAGACAATTACTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                            
#> [11] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGCATTCCGAGGGGCACGCCTGTTCGAGTGTCGTGAACTCCTCCACTCCTTACCTTTTTCGAAGGGTACGGGGCTGGGATTTGAGAGTTGCGGGTCCCTGGCCGATCCGCTCTCTTTGAATGCATTAGCGAAGCCCTCGCGGCCTTGGTGTGATAATTATCTACGCCTCGGCTTAGCGCATGACCGGGCTATGCTTCGAACCGTCTCGCAAGAGACAACGACTACCAACTTGACCTCGAATCAGGCGGGACTACCCGCTGAACTTA"                                                  
#> [12] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCATTAAATTCTCAACCTTTCCAGCTTTTACTAGCTTGGTCAGGCTTGGATCTGGGGGCTGCGGGCTTTTATTCTAAAGTCGGCTCTCCTTAAATGCATTAGCGAAACCTTTGGTTGACCAATCCCTGGTGTGATAATTATCTACGCCACCGATAAACCTTACAAGGCCTTTTTGCTTCTAACCGTCCCTTGGGACAACACTTCATGACATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                     
#> [13] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCATGGTATTCCGTGGAGCATGCCTGTTTGAGTGTCACGTCAACCCTCATTCATGTGATGTTACAGTCGCATGAATGGATTTGGACCGTGTCGTAATCGACTCGTCCCGAAATGCATTAGCTTGCGCTTTTAGAGTGCTGGGCACCGGTGTGATAATTATCTGCGCCAACGCCTTAAGCCTCTTCAGCGGCGCTGCTTACAGCCGTCCTGACGGACAACTATTTTAACGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                    
#> [14] "AAATGCGATAAGTAATGTGAATTGCAGTATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCCCTTTGGTATTCCGAAGGGCATGCCTGTTTGAGTGTCATGAAAATCTCATCCCATCGGGTTTCCGATGACGATGGACTTGGGTGTTGCCGTCCCGGCTCGCCTTAAAAGTGTTAGTGGTTCCGAGAGGGAATCTGGCGTAATAAGTTTCGTCGGTGAATCTCGAAGGCTGCTCATAACCCGTCTTCTGACATCTTTTTGCTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTAA"                                                                                  
#> [15] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAACATATCAACCCCGACGGGATTTTAATCTCGACGGTGGCTTGGACTTGGGAGCCTGCTGGCCGTTGTATTTCGGTCGGCTCTCCTTGAATGCATTAGCTTGGACTTGTTTCATGGGCTGCATAGCGTGATAACTATCTACGCTTTGTACTGTCTGTGGGATACTCATGAGTCTGGCTTCTAATCGTCCTCTGGACAATTACTTTACGTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                  
#> [16] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCATCTTGCGCTCCTTGGTACTCCATGGAGCATGCCTGTTTGAGTGTCACGTAAACCCTCACCCTTGCGATGTAACAGTCGTTCGCGGTGGACTTGGACCGTGCCGTAGTTGGCTCGTCCTGAAATGCATTAGCTGGTGCTTTTAGAGTGCTGGGCGACGGTGTGATAATTATCTGCGCCAATGCCCTGGGCCTCTTCAGCGGTGCCGCTTACAGTCGTCCCCAGAGGACAACGTATCTAAAGCTTTGGCCTCAAATCAGGTAGGATTACCCGCTGAACTTA"                                                  
#> [17] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCGTCGAATCTTTGAACGCACCTTGCGCCCTTTGGTGTTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCGCTTTCTTGTGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCAACTTTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTTGCTTCTGTGGAGACACGCTTTCAAACCGTCCGAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                      
#> [18] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCTTTGGTATTCCGAAGGGCACACCTGTTTGAGTGTCGTGAAATTCTCAACCCTCTTCGCTTTCTTGTGAACGTAGGGATTGGACTTGGAGGCTTTGCCGGGCAACTTTTGCTCGGCTCCTCTCAAATGCATTAGTGCGTCTTGTTGCGACGTGCGCCTCGGTGTGATAATTATCTACGCTGTGGTGCGCTCGCTTCTGTTGGAGACACGCTTTCTAACCGTCCCGTAAAGGACAGCTTTCATCGAACTTTGACCTCAAATCAGGTGGGACTACCCGCTGAACTTA"                                   
#> [19] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTCTCAACCTATAAACTTTTGCGGGTTTGTAGGGTTGGCTATTGGAGGCTTTTGCTGGCCCTCGTTTGAGTCAGCTCCTCTCAAATGCATTAGCTTGGTTCCTTGCGGATCGGCTGTCGGTGTGATAATGTCTACGCCGCGACCGTGAAGTGTTTGGAGAGCTTCTAATGGTCTCGTCAGAGACAGCTTTTATGAACTCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                           
#> [20] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTTTCCGGTATTCCGGAAAGCATGCCTGTTTGAGCGCCTTATCATAGTTGTGCAATTGCATTTTTTGTAGTTGCACAGTATTGGGATTTGCTGCGTGAGCAGCTCACCCTGAATGCATTAGCAATTCGACACTTGAAACTAACGTGATAAGTATCTTCGTTGGCTTGAAGTGAAGACTTGCTTCGAATCGTCTTCGGACAATACCTAATGCTAGGCATCAAATCAGGTAGGATTACCCGCTGAACTTA"                                                                          
#> [21] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCTTTGGTATTCCGAAGAGCATGCCTGTTTGAGTGTCATGAATCTCTCAACTCTTATGATCTGTAATGGATCGTATTGAGCTTGGACTTGGAGGCTTGCCGGCGTTCGTCGGCTCCTCTTAAATGCATTAGCTATGGATCACGATGTGTATGATCTACGGTGTGATATATGTCTACGCCTATGGATCTGCGCGCGAGTTACCGCCTGCAAAGGTAGCGCTCTCCGGAGCCGTGTAGCGAGTGAATACGAAGCTTCTAATCGTCTTCGGACAACTTCTTTGACATCTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"
#> [22] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCCCCTTGGTATTCCGAGGGGCACGCCTGTTCGAGTGTCGTGAAATTCTCAACTCAATCCTCTTGTTATGAGAGGGCTGGGCTTGGACTTGGAGGTCTTGCCGGTTGTGTTCCCTCGGGACGTCGGCTCCTCTTGAATGCATGAGTGGATCCCTTTTTGTAGGGTTTGCCCTTGGTGTGATAATTATCTACGCCGCGGGTAGCCTTGCGCGTTGGTCTGCTTCTAACCGTCCTTCAGGACATTTCATCTCAACTTGACCTCGAATCAGGCGGGACTACCCGCTGAACTTA"                                 
#> [23] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCTATTCTTTATTGAATGGTGAGCTTGGATTTGGAGATTGCTGGTGAAAATCGGCTTCTCTTGAATGAATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACTGTGTAGCTCGGATTGTCTGGCTTCTAATCGTCTTTCATGGACAATTTGATCATTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                          
#> [24] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCTTGGTATTCCGAGGAGCATGCCTGTTTGAGTGTCGTGAAATTATCAACTCCCTTTCTTTATTGAAGGGTGAGCTTGGACTTGGGGATTGCTGGTGCAAATTGGCTTCCCTTGAATGTATTAGCTGGAATTTGATTCGCAGCATATATGCGGTGTGATAATGTCGTCACTGTGTAGCTCGGATTGTCTGGCTTATAATCGTCCTTCACGGGACAATTTGATCTTTTTGACCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                         
#> [25] "AAATGCGATAAGTAATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACCTTGCGCTCCATGGTATTCCGTGGAGCATGCCTGTTTGAGTGTCACGTCAACCCTCATTCATGTGATGTTACAGTCGCATGAATGGATTTGGACCGTGTCGTAATCGACTCGTCCCGAAATGCATTAGCTTGCGCTTTTAGAGTGCTGGGCACCGGTGTGATAATTATCTGCGCCAACGCCTTAAGCCTCTTCAGTGGCGCTGCTTACAGCCGTCCTGACGGACAACTATTTTAACGCTTTGGCCTCAAATCAGGTAGGACTACCCGCTGAACTTA"                                                    
#> 
#> $unique_seqs_2
#> character(0)
#> 
#> $same_name_diff_seq
#> [1] taxa_name   seq_physeq1 seq_physeq2
#> <0 rows> (or 0-length row.names)
#> 
#> $same_seq_diff_name
#> [1] sequence     name_physeq1 name_physeq2
#> <0 rows> (or 0-length row.names)
#> 
#> $mean_nn_kmer_dist_1
#> [1] 0.1210417
#> 
#> $mean_nn_kmer_dist_2
#> [1] NA
#> 
#> attr(,"class")
#> [1] "compare_refseq"

# From a list_phyloseq (uses the first two objects and their names)
lpq <- list_phyloseq(list(full = data_fungi_mini, subset = sub))
#> ℹ Building summary table for 2 phyloseq objects...
#> ℹ Computing comparison characteristics...
#> ℹ Checking sample and taxa overlap...
#> ℹ Detected comparison type: REPRODUCIBILITY
#> ℹ 137 common samples, 20 common taxa
#> ✔ list_phyloseq created (REPRODUCIBILITY)
res3 <- compare_refseq(lpq)
#> Error in compare_refseq(lpq): list_phyloseq must contain at least 2 phyloseq objects.
res3
#> Error: object 'res3' not found