Output Files¶

This guide documents all output files produced by the pipeline.

Output Directory Structure¶

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12

{output_directory}/
├── pod5/                    # Merged POD5 files
├── bam/                     # BAM files at each stage
├── fq/                      # Extracted FASTQ files
├── summary/                 # Analysis outputs
│   ├── tables/             # Tabular summaries
│   ├── modkit/             # Modification calling
│   └── qc/                 # Reference QC metrics
├── reports/                 # Rendered QC reports
├── demux/                   # Demultiplexing outputs (if enabled)
├── logs/                    # Rule execution logs
└── squiggy-session.json     # Squiggy session file for Positron

Data Flow and Outputs¶

flowchart TB
    subgraph Input
        A[POD5 files]
    end

    subgraph Processing
        B[pod5/{sample}/{sample}.pod5<br/>Merged POD5]
        C[bam/rebasecall/{sample}/{sample}.rbc.bam<br/>Basecalled]
        D[fq/{sample}/{sample}.fq.gz<br/>FASTQ]
        E[bam/aln/{sample}/{sample}.aln.bam<br/>Aligned]
        F[bam/charging/{sample}/{sample}.charging.bam<br/>Classified]
        G[bam/final/{sample}/{sample}.bam<br/>Final BAM]
    end

    subgraph Outputs
        H[summary/tables/<br/>Charging & Stats]
        I[summary/modkit/<br/>Modifications]
        J[summary/qc/<br/>Reference similarity]
        K[summary/tables/<br/>Odds ratios]
        L[reports/<br/>QC report]
    end

    A --> B --> C --> D --> E --> F --> G
    G --> H
    G --> I
    G --> J
    I --> K
    H --> K
    H --> L

Core Outputs¶

Final BAM¶

bam/final/{sample}/{sample}.bam

The final BAM file with charging classification and adapter position tags.

Tags:

View tags:

1

samtools view results/bam/final/sample1/sample1.bam | head -1 | tr '\t' '\n' | grep -E "^(CL|CM|PT):"

Charging Probability Table¶

summary/tables/{sample}/{sample}.charging_prob.tsv.gz

Per-read charging likelihood scores.

Interpretation:

• Score ≥ 200: Charged (aminoacylated)
• Score < 200: Uncharged

Example:

1

zcat results/summary/tables/sample1/sample1.charging_prob.tsv.gz | head

1
2
3

read_id                                 tRNA                    charging_likelihood
00a1b2c3-4567-89ab-cdef-0123456789ab   tRNA-Ala-AGC-1-1        245
00a1b2c3-4567-89ab-cdef-0123456789ac   tRNA-Gly-GCC-2-1        87

Charging CPM Table¶

summary/tables/{sample}/{sample}.charging.cpm.tsv.gz

Per-tRNA aggregated charging counts, normalized to CPM (counts per million).

Example:

1

zcat results/summary/tables/sample1/sample1.charging.cpm.tsv.gz | column -t | head

1
2
3

tRNA              counts_charged  counts_uncharged  cpm_charged  cpm_uncharged
tRNA-Ala-AGC-1-1  1523            234               15230.5      2340.2
tRNA-Gly-GCC-2-1  892             1456              8920.3       14560.8

Quality Control Outputs¶

Alignment Statistics¶

summary/tables/{sample}/{sample}.align_stats.tsv.gz

Read counts through pipeline stages.

Base Calling Errors¶

summary/tables/{sample}/{sample}.bcerror.tsv.gz

Per-position base calling error metrics.

Coverage Tracks¶

summary/tables/{sample}/{sample}.{cpm,counts}.bg.gz

BedGraph coverage tracks for visualization.

• .cpm.bg.gz - CPM-normalized coverage
• .counts.bg.gz - Raw count coverage

Load in IGV:

1
2

gunzip -c results/summary/tables/sample1/sample1.cpm.bg.gz > sample1.cpm.bg
# Load sample1.cpm.bg in IGV

Modification Outputs¶

Modification Pileup¶

summary/modkit/{sample}/{sample}.pileup.bed.gz

Per-site modification consensus in BED format.

Per-Read Modification Calls¶

summary/modkit/{sample}/{sample}.mod_calls.tsv.gz

Individual modification calls per read.

Full Modification Export¶

summary/modkit/{sample}/{sample}.mod_full.tsv.gz

Comprehensive modification information including all modkit fields.

Reference Similarity Matrix¶

summary/qc/reference_similarity.tsv

Pairwise sequence similarity matrix for the reference FASTA, useful for identifying potential cross-mapping issues.

1

pixi run snakemake compute_reference_similarity --configfile=config/config.yml

Format: Square TSV matrix with sequence names as row and column headers, values are percent identity (0-100).

Modification Odds Ratios¶

summary/tables/{sample}/{sample}.odds_ratios.tsv.gz

Per-tRNA pairwise modification odds ratios testing whether modification at one position is correlated with modification at another position (or with charging status).

1

pixi run snakemake compute_odds_ratios --configfile=config/config.yml

QC Report¶

reports/qc_report.html

A combined Quarto HTML report with per-sample QC tabs, including alignment statistics, charging distributions, and basecalling error metrics.

1

pixi run -e report snakemake render_combined_qc_report --configfile=config/config.yml

Squiggy Session File¶

squiggy-session.json

A JSON session file generated at the root of the output directory for loading pipeline outputs in the Squiggy extension for Positron IDE.

Contents:

• Relative paths to POD5, BAM, and reference FASTA files for each sample
• MD5 checksums and file metadata for integrity verification
• Default plot options (eventalign mode, z-normalization)

JSON structure:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14

{
  "version": "1.0.0",
  "timestamp": "...",
  "sessionName": "aa-tRNA-seq: ...",
  "samples": {
    "sample1": {
      "pod5Paths": ["pod5/sample1/sample1.pod5"],
      "bamPath": "bam/final/sample1/sample1.bam",
      "fastaPath": "../path/to/reference.fa"
    }
  },
  "plotOptions": { ... },
  "fileChecksums": { ... }
}

Usage:

Open the squiggy-session.json file in Positron to load all samples with their associated POD5, BAM, and reference files.

Intermediate Files¶

These files are produced but typically not used directly:

Merged POD5¶

pod5/{sample}/{sample}.pod5

Merged POD5 file containing all raw signal data for the sample.

Rebasecalled BAM¶

bam/rebasecall/{sample}/{sample}.rbc.bam

Dorado output with basecalls and move tables.

Aligned BAM¶

bam/aln/{sample}/{sample}.aln.bam

BWA MEM alignment output.

Charging BAM¶

bam/charging/{sample}/{sample}.charging.bam

Remora classification output with ML/MM tags (before renaming).

FASTQ¶

fq/{sample}/{sample}.fq.gz

Extracted reads for alignment.

Demultiplexing Outputs¶

When demultiplexing is enabled:

Barcode Mapping¶

demux/read_ids/{run_id}/barcode_mapping.tsv.gz

Read ID to barcode assignments.

Per-Sample Read Lists¶

demux/read_ids/{sample}.txt

Read IDs belonging to each sample.

Split POD5¶

demux/pod5/{sample}.pod5

Per-sample POD5 files after demultiplexing.

Signal Metrics (Optional)¶

If remora_kmer_table is configured:

summary/tables/{sample}/{sample}.remora.tsv.gz

Remora signal metrics per read per position.

Log Files¶

logs/{rule}/{sample}.log

Standard output and error for each rule execution.

File Sizes¶

Approximate file sizes for a typical sample:

Cleanup¶

Remove intermediate files to save space:

1
2
3
4
5

# Remove intermediate BAMs (keep final)
rm -rf results/bam/rebasecall results/bam/aln results/bam/charging

# Remove FASTQ (can be regenerated)
rm -rf results/fq

Next Steps¶

• Workflow Overview - Understand the pipeline stages
• Rules Reference - Detailed rule documentation