MATPredict · mating-type report

AvariabilisNCCPF102052

organism not recorded

Run by MATPredict 0.8.1 (v0.8.1, 3a62e6302b) · database 86f18d232e4f · local NCBI taxonomy snapshot 2026-10-01 · rendered 2026-11-07 02:40 UTC by MATPredict 0.8.1

Result

Mating type

Called loci 1

Minus MAT family · Mucoromycota

high confidence
Gene order of the Minus locus on MZZL01000138Minus locus: genes from left to right on MZZL01000138, 192,120 to 199,388: sexM (core MAT gene, minus strand); sexP (core MAT gene, minus strand); tptA (conserved flank, plus strand); algA (variable flank, plus strand).192 kb194 kb196 kb198 kb← 191.8 kb more to contig endcalled locussexM: core MAT gene, 192,375–192,626 (-), identity 65.48%, the two modelling tools disagreesexMsexP: core MAT gene, 192,375–192,659 (-), identity 30.53%, gene modelsexP (lost to sexM)tptA: conserved flank, 195,578–196,766 (+), identity 73.64%, the two modelling tools disagreetptAalgA: variable flank, 197,921–199,388 (+), identity 49.73%, the two modelling tools disagreealgA
Gene order on MZZL01000138. Scroll sideways to see the whole figure.
  • core MAT gene
  • conserved flank
  • variable flank
  • gene model
  • tools disagree
  • alternate model
  • sexP is a weak hit to the other idiomorph's gene and was not counted against the call.
Location
MZZL01000138:192,120–199,388 · 7.3 kb
Own genes
MZZL01000138:192,375–199,388 · 7.0 kb · the span runs 255 bp beyond these genes
Locus
Complete locus: MAT gene(s) with their flanking genes.
Genes modelled
3 of 4 gene hits
Nearest contig end
192.1 kb
Region (samtools)
MZZL01000138:192120-199388
Idiomorph scores
  • Minus 173.8
  • Plus 83.2

Margin over the next idiomorph: 90.5

Score: HMM bitscore of the MAT protein against each idiomorph's profile.

Profile-HMM classifier
  • Minus 173.8
  • Plus 83.2

Margin 90.5 (a call needs 25.0 or more). Scored: sexM, sexM (model).

Overlapping hits resolved

Where genes of two idiomorphs hit the same place, the stronger hit was kept.

keptidentity %set asideidentity %overlap %basis
sexM65.5sexP30.5100.0profile bitscore of its own idiomorph

Gene evidence

GeneRoleModel Identity %E-value Position (strand)Exons Reference record
sexMcore MATmodels differexonerateminiprot: 192,120–192,719, 39.6%65.5cov 39.8%3.8e-38140 bits192,375–192,626 -14837_nrrl1555_MAT_Minus
sexPcore MATother idiomorph, not countedlost to sexM, not countedgene modelexonerate30.5cov 32.0%1.8e+0031 bits192,375–192,659 -14837_ubc21_MAT_Plus
tptAconserved flankmodels differexonerateminiprot: 195,419–196,799, 67.2%73.6cov 79.8%8.5e-157361 bits195,578–196,766 +23108396_ai-2023a_MAT_Plus
algAvariable flankmodels differexonerateminiprot: 197,921–199,388, 48.4%49.7cov 98.4%8.1e-93182 bits197,921–199,388 +529922_cbs109-16_MAT_Minus
rnhAexpected, not found
glrAexpected, not found

Withheld candidate loci 6

Places with some MAT-gene evidence that did not pass a quality gate. Most genomes have several. They are listed so a known locus can be checked against them, not because they are likely calls.

Reasons: fewer than 2 genes could be modelled (4); core protein closer to a known non-MAT paralog (1); too few of the family's genes found (1).

Show the 6 withheld loci
familylocationidiomorph genesmodelledreason
MATMZZL01000137:511,352–532,286undeterminedsexP, sexM, rnhA0fewer than 2 genes could be modelled
MATMZZL01000409:162,681–205,772undeterminedtptA, sexM, rnhA0fewer than 2 genes could be modelled
MATMZZL01000389:353,473–378,681undeterminedsexP, rnhA, algA0fewer than 2 genes could be modelled
MATMZZL01000014:1,879–9,841undeterminedsexP, sexM, rnhA0fewer than 2 genes could be modelled
MATMZZL01000384:34,749–118,488undeterminedtptA, sexP, sexM, glrA2core protein closer to a known non-MAT paralog
MATMZZL01000384:386,367–432,643undeterminedsexM, btbA1too few of the family's genes found

Provenance

Genome file
AvariabilisNCCPF102052.genome.fasta
Assembly
411 contigs · 39.3 Mb · N50 647.5 kb
Genome SHA-256
bf6153c3565c131380bb26691a0ec901f23bbd0a40465f2f16c0251125eccf6f
MATPredict
0.8.1 (v0.8.1, 3a62e6302b)
Database
86f18d232e4fe95d · 185 records
Taxonomy
local NCBI taxonomy snapshot 2026-10-01
Parameters
phylum=Mucoromycota, phylum_guess={'requested': 'auto', 'status': 'ok', 'usable': True, 'reason': None, 'phylum': 'Mucoromycota', 'phylum_taxid': 1913637, 'deepest_rank': 'species', 'routing_taxid': 34489, 'routing_taxid_rank': 'family', 'flags': ['identical_to_reference'], 'gmg_version': 'gmg 0.3.2 (66055154a)', 'index': 'bfd-reps-2026.10.5-s20000.gmgi', 'index_sha256': '0f319567f1cab17643e1daf7ffdc0e3b312bd43913e351568108d67a4a19702c', 'index_spec': 'bfd-reps-2026.10.5-s20000.gmgi', 'note': 'GmG called phylum Mucoromycota', 'action': 'routed', 'weak_call_accepted': False, 'genetic_code': 1, 'genetic_code_from': 'family taxid 34489', 'routed_by': 'family taxid 34489 (lineage)'}, genetic_code=1, min_hits=2, require_core_role=True, max_polished_clusters_per_family=6, polish_strong_identity=50.0, polish_core_pair=True, partial_core_policy=off, homothallic_policy=off, confidence_cap_policy=off, cross_family_policy=off, fallback_after_withheld_policy=off, strong_single_gene_policy=off, scoped_roster_policy=off, supported_min_bitscore=33.0
Run started
2026-11-07T02:39:47Z
Run time
51 s
Files
detection_report.yaml, detected_loci.gff3

Cite MATPredict as given in CITATION.cff at github.com/stajichlab/MATPredict. The reference records named in the gene tables are curated database entries; each record's metadata cites its source publication.

How to read this report

Idiomorph
The form of the mating-type locus a strain carries (for example MAT1-1 or MAT1-2, Plus or Minus). Basidiomycota HD and PR alleles are named after the reference they match.
High confidence
Every expected core MAT gene was found and modelled, with a conserved flanking gene where the family has one, on one contig.
Medium confidence
The core genes were found, but a gene could not be modelled, no conserved flank was found, or the locus is split (one step down).
Low confidence
A single isolated gene hit, or no core MAT gene.
Model status
Gene model: the reference protein was aligned to the genome to build the gene (exonerate, miniprot). Models differ: the two tools disagree on the structure; the alternate model is shown. Hit only: a search hit with no model, weaker evidence.
Identity
Percent amino-acid identity between the gene (or hit) and the closest reference protein.
Reference record
The curated database record whose protein matched best. Each record's metadata cites the publication or deposit it comes from.

Limits