Bio::T...gsToSeq(3pm)

πŸ“¦
version 1.4.5 Β· 3 versions
, in debian bookworm, ubuntu jammy

🌸 Bio::Tradis::AddTagsToSeq πŸ’…

okay bestie, so this whole page is about Bio::Tradis::AddTagsToSeq (the manpage header calls it Bio::T...gsToSeq(3pm) because manpage names get truncated, it's giving mysterious but it's literally just this module). It takes a BAM file and makes a brand new BAM with tr and tq tags added onto the sequence and quality strings. Iconic behaviour, no cap.


✨ VERSION ✨

version 1.4.5 β€” she's a specific girl, respect her.


πŸ¦‹ SYNOPSIS πŸ¦‹

So Bio::Tradis::AddTagsToSeq parses BAM files, adds the tags you give it to the start of the sequence, and creates a temporary SAM file, which then gets converted to BAM. That's the whole pipeline, bestie.

use Bio::Tradis::AddTagsToSeq;

my $pipeline = Bio::Tradis::AddTagsToSeq->new(bamfile => 'abc');
$pipeline->add_tags_to_seq();

πŸ’– NAME πŸ’–

Bio::Tradis::AddTagsToSeq β€” that's her, that's the name, periodt.


πŸŽ€ PARAMETERS πŸŽ€

Required (you HAVE to bring these, fr fr)

  • bamfile β€” path to / name of the file containing the reads and the tags. Not optional, don't try it.

Optional (nice to have 🫢)

  • outfile β€” the name you want to give the output BAM. If you don't say anything, it defaults to file.tr.bam for an input file named file.bam. She'll name it for you, she's thoughtful like that.

πŸ”₯ METHODS πŸ”₯

  • add_tags_to_seq β€” adds the TraDIS tags to the reads.

Here's the tea though, because the two cases are not the same:

  • For unmapped reads, the tag is added to the start of the read sequence and quality strings.
  • For reads where the flag indicates that it is mapped and reverse complemented, the reverse complemented tags are added to the end of the read strings.

Why, you ask? Because many conversion tools (e.g. picard) take the read orientation into account and will re-reverse the mapped/rev comp reads during conversion β€” which leaves all the tags in the correct orientation at the start of the sequences in the resulting FastQ file. It's giving galaxy brain. πŸ’…


🌈 AUTHOR 🌈

Carla Cummins <path-help@sanger.ac.uk> β€” shoutout to the queen who wrote this. πŸ‘‘


This software is Copyright (c) 2013 by Wellcome Trust Sanger Institute.

This is free software, licensed under:

The GNU General Public License, Version 3, June 2007

and that's the whole page, bestie πŸ’œ this was perl v5.30.0, dated 2020-01-11, and it's filed under Bio::T...gsToSeq(3pm). go add some tags, you got this βœ¨πŸ¦‹πŸŽ€

πŸ“¦ where this page comes from