bold coordinates question

Hi all,

Is there a way to get longitude and latitude information related to each BOLD accesion? the function bold_seq gives back four entries; id, name, gene, and sequence but not lon and lat.

Thanks in advance!

2 Likes

Thanks for your question about the bold package @nereal

Checkout the function bold_seqspec(), which gets sequence data and specimen data, e.g.,

res <- bold_seqspec(taxon = "Coelioxys")
tibble::as_tibble(res[, c("lat", "lon", "nucleotides")])
#> # A tibble: 315 x 3
#>      lat      lon nucleotides
#>    <dbl>    <dbl> <chr>
#>  1  47.4   16.3   ---------------------------------------------------------------------------------------------------------------------------------------------AATAATGATATAATTTATAACā€¦
#>  2  52.9 -118.    TATAATATATATAATTTTTGCAATATGATCAGGAATAATTGGATCTTCTATAAGTATAATTATTCGAATAGAATTAAGAATCCCAGGATCATGAATTAATAATGATCAAATTTATAACTCTTTTATTACAGCACATGCATTTTTAATAATTTTTTTTTTAGTā€¦
#>  3  NA     NA     -------------------------------GGTATAATTGGATCATCTTTAAGAATAATTATTCGCATAGAATTAAGAATCCCGGGTTCTTGAATTAACAATGATCAAATTTATAATTCTTTTATTACAGCTCATGCCTTTTTAATAATTTTTTTCCTAGTā€¦
#>  4  NA     NA     CATCCTATATATAATTTTTGCCATATGGTCAGGAATAATTGGATCTTCATTAAGAATAATTATTCGTATAGAATTAAGAATCCCAGGCTCTTGGATTAGTAATGACCAAATTTATAATTCTTTTATTACTGCTCATGCATTTTTAATAATTTTTTTTTTAGTā€¦
#>  5  51.5   -0.967 ATAAAGATATTGGTATAATATATATAATTTTTGCTATATGATCAGGAATAATTGGATCTTCATTAAGAATAATTATTCGAATAGAATTAAGAATTCCAGGGTCTTGAATTAATAATGACCAAATTTATAATTCTTTCATTACAGCCCATGCATTTTTAATAAā€¦
#>  6  NA     NA     ----------------------------------ATAATTGGATCTTCTTTAAGAATAATCATTCGAATAGAATTAAGAACTCCAGGATCTTGAATTAATAATGATCAAATCTATAACTCATTTATTACCGCCCATGCATTTCTAATAATTTTTTTTTTAGTā€¦
#>  7  NA     NA     AATATTATATATAATTTTTGCTATATGATCTGGAATAATTGGATCTTCATTAAGAATAATTATTCGAATAGAATTAAGAATTCCAGGGTCATGAATTAATAATGATCAAATTTATAATTCTTTCATTACAGCTCATGCATTTTTAATAATTTTTTTTTTAGTā€¦
#>  8  18.8  -96.9   -ATATTATATATAATTTTTGCTATATGATCTGGAATAATTGGATCTTCATTAAGAATAATTATTCGAATAGAATTAAGAATTCCAGGATCATGAATCAATAATGATCAAATTTATAATTCTTTTATTACAGCTCATGCATTTTTAATAATTTTTTTTTTAGTā€¦
#>  9  NA     NA     ---------------------------------AATAATTGGATCATCAATTAGAATATTAATTCGAATAGAATTAAGAATTCCTGGAACATGAATTAATAATGATCAAATCTATAATTCATTTATTACCGCTCATGCATTTTTAATAATTTTTTTTTTAGTā€¦
#> 10  52.2   13.5   AATAATATATATAATTTTTGCTATATGATCAGGAATAATTGGTTCTTCATTAAGAATAATTATTCGTATAGAATTAAGAATTCCTGGAGCTTGAATTAATAATGATCAAATTTATAATTCATTTATTACAGCTCACGCATTTATAATAATTTTTTTTTTAGTā€¦
#> # ā€¦ with 305 more rows

Does that get what you want?

1 Like

Yes, fantastic! sorry i missed that funtion out

Thank you Sckott

1 Like