Advanced Strain Search

ColumbiaRust-1-SW

PSP 10021

DNA Source: SporeWorks

General Information

Accession Date
May 11, 2021
Report Type
Whole-Genome Sequencing
DNA Source Website
sporeworks.com

This visualization shows how genetically distinctive this sample is compared to every other sample in this database. Distinctiveness is the average genetic distance from this sample to its five nearest samples: a sample with no close relatives in the collection scores high (rare), while one with several near-identical samples, for example repeated submissions of the same cultivar, scores low (common). The curve is the distribution of that score across the whole collection and the marker shows where this sample falls, so you can read, as a single percentile, how unique it is within this database rather than just how far it is from its nearest match.

Rarity: Rare Rare

More genetically distinct than 99 of 110 cultivars (91st percentile).

Clonal or not?

Assessing…

Interactive 3D Psilocybe Atlas See ColumbiaRust-1-SW in the tree of life Spin, zoom, and explore exactly where ColumbiaRust-1-SW sits among its closest genetic relatives. Launch 3D tree

Clone Test

Is ColumbiaRust-1-SW a resequencing of a culture already in Psilocydia? Every assembled genome was aligned directly to its closest genome in the collection and the divergence measured. Pairs that are one culture sequenced twice sit at or below about 0.03%; genuinely different cultures sit well above it. Where ColumbiaRust-1-SW falls on that scale answers the question, calibrated against the collection's own resequencings rather than a rule of thumb.

Loading the clone test…

Genetic Information

Loading read coverage…

What is in the assembly

Loading the contigs…

Variants

Changes in the psilocybin genes were found by two independent analysis methods applied to this sample's sequencing reads: aligning the reads to the reference genome, and assembling the genome from them. Each change is labelled with the method that found it, or Both; hover a label for the read counts.

How we found these
Both
In the assembly and in the reads, so confirmed twice.
Reads only
In the reads, not in the assembly. The hover says why where we know: the assembly has the reference base there, does not cover the site, or this sample has no assembly. Many of these are changes carried by only part of the reads; an assembly keeps one base at each position.
Assembly only
In the assembly, but the reads show the reference base there, or too few reads cover the site to say.

Read counts. At every site we count the reads carrying each base (reads with mapping quality of at least 20, bases with quality of at least 13; 8 reads or more to say anything). The change is in most reads at 80% or more, mixed between 20% and 80%, in a minority from 5% to 20% (and at least 2 reads), and essentially absent below that. These counts say what the reads show; they say nothing about how many copies of the gene the sample carries. A change counts as found in the reads when these counts show it (in most, mixed or a minority of reads) or when the lab's variant report lists it. That report lists only changes its gene models rate HIGH or MODERATE, so many changes the reads show are not on it.

Names. Protein changes are named on our curated gene models for the psilocybin cluster. The read-alignment report names them on the public reference's gene models, which in places differ (for example, the public models join psiT1 and psiK into one gene); where the names differ, the hover shows both.

What is listed. Changes that alter the protein (rated HIGH or MODERATE), and every change on the lab's variant report for the cluster. Positions are on chromosome 10 of the reference genome.

psiT1 2 variants — 2 MODERATE
psiT1
MODERATE E121G p.Glu121Gly missense A/G Both
chr10:1663786 54 of 54 reads
psiT1
MODERATE I214L p.Ile214Leu missense A/C Both
chr10:1664184 39 of 39 reads
between psiT1 and psiK 2 variants — 2 MODERATE
between psiT1 and psiK
MODERATE G>A outside the coding sequence G/A Both
chr10:1666040 55 of 55 reads
between psiT1 and psiK
MODERATE T>C outside the coding sequence T/C Both
chr10:1666041 56 of 56 reads
psiK 1 variant — 1 MODERATE
psiK
MODERATE R173G p.Arg173Gly missense A/G Both
chr10:1666622 46 of 46 reads
psiH 2 variants — 2 MODERATE
psiH
MODERATE T502S p.Thr502Ser missense A/T Both
chr10:1669642 54 of 54 reads
psiH
MODERATE L508P p.Leu508Pro missense T/C Both
chr10:1669661 45 of 45 reads
psiM 1 variant — 1 MODERATE
psiM
MODERATE R102K p.Arg102Lys missense C/T Both
chr10:1674573 40 of 40 reads
psiD 1 variant — 1 MODERATE
psiD
MODERATE L22F p.Leu22Phe missense C/A Both
chr10:1677485 50 of 50 reads

Genetic Relatives

Nearest

  1. WFB-UM-WM (PSP10121) — 0.240
  2. WFB-UM (PSP10112) — 0.240
  3. Roatan-Honduras-4-Mush (PSP10079) — 0.314
  4. Arenal-Volcano-1-Mush (PSP10059) — 0.315
  5. Luminous-lucies-2-PS (PSP10032) — 0.316
  6. Luminous-lucies-3-PS (PSP10033) — 0.316
  7. Luminous-lucies-1-PS (PSP10034) — 0.316
  8. Costa-Rica-3-redo-SW (PSP10024) — 0.321
  9. Bplus-ITW-WM (PSP10129) — 0.324
  10. Guapile-ITW-col-WM (PSP10132) — 0.325
  11. Treasure-Coast-2-Mush (PSP10082) — 0.327
  12. Treasure-coast-3-PS (PSP10056) — 0.327
  13. Treasure-coast-2-PS (PSP10055) — 0.327
  14. Treasure-coast-1-PS (PSP10054) — 0.327
  15. Costa-Rica-1-redo-SW (PSP10022) — 0.331
  16. Penvy-Col5-PS (PSP10044) — 0.333
  17. Penvy-1-PS (PSP10045) — 0.333
  18. Penvy-Col4-PS (PSP10043) — 0.333
  19. Penvy-2-SW (PSP10075) — 0.333
  20. Costa-Rica-2-redo-SW (PSP10023) — 0.337

Most Distant

  1. Puerto-Rican-ITW (PSP10078) — 0.611
  2. Blue-Meanie-ITW-WM (PSP10128) — 0.610
  3. GoldenCoast-BAF-B1-UM (PSP10107) — 0.609
  4. Bplus-3-PS (PSP10001) — 0.609
  5. GoldenCoast-BAF-B2-UM (PSP10108) — 0.608
  6. GoldenCoast-BAF-A1-UM (PSP10105) — 0.607
  7. Gulf-Coast-BAF (PSP10090) — 0.607
  8. Huautla-2-redo-SW (PSP10030) — 0.607
  9. GoldenCoast-BAF-A2-UM (PSP10106) — 0.607
  10. Gulf-Coast (PSP10091) — 0.607
  11. Bplus-1-redo-SW (PSP10014) — 0.606
  12. True-Albino-Teacher-ITW-WM (PSP10140) — 0.606
  13. Bplus-2-redo-SW (PSP10015) — 0.606
  14. Bplus-3-redo-SW (PSP10016) — 0.606
  15. Burma-1-SW (PSP10020) — 0.606
  16. Cambodian--1-PS (PSP10060) — 0.605
  17. Bplus (PSP10088) — 0.605
  18. Huautla-3-PS (PSP10005) — 0.605
  19. Orissa-India-ITW (PSP10071) — 0.605
  20. Huautla-3-redo-SW (PSP10031) — 0.605

Sequence Data Downloads

The blockchain-registered VCF above can be downloaded by anyone, so its hash can be checked independently against the blockchain record. The remaining sequence data files — FASTQ, BAM and assemblies — are available to authorised accounts. Log in to download, or contact us for access.

Psilocybin

The psilocybin biosynthesis genes sit together in one cluster on chromosome 10 (psiD, psiK, psiM, psiH, psiR, psiT1, psiT2). Every assembled genome in Psilocydia was compared with the MGC_Penvy_1 reference across this window; below are the variants ColumbiaRust-1-SW carries there, and which other samples carry the same ones. This section compares assemblies only; the psilocybin gene cards under Variants combine the assembly with the reads, so they can list changes that are not here.

Loading the cluster variants…

The enzymes in 3D

This sample's version of each enzyme next to the P. Envy reference's. Flip between the two models; the track underneath shows every change seen in the collection and how many samples carry it, with this sample's changes filled in.

Loading the enzymes…

Blockchain Registration Information

Transaction ID
b4f743fb38a1c8ea62ff86b7de2a895cca847d7313ad5574a8d07ba22bc7438a
Stamping Certificate
Download PDF (39.7 KB)
SHASUM Hash
25cfae81d6087d08790d0d1e23e2970582f2cd4547a03c5307d1f81e9d676c19

Psilocydia uses cookies

By clicking Allow All, you agree to the storing of cookies on your device to enhance site navigation, analyze site usage, and assist in our marketing efforts.

Customize Settings