Skip to content

Commit

Permalink
Merge branch 'jdaw/add-bc-to-summary' into 'master'
Browse files Browse the repository at this point in the history
Add barcode to summary file

Closes DOR-396

See merge request machine-learning/dorado!737
  • Loading branch information
iiSeymour committed Nov 29, 2023
2 parents 1480b00 + 39d667c commit 46bbfdd
Showing 1 changed file with 7 additions and 1 deletion.
8 changes: 7 additions & 1 deletion dorado/cli/summary.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -54,6 +54,7 @@ int summary(int argc, char *argv[]) {
"template_duration",
"sequence_length_template",
"mean_qscore_template",
"barcode",
};

std::vector<std::string> aligned_header = {
Expand Down Expand Up @@ -122,6 +123,11 @@ int summary(int argc, char *argv[]) {
auto num_samples = reader.get_tag<int>("ns");
auto trim_samples = reader.get_tag<int>("ts");

auto barcode = reader.get_tag<std::string>("BC");
if (barcode.empty()) {
barcode = "unclassified";
}

float sample_rate = num_samples / duration;
float template_duration = (num_samples - trim_samples) / sample_rate;
auto exp_start_dt = read_group_exp_start_time.at(rg_value);
Expand All @@ -131,7 +137,7 @@ int summary(int argc, char *argv[]) {
std::cout << filename << separator << read_id << separator << run_id << separator << channel
<< separator << mux << separator << start_time << separator << duration
<< separator << template_start_time << separator << template_duration << separator
<< seqlen << separator << mean_qscore;
<< seqlen << separator << mean_qscore << separator << barcode;

if (reader.is_aligned) {
std::string alignment_genome = "*";
Expand Down

0 comments on commit 46bbfdd

Please sign in to comment.