Barcode Statistics Documentation
isoseq bcstats
emits statistics for each barcode:
- Barcode sequence
- Number of reads matching the barcode
- Frequency Rank (within barcodes)
- Number of unique molecular barcodes matching this barcode
- Whether the barcode is Group/Cell barcode or a Molecular Barcode/UMI
This tool should be run after isoseq correct
.
# Example:
isoseq bcstats --json sample.bcstats.json -o sample.bcstats.tsv corrected.bam
If --json
is unset, JSON summary information is written to stderr (“/dev/stderr”). Similarly, if ‘-o’ is unset, output TSV information is written to stdout (“/dev/stdout”).
By default, the program only emits stats on group barcodes. Adding --umi
will cause stats for the full molecular barcodes to be emitted as well.
Note that, if isoseq correct
was called using the percentile method (see Cell Calling), then isoseq bcstats
need to use the same parameters as well. For example:
isoseq correct --method percentile --percentile 90 ...
isoseq bcstats --method percentile --percentile 90