Costa-Rica-3-redo-SW

PSP 10024

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: Ultra-rare Ultra-rare

More genetically distinct than 109 of 110 cultivars (100th percentile).

Clonal or not?

Assessing…

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

Clone Test

Is Costa-Rica-3-redo-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 Costa-Rica-3-redo-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 7 variants — 2 HIGH   5 MODERATE
psiT1
MODERATE E121G p.Glu121Gly missense A/G Reads only
chr10:1663786 10 of 10 reads
psiT1
MODERATE I214L p.Ile214Leu missense A/C Reads only
chr10:1664184 9 of 9 reads
psiT1
HIGH S223* p.Ser223Ter stop gained C/G Reads only
chr10:1664264 3 of 6 reads
psiT1
MODERATE T246I p.Thr246Ile missense C/T Reads only
chr10:1664333 4 of 4 reads
psiT1
MODERATE R266S p.Arg266Ser missense A/T Reads only
chr10:1664394 6 of 6 reads
psiT1
MODERATE T386A p.Thr386Ala missense A/G Reads only
chr10:1664803 11 of 11 reads
psiT1
HIGH L461* p.Leu461Ter stop gained T/A Reads only
chr10:1665085 2 of 15 reads
between psiT1 and psiK 9 variants — 1 HIGH   8 MODERATE
between psiT1 and psiK
MODERATE A>G outside the coding sequence A/G Reads only
chr10:1665277 8 of 8 reads
between psiT1 and psiK
MODERATE A>G outside the coding sequence A/G Reads only
chr10:1665560 9 of 9 reads
between psiT1 and psiK
MODERATE G>A outside the coding sequence G/A Reads only
chr10:1665570 10 of 10 reads
between psiT1 and psiK
MODERATE A>C outside the coding sequence A/C Reads only
chr10:1665631 8 of 8 reads
between psiT1 and psiK
MODERATE G>A outside the coding sequence G/A Reads only
chr10:1665634 8 of 8 reads
between psiT1 and psiK
MODERATE A>G outside the coding sequence A/G Reads only
chr10:1665658 7 of 7 reads
between psiT1 and psiK
HIGH G>T outside the coding sequence G/T Reads only
chr10:1665665 10 of 10 reads
between psiT1 and psiK
MODERATE G>A outside the coding sequence G/A Reads only
chr10:1666040 8 of 8 reads
between psiT1 and psiK
MODERATE T>C outside the coding sequence T/C Reads only
chr10:1666041 9 of 9 reads
psiK 5 variants — 5 MODERATE
psiK
MODERATE S101A p.Ser101Ala missense T/G Reads only
chr10:1666406 12 of 12 reads
psiK
MODERATE L110V p.Leu110Val missense T/G Reads only
chr10:1666433 8 of 8 reads
psiK
MODERATE E167G p.Glu167Gly missense A/G Reads only
chr10:1666605 18 of 18 reads
psiK
MODERATE R173G p.Arg173Gly missense A/G Reads only
chr10:1666622 20 of 20 reads
psiK
MODERATE Q279L p.Gln279Leu missense A/T Reads only
chr10:1666941 6 of 6 reads
psiH 1 variant — 1 MODERATE
psiH
MODERATE K236N p.Lys236Asn missense G/T Reads only
chr10:1668429 4 of 4 reads
psiT2 3 variants — 3 MODERATE
psiT2
MODERATE E291Q p.Glu291Gln missense C/G Assembly only
chr10:1671235 6 of 6 reads
psiT2
MODERATE S139T p.Ser139Thr missense C/G Both
chr10:1671801 11 of 11 reads
psiT2
MODERATE G63D p.Gly63Asp missense C/T Both
chr10:1672192 9 of 9 reads
psiM 2 variants — 1 HIGH   1 MODERATE
psiM
MODERATE R102K p.Arg102Lys missense C/T Reads only
chr10:1674573 14 of 14 reads
psiM
HIGH T>A splice acceptor T/A Reads only
chr10:1674655 4 of 10 reads
psiD 1 variant — 1 MODERATE
psiD
MODERATE L22F p.Leu22Phe missense C/A Reads only
chr10:1677485 7 of 7 reads
Fvcpc2 1 variant — 1 MODERATE
Fvcpc2
MODERATE p.Ser201Pro c.601T>C missense variant T/C
CM038999.1:1973492

Genetic Relatives

Nearest

  1. Costa-Rica-2-redo-SW (PSP10023) — 0.128
  2. Costa-Rica-1-redo-SW (PSP10022) — 0.134
  3. Costa-Rica-1-redo-SW (PSP10003) — 0.202
  4. Costa-Rica-2-redo-SW (PSP10004) — 0.210
  5. Penvy-2-SW (PSP10075) — 0.296
  6. Penvy-Col5-PS (PSP10044) — 0.296
  7. Penvy-1-PS (PSP10045) — 0.296
  8. Penvy-Col4-PS (PSP10043) — 0.296
  9. Roatan-Honduras-4-Mush (PSP10079) — 0.310
  10. Arenal-Volcano-1-Mush (PSP10059) — 0.314
  11. Luminous-lucies-2-PS (PSP10032) — 0.316
  12. Luminous-lucies-3-PS (PSP10033) — 0.316
  13. Luminous-lucies-1-PS (PSP10034) — 0.316
  14. Treasure-coast-1-PS (PSP10054) — 0.317
  15. Treasure-coast-3-PS (PSP10056) — 0.317
  16. Treasure-Coast-2-Mush (PSP10082) — 0.317
  17. Treasure-coast-2-PS (PSP10055) — 0.317
  18. ColumbiaRust-1-SW (PSP10021) — 0.321
  19. Bplus-ITW-WM (PSP10129) — 0.328
  20. Orissa-India-4-Mush (PSP10073) — 0.331

Most Distant

  1. Bplus-3-PS (PSP10001) — 0.639
  2. GoldenCoast-BAF-A1-UM (PSP10105) — 0.637
  3. GoldenCoast-BAF-B1-UM (PSP10107) — 0.636
  4. GoldenCoast-BAF-B2-UM (PSP10108) — 0.635
  5. GoldenCoast-BAF-A2-UM (PSP10106) — 0.635
  6. Puerto-Rican-ITW (PSP10078) — 0.635
  7. Orissa-India-ITW (PSP10071) — 0.634
  8. Bplus-2-redo-SW (PSP10015) — 0.633
  9. Huautla-3-PS (PSP10005) — 0.632
  10. Bplus-1-redo-SW (PSP10014) — 0.632
  11. Bplus-3-redo-SW (PSP10016) — 0.632
  12. Maui-2-PS (PSP10041) — 0.631
  13. Huautla-3-redo-SW (PSP10031) — 0.631
  14. Huautla-2-redo-SW (PSP10030) — 0.631
  15. Ecuador-3-SW (PSP10025) — 0.631
  16. GoldenTeacher-3-SW (PSP10029) — 0.631
  17. McKennaii-1-PS (PSP10069) — 0.630
  18. Maui-3-PS (PSP10042) — 0.630
  19. True-Albino-Teacher-ITW-WM (PSP10140) — 0.630
  20. Maui-1-PS (PSP10040) — 0.630

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 Costa-Rica-3-redo-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…

Blockchain Registration Information

Transaction ID
8c50e820eef1c5e91e17155e896dc6f9265bce833f924b443d238e95fad3096f
Stamping Certificate
Download PDF (39.6 KB)
SHASUM Hash
b6d4497f3c8e7e714c77dd267d999a13dfb7999a15e42678dc5966a05392361c

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