AB Reverse transcribe

--reverse_transcribe, -r2d

Description

Convert RNA to cDNA.

Example

Input file: Alignments_cds.phy

 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

$: sb Alignments_cds.phy -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