我有一个 bam 文件并使用 bioperl (Bio::DB::Sam) 来处理它。现在我想问一下是否有可能在这个文件的对齐中添加标签?
我用
my $iterator = $bam->features(-iterator => 1,
-flags => {M_UNMAPPED=>0});
while (my $align = $iterator->next_seq) {
...
}
循环通过对齐的读取。现在我正在搜索类似的东西
$align->addTag(key=>value)
再见