AB Back transcribe - mendessoares/BuddySuite GitHub Wiki

--back_transcribe, -r2d

Description

Convert RNA alignments to DNA.

Examples

Input file: Alignments_rna.physr

 3 100
Mle-Panxα9  ---auguuagacauacuuucaaaguuuaaaggaguuacuccuuuuaaagguauaacgauagaugacgggugggaucaacucaaucggaguuuuauguucg
Mle-Panxα4  augguuauugagcugcuagcuggauacaaaggucuguccccguuuaaagacgcgacuguugacgacucaugggaccaaauaaaccgauguuacguguuca
Mle-Panxα6  auguuauuggagauauuagcgaacuucaaaggagcgacaccuuucaaagaaauaguucuagaugacaagugggaccagauuaaccgauguuacauguucc

 3 100
Ael_PanxβA  ---------------------------------------------------------------------------------------AUGGUAGUCAUUC
Ael_PanxβB  ---------------------------------------------------------------------------------------AUGGUUGUCAUAC
Ael_PanxβC  AUGCCCAACAACAUAUACCCAAACAGACUAUUCGUGAAGACUAAUGAUAUCCCGGAAAAAUUAAACACUCCGUGGUCAUACGAAAAAAUGGUUGUAGUGC

Usage example

$: alb Alignments_rna.physr -r2d

Output

 3 100
Mle-Panxα9  ---atgttagacatactttcaaagtttaaaggagttactccttttaaaggtataacgatagatgacgggtgggatcaactcaatcggagttttatgttcg
Mle-Panxα4  atggttattgagctgctagctggatacaaaggtctgtccccgtttaaagacgcgactgttgacgactcatgggaccaaataaaccgatgttacgtgttca
Mle-Panxα6  atgttattggagatattagcgaacttcaaaggagcgacacctttcaaagaaatagttctagatgacaagtgggaccagattaaccgatgttacatgttcc

 3 100
Ael_PanxβA  ---------------------------------------------------------------------------------------ATGGTAGTCATTC
Ael_PanxβB  ---------------------------------------------------------------------------------------ATGGTTGTCATAC
Ael_PanxβC  ATGCCCAACAACATATACCCAAACAGACTATTCGTGAAGACTAATGATATCCCGGAAAAATTAAACACTCCGTGGTCATACGAAAAAATGGTTGTAGTGC