SB Reverse complement - mendessoares/BuddySuite GitHub Wiki
--reverse_complement, -rc
Description
Reverse compliment sequences and features.
Example
Input file: Mle-Panxα2.gb
LOCUS Mle-Panxα2 1314 bp DNA UNA 02-JAN-2015
DEFINITION cDNA - ML25998a.
ACCESSION Mle-Panxα2
VERSION Mle-Panxα2
KEYWORDS .
SOURCE
ORGANISM . . .
.
FEATURES Location/Qualifiers
CDS order(1..144,145..307,308..555,556..688,689..810,811..1314)
TMD1 85..147
TMD2 421..483
TMD3 682..744
TMD4 931..1023
ORIGIN
1 atggtattgg atctcatttc tggaagcttg aatggctttt taaagatcaa gtcagttagc
61 atcgacgatc agtgggacca gattaacaga acctatttgg tcatgttttg tattttatct
121 ggtacaatca tgacctttaa acagaattta ggatcaataa tacactgtat atcggatgca
181 agaggcgacg acagttcgtt tgcggatgct catgcgacat ttgtgcaaga ctattgtgct
241 gctcaagggc tgtacacttt aaaagaagtg tatgacaagt cttggccaga tgaaattcct
301 tacccaggta ttctccaaat gaaaacaatc ggttgtttcc cggggagaca gttcaaaaac
361 ggaaccccca tccagtgccc ggacgagaaa gatctgaaac ccttcacaac ggtctatcat
421 gtctggtaca tgttcgtacc gttctacttc tgcgctgttg gcatcgcttt ttacttcccc
481 tacacggttt tcagacacct cagcggcatc tacgacatca agcctatgtt gaacagcctt
541 gccctcgaca ttggggccta cacggaggag gacataagtc gacgtataga caatgtctcg
601 aggtggttgt acatcaagtt ggatccctac atgaacaaca tgcttcctta tactcagata
661 gttcacaaac attccatctt ttacacggtg atgttggtga aggtgatgta cctagctacc
721 agtgtttcta ttttttacgc cactcaccgg atattcgacc aaggaaactt tgcactctac
781 ggatacgatg ttctaatgag cataccacag gaaacaagct ataaagtgat ggacacaatc
841 ttccctaaaa tggttggctg tgagatcaac atgtggggcc ggactggcga acagagcgaa
901 tctcttctgt gtgtcctccc tcaaaacatc ggcaaccaat acttcttcct tatattctgg
961 tttctcctga ttctcaccat actttccaac tgtatctctg taatagtgac catattcaga
1021 tttatattcg ttagtgggag ctacaaaagg ttcctggcta ccagcctctt gaatcacgaa
1081 gaacgataca agctggtgtt tacacatgtc ggcacgactg gaagatacat tttactgctc
1141 tgtgccgatc atagcaaccc caaaatattc gaggatcttc tagagatcgt ctgttccctt
1201 ctcatagcaa actatcacaa aagaaagagg agtcgggata agggacacag tcgagcggag
1261 ggggtaggga ctaaagggcg acacggtctg tcttttgtgg actcaaccgt gtga
//
Usage
$: sb Mle-Panxα2.gb -rc
Output
LOCUS Mle-Panxα2 1314 bp DNA UNA 02-JAN-2015
DEFINITION cDNA - ML25998a.
ACCESSION Mle-Panxα2
VERSION Mle-Panxα2
KEYWORDS .
SOURCE
ORGANISM . . . .
.
FEATURES Location/Qualifiers
CDS complement(order(1..504,505..626,627..759,760..1007,
1008..1170,1171..1314))
TMD1 complement(1168..1230)
TMD2 complement(832..894)
TMD3 complement(571..633)
TMD4 complement(292..384)
ORIGIN
1 tcacacggtt gagtccacaa aagacagacc gtgtcgccct ttagtcccta ccccctccgc
61 tcgactgtgt cccttatccc gactcctctt tcttttgtga tagtttgcta tgagaaggga
121 acagacgatc tctagaagat cctcgaatat tttggggttg ctatgatcgg cacagagcag
181 taaaatgtat cttccagtcg tgccgacatg tgtaaacacc agcttgtatc gttcttcgtg
241 attcaagagg ctggtagcca ggaacctttt gtagctccca ctaacgaata taaatctgaa
301 tatggtcact attacagaga tacagttgga aagtatggtg agaatcagga gaaaccagaa
361 tataaggaag aagtattggt tgccgatgtt ttgagggagg acacacagaa gagattcgct
421 ctgttcgcca gtccggcccc acatgttgat ctcacagcca accattttag ggaagattgt
481 gtccatcact ttatagcttg tttcctgtgg tatgctcatt agaacatcgt atccgtagag
541 tgcaaagttt ccttggtcga atatccggtg agtggcgtaa aaaatagaaa cactggtagc
601 taggtacatc accttcacca acatcaccgt gtaaaagatg gaatgtttgt gaactatctg
661 agtataagga agcatgttgt tcatgtaggg atccaacttg atgtacaacc acctcgagac
721 attgtctata cgtcgactta tgtcctcctc cgtgtaggcc ccaatgtcga gggcaaggct
781 gttcaacata ggcttgatgt cgtagatgcc gctgaggtgt ctgaaaaccg tgtaggggaa
841 gtaaaaagcg atgccaacag cgcagaagta gaacggtacg aacatgtacc agacatgata
901 gaccgttgtg aagggtttca gatctttctc gtccgggcac tggatggggg ttccgttttt
961 gaactgtctc cccgggaaac aaccgattgt tttcatttgg agaatacctg ggtaaggaat
1021 ttcatctggc caagacttgt catacacttc ttttaaagtg tacagccctt gagcagcaca
1081 atagtcttgc acaaatgtcg catgagcatc cgcaaacgaa ctgtcgtcgc ctcttgcatc
1141 cgatatacag tgtattattg atcctaaatt ctgtttaaag gtcatgattg taccagataa
1201 aatacaaaac atgaccaaat aggttctgtt aatctggtcc cactgatcgt cgatgctaac
1261 tgacttgatc tttaaaaagc cattcaagct tccagaaatg agatccaata ccat
//