Mercurial > repos > galaxyp > translate_bed_sequences
changeset 0:472ba4c3f53d draft
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/translate_bed_sequences commit e04ed4b4960d6109a85c1cc68a2bf4931c8751ef
| author | galaxyp |
|---|---|
| date | Wed, 13 Jan 2016 07:38:29 -0500 |
| parents | |
| children | 6c8e30f2e1f0 |
| files | test-data/Extract_Genomic_DNA.bed test-data/translated_bed_seqs.bed test-data/translated_bed_seqs.fa test-data/translated_bed_sequences.bed test-data/translated_bed_sequences.fa tool_dependencies.xml translate_bed_sequences.py translate_bed_sequences.xml |
| diffstat | 8 files changed, 596 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/Extract_Genomic_DNA.bed Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,7 @@ +15 40902459 40907576 JUNC00019210 1 + 40902459 40907576 255,0,0 2 36,38 0,5079 CTCATATGAAAATAGTGAGAAAGTCAGAAATGGAAGGTAAGTATGTTACTATAATTCCTGCTAAATCTGCTTTTTAAAAAATTTTACTGTAGCTATCTAACCAGACAACTAATATAACTTATTGTCTATCTTATTGCTTAAAAAATATATGATGGGCCATGCATAGTGGTTCACACCTGTGATCCCGACACTTTGGCAGTCTGAGGTGGGAGGATCGCTTGAGCCCAGGAGTTCAAGACCAACCTGGATAGCATAGTGAGATTCTGTCTCTACAAAAAATAAAAAATTAGCCAGGTGTGGGGGTGCATTCCTGTAGCCCCAGCTACTCAGGAGGCTGAGGTGAGAGGATCACTTGAGTCCAGAAGATTGAGACTGCAGTGATTGCACCAGTGCACTCCAGCCTGGGCAACAGAGCAAGACTCTGTCTCAAAAACAAAAAAAGTATACCTTGAGTCTTGGGAGTAGAGAGTAGAGTGGATAGTAAGTAAATTTATTTATTTATTTATTTTGAGACGGATTTTCACTGTTGTTCCCCAGGCTGCAGTGCAATGGTGCCATCTCAGCTGCAACTTCTGCCTCCTGGGTTTAAGCGATTCTCCTGCCTCAGCCTCCCAGGTAGCCGGGATTACAGGCATGCACCACCACACCCAGCTAATTTTTGTATTTTTAGTAGAGACAGGGTTTCACCATATTGGCCAGGCTGGTCTCGAACCCCTGGCCTCAAGTGATCCGCCTGCCTTGGCCTCAAAATGCTAGGATTACAGGCGTGAGCCACCACGTCCAGCAGTTAATAATTTTTTTAAATGTTGTTAATGTCCTTAAGTAATTTTTTAAAAATTTTATTGTATTTTTTTTTTTTTTGAGACAGAGTCTTGCTTGCCCTGTCGCCCAGGCTAGAGTGCAGTGGTGCGATCTCGGCTAACTGCAATCTCTGCCTCCCAGGTTCAAGCAATTCTCCTACCTCAGCCTCCTGAGTAGCTGGGATTACATGTGCACGCCACCACACCTGGCTAATTTTTGTATTTTTAGTAGAGACAGGGTTTCACCGTGTTGGTCAGGTTGATCTCGAACTCCTGACCTCGTGGTCCACCCGTCTCGGCCTCCCAAAGTTCTGGGATTACAGGCATGAGCCACCGGGCCCAGCTGTCCTTAAGTAATTTTTCTAATGCTTTATCAGTGTATCCTTGGACAAATTTAATTTTAATTTTATTTTTAGAAACAGAAACAGGAGAAAATCTTCTTTTGATACAGTAAGTAATATTCAAATTATATTTAAAAAGTTATAAATTAATATTATTTTTATTTCTTATTTTTGTGTTCCTTTTGTGGAATACTTAGGATTAATTCAGAAAACCAGTTATTTGCAGAATTTAGTATACTGCTATATCTGATTAATATGGCCATTGTACTGTTTATAGTTTTAATAGTGATGAAACAATATAATATATAAACCTTTATATTTATTAACTAATATAAAACAATAGATAGCAAGTCTTTAGTTTTTTTATTAAATTAATCTTTCTTATAATATTATTACATAGAAAGGCAAATACCAGTAGCTCTTCAACAAATGAAAAGATGGTTATCTTTATTCATCATAAGAGAAGTACAGATTTGATTATATTGAGATATTTTTGGCTGGACGTGGTGGCTCATGCCTGTAATCCGAGCACTTTGGGAGGCCGAGGCTGGTGGATCACCTGAGGTCAGGAGTTCGAGACCAGCCTGGCCAACATGATGAAACCCCATCTCTACTAAAAAAATACAAAAAATTGGCCGGGCACAGTGGCTTACGCCTGTAAACCAAGCACTTAGGGAGGCTGAGGCGGGTAGATCACGAGGTCAGGAGATCGAGACCATCCTGGCCAACGTGGTGAAACCCCATCTCTACTAAGAATACAAAAATTAGCCAGGTGTTGTGGTGGGCGCCTGTAGTCCCAGCTACTCAGGAGGTTGAGGTAGGAGAATCACTTGAATCTGGGAGGTAGAGGTTGCAGTGAGCCGAGATCATGCCACTGCACTCCAGCCTGGGCGACAGAGCGAGACTCCGTCTCAAAAAAACAAAAAAAACCCCAAAAAATTAGCCAGGTATGGTGGCACTGTGCCTGTAATCCCAGCTACTCAGGAGGCCAAGGCACGAGAATCGCTTGAACCCGGGAGGCGGAGGTTGCAGTGAGCTGAGATCATGCCCCTGCACTCCAACCTGGGCAGCAGAGCGAGACTCTGTCTCCAAAAAAGAGGTTTTTTTGCTGTTGTTGTTGTTTGTTTGTTGTTTAGACGGCATCTTCCTCTGTCGCCAGGCTGGAGTGCAATGGTACGATCTTGGCTCTCTGCAACCTCCGCCTCCTGGGTTCAAGTGATTCTCCTGCTTCAGCCTCCCGAGTAGCTGGGACTACAGGCGCCTGCCACCACGCCTGGCTAATTTTGTTGTATTTTTAGTAGACACGGGGTTTCACCATGTTGGCCAGAATTGTCTTGGTCTCTTGACCTCGTGATCCGCCCACTTTGGCCTCCTAAAGTGCTGGGATTACAGGCGTGAGCCACTGTGCCTGGCCTAAAAAAGAGATATTTTTATCCAGTCAGTTAAGGAAATAGGTAAAAAGTCACTCCGTGCTAGTGGGAGTATAATTTGTAGAAGGCAACTTGGTAATATCGTTCAGAATTTTAAATGCACATGCCCTTTGACTTCAGTTTTATGTCTGGAAATTATCTCACCAATATACGTGTAAATGTTTGAAATGATGAATAGATAAGGTTATTCACAGCAGCGTTTATTATAATAGTAAACGTTGTGAAAGGTCCTGAAGGTCCATCAGTAAGGGCCTGGTTAAATAAATTGTACATTTTCCATATAATTGACTATGCAGCCGGAAAAAAAATTTTATGCTAAATAGAGAAATCTCTAAGCTATATTACTTAAGATAAAAAGCAAGACGCAAAACAAAGTATATCATGTATTACCACTTGGGTAAAACAGGGAAGAAGAAAAGCGTGTGTGTGTGTGTGTATCTTCTCGTGTGTGTACATGTATAATGTATGTATATGTAGTGACTTATGCATATATATGCATAATATCTCTGGAAGAAACTAGTAAGATTGCCTACAGAGAAAAGAACCTTGCTATAGGTTACATGTTGGCTTTGTGGTTTATTAATTTAAAACCATATGAAAGTATTACCTATTACATTCTTAAAATTGTTGAAACGAAGAAATGGAAGTAATAGTAGTCATGGTTTAGAAAAGGTACAGAATCAACTCCAACAGAACTTAAAGGAAAAATAATTTTTTAAATAAAAAGGTACAGAAGTATATACATTGAAAAGTAAAAATTATTCCCTCTTTCCCTCTTCTTTCTTTTTTGTGAGACGGGGTTTCGCTCTTGTTGGCCAGACTGGAGTATAATGGCATGATCTTGGCTCACTGCAACCTCAGCCTCCTGGGTTCAAGCAATTCTCCTGCCTCAGCCTCCTGAGTAGCTAGGACTACAGGTGTGCACCACTACACCTGGCTAATTTATTTTATTTTTATTATTTTATTTTATTTTATTTTATTTTGAGATAGAGTCTCGCTCCGTCACCCAGGCTGGAGTGCAGTGGCGCCATCTCGGCCCACTGCAAGCTCCGCCTCCTGGTTTCACGCCATTCTCCTGCCTCAGCCTCCCTAGTAGCTGGGACTACAGGCGCCTGCCACCACGCCTGGCTCATTTTTTTGTATTTTTAGTAGAGACAGGGTTTCACTGTGTTAGCCAGGATGGTCTCAATTTCCTGACCTTGTGATCCACCCGCCTCGGCCTCCCAAAGTGCTGGGATTACAGGCGTGAGCCACCACGCCTGGCCCCCCTTTTTTTTTTTTTTGAGATGGAGTCTTGCTCTGTTGTGCAGGCTGGAGTGCAGTGGCCTAATCTTGGCTCACTGCAACCTCCGCCTCATGGGTTCAAGCAATTCTCCTGCCTCAGCCTCCCGAGTACCTGGGACTACAGGTGCACACCACCACGCCCAGCCAATTTTTGTATTTTTAGTAGAGATGGGGTTTCACCACGTTGGCCAGGCTGGTCTCGAACTCCTGACCTCAGTTTATCTGCCTGCCTTGGCTTCCCAAAGTGCTGGGATTTTAGGCATGAGCCACCACACCCGGCCTAATTTTGTATTTTTAGTAGAAGCGGTGTTTCACCATGTCGGCCAAGCTAGTCTCAAACTCCTGACCTCGGGTGATCTGCCTGCCTTGGCCTCCCAAAGTGATGGGATTACAGACGTGAGCCACCACACCCAGCCTATTCCCTCTTCTTCAGGTGACTCACTCCTTACAGATAACCACTGTTGACAATTTTTGGTGTACTTGCAGAAATTTTTTTTGTATGTCCAAGTGCATGTTTAGCTTTTATTTTTTTAAACAAGACAAAAATGGGATTGCCTAGTAATATGTCTTATAGCACTTTCTATGTTATTACATAAAGATTTACCTCAGTCTGTTACAATCTCTTGTATGGCTGCACCAACATTTAGTGAACTCCCTATTGATGGACACTTGAGTGGTTTCTTTCATGTACTTTTTGTTTTAAGGCTCAAGTTTTAGAAAAGGAATTGCTTTATAGGTCATTTTTTGAAAAAATAAGCAAGCAGTATATATTAAGCATATTTTTCTAGTTTATTTTTTTGAGATGGGGTCTCACTCTGTCACCCAGGCTGGAGTGCAGTGACATGATTATAGTTCATTGCAGCCCTGGACTTCAGGTCTCAAGTGATCCTCCCACTTCACGCTCCCAAAGTGTTGGGATTACAGATGTGAGCCACCACGGCTGGCTGATAATATTCTTGAAATTAGCAAACCATTAAGACCAATCTGAAAAGTATCCTTTAGTTGCCTTGTAATTCCTTCAGCCAGTATTAGATGAGTCTTTTTTTAAAAAGATGTCTTTATAATTTTCTGCATTAGATGAGTCTTTTTTAAAAAGATGTCTTTATAATTTTCTGCTAGTCATTCTATGAAAGTCTTTCACCTTCTGTACTCTTTCAGCTGTTATTGAAAGATTCTTGGTAATTGTTTGGGGTATAGATTATAAAAATACTTTTTAATTTTTAGGAATAAGAAATTAGAAGATAATTACTGTGAAATTACTG +9 17404 18053 JUNC00000003 1 - 17404 18053 255,0,0 2 75,26 0,623 GTCTCTGCCCCCAGCATCACATGGGTCTTTGTTACAGCACCAGCCAGGGGGTCCAGGAAGACATACTTCTTGTACCTACGGAGGCGACATGGGGGTCAGGCAAGCTGACACCCGCTGTCCTGAGCCCATGTTCCTCTCCCACATCATCAGGGGCACAGTGTGCACTGTGGGGTCCCAGGCCTCCCGAGCCGAGCCACCCCAGTCACCCCCTGGCTCCTGGCCTATGTGCTGTACCTGTGTCTGATGCCCTGGGTCCCCACTAAGCCAGGCCGGGCCTCCCGCCCACACCCCTCGGCCCTGCCCTCTGGCCATACAGGTTCTCGGTGGTGTTGAAGAGCAGCAAGGAGCTGACAGAGCTGATGTTGCTGGGAAGACCCCCAAGTCCCTCTTCTGCATCGTCCTCGGGCTCCGGCTTGGTGCTCACGCACACAGGAAAGTCCTTCAGCTTCTCCTGAGAGGGCCAGGATGGCCAAGGGATGGTGAATATTTGGTGCTGGGCCTAATCAGCTGCCATCCCATCCCAGTCAGCCTCCTCTGGGGGACAGAACCCTATGGTGGCCCCGGCTCCTCCCCAGTATCCAGTCCTCCTGGTGTGTGACAGGCTATATGCGCGGCCAGCAGACCTGCAGGGCCCGCTCGTCCAGGGGGC +8 27369375 27370080 JUNC00000874 1 + 27369375 27370080 255,0,0 2 52,49 0,656 GCCGACCTCTTGCAATCCAAGTGACATGAGCCATGGGTACGTGACAGTAAAGGTGAGTCAGTTTTGTCTCTCAGTCGGCTAAGTGCTCTCCCACCAGGTCACCTAAAACGACCAGCAGAGACACCCAAGAGGCTGAGCTGTGAGGATCACCTGAACCTGAGCCTGGGAAGTGGAGGTTGCAGTGAGCTGTGATCACACCACTGTGCTCCAGCCTGGGCAACGGAGTGAAACCCTGTCTCAAGAAAAAAAAAAAAAAAAAAGGACCAGCAGTGACATTTGTTAAATATCGAGGGTGGTTGAACATCCACTATTTATAAGGAAATGTTATTTCCCACAAATCTCATTCCTCAGAAATCAGTGAAAGACAGACCCTGTCTCGGATTCTATAAAGCAGTGTGACTGATGTGGCCAAACCCCCATAGAACAGCAGGCTAGGGAGCTTGGGCAGCCTGGGAGGCTAAAGCATGCACACAGCTGGCCTGTGTTGGAGGAGAGCATGTTTTCATTCAACAGGCATTTATCGAACACTTCCTCTGTGTCAGGCACAGGCCAGATACAAGGGGTGCAGTGACGATCAAGACAGACCTGTCCTCAAGGAGCCCACATTGAGGGCAAAGAGAGGGAAAGGTCATTAAAATGTTCTTTTTGTTTTAAAGGGTTATCATAAAGCCAAGGCCACACACAGAGGGCCCTGGCTCGTGGCCT +7 99077597 99079394 JUNC00001897 1 + 99077597 99079394 255,0,0 2 64,37 0,1760 CATTTATTGTCTTTATTTTTCCAACTTGGTGTGCATATCTGTATACCCAACTTCAGAAATAACAGTGTAGGCCAGGTGCGGTGGCTCACACCTGTAATCCCAGCACTTTGGGATTCCGAGGCGGGTGGATCACATGAGGTCAGTAGTTTGAGACCAGCCTGGCCAATATGGTGAAACCCAGCTCTACTTAAGATGCGAAAATTAGCCTGGTGTGGTGGTGCACGCTTGTAGTCCCAGTTACTTGGGAGGCTGAGGCAGGAGAATCGCTTGAACCCGGGAGGCAGAGGTTGCATTGAGCCGAGATCATGCCGCTGCATTCCAGCCTGGGTGACAGAGCGAGACTCCGTCTCAAAACAAAACAAAACAACAACAAAAATAACAGTGTATTTGAGGGCACCACTCCAAACCCTACTGGGGATGATACGTTATGTAAAATATTTACACCAGGTCACCTTTCTGAAATGGGAAAAATTCTGACTCCTGAAACACATCTGGTCCTAAGGGTTTCACATAAGGGACTGGGGACAGTGAATAGTCGTGTTGTGTTAGGGCTCAGACCGAGACACTTCATGCATTCTTCACTCTCACCTCGTGTTCCTCATCTGTAAGACCTCCAGGAGTTGAGAATTATGGTAGATTTACATTCCAAAGGCTCACTTCTTATGCTCTTGAAAGTCTTTTTCTTTCTTTTCTCTGGAAATGGGTAACTTCTTGGTGTGTAACATTAGTGGGTGCTCAGGAAATGATTTTCAACCAACTTCATTTCTGAACCTTCACTGCTTCCATTCTTTTCAGGAGCCTACCCTGTCCTAGTCCACTGTAGTCTCCCACTGTGGGAGTTGCCTTGCTAAATAGTTTTCTTTCCGGGGTCACATGTTTAACCCTTGCACTGACACTCTTAGCACAGAGTCTTGACCACCTCTGTGGGACTGTCCACCCACCTCTGGATTTTTATAATTTTAATTTTCCCTTTAAAAGCATTACAGGAGCATACATAAGTACATGGAAGGCATAGAAAAGGAGGAGAAGAGGAAAGTGTTCATGGGTAACACTGTTAATATTTTGGAAGGTTTCCTGCCAGTCTTTTTTTCCTATGCCTTTACAAAAGTAGCTGAGCATATTCTATGTGCAATCCTGTAATCTGTTTTAAATGTAACATTATTGTATCAGCATTTCCCCACATAATCAGTAAACGAGCATAATATAATTCCATTGTGTGGTCCATAAAAATATTTCATTGTTTTCCCATGTGGGGACATTGAGCTTGCCTTTGTTTGTTCGTTGGTTTGCTTGTGTATTTATAATGTTTGAATAGATGTTTGAGTATAAATCTTTGCATGAATTTTTGATCATTTCCCAAGATACTTTCCCAAAATCAGAATTTTGAAGTTAAATGAGATGGCACTTTAATCTTTATTGATTACTGAATGATTCTTCAACAAACTTTGTACTGTTTTGCATATCCAATAGTAACATTTTAGAGTGCTTCTTTTTTTTTTTGAGACAGAGTTTCGCTCTTGTTGCCCAGGCAGGAGTACAATGGCAGGATCTCAGCTCACCACAACCTCTGCCTCCCCGGTTCAAGTGATTCTCCTGCCTCAGACTCCTGAGTAGCTGGGATTACAGGCATGTGCCACCACACCTAGCTAATTTTGTATTTTTAGTAGAGACGGGGTTTCGCCATGTTGGTCAGGCTGGTCTCAAACTCCCGATCTCGGGTGACCCGCCCACCTCAGCCTCCCAAAGTGCTGGGATTACAGGTGTGAGCCACTGCACCCAGCTGAGAGTGCTTCTTTA +7 148909514 148910833 JUNC00002152 1 + 148909514 148910833 255,0,0 2 60,41 0,1278 GATCAGCAGGATTTGGCAGACAGAGATATTCCCACGGATCCCAATTCAGGTGAGAACAAGGTCAGAATGAATCTTGAGGGCAACAAGTGTTTTTCAAAAGGAATCATTAAATTAATAATTTATTGTGTTTATGCAGCTGATACTAAAATTACTTCAGCACTGATTTTTTGTTCTTTGAGACGGAGTCTTGCTCTGTTGTTCAGGCTGGAGTGCAGTGGTGTGATCTCAGCTCACCGCAACCTCCGTTTCCCGGGTTTAAGCGATTCTCCTGCCTCAGCCTCCTGAGTAGCTGAGATTACAGGCATGTGCCACCACACCTGGCTCATTTTTTGTATTTTTAGTAGAGATGGGGTTTCACCATGTTGGTCAGGCTGCTTCTGAACTCCTGGCCTCAAATGATCTGCCCGCCTCAGCCTCCCAAAGTGCTGGGATTACGGGCATGAGCCACCGCGTCCGGCCATTCAGCACTGATTTTTGAACAATAAAAACAACCAAAAAGCATCTAAAAGCTGCATATATTTCTGTGTCAGAAAATTATTTACATTTTCTGGAATTTCTCTTGAGTGAATACTAGTTAGAAGTTTCCCAAGAAACTTGTAGGGAGGGTCTAGGCATTCAGGGAGGCATAAAAGAACAATGGAAGCCCCACCTGCCCACGACCTATGTGTGTGTCCATGTGATTGCAGAGCTGAGAAGTCCAGTGTCATTGCCTGGTAGTTAGGACACAAGTCTCATCCTACCCACCATACCAGGGAGAGAGGATGAGAAGTTGAATCTAGGATTCCTTGTACTTTGTCAACTTCACCACAGAGAAATTGAGGGGACCCAGGGCAGAGGGAAAGCGGGGCCTTCCTGCTTCACTGTGGTTCTTATCAATGACCTCTGGTTGCAGCCTCTGCCATCGCTGGTTTCTAAATCATACCTGTTAGCCATGGCCCTGAGGGCCGTTATCTGGGGCTTTCTTCTGGCCTGTCCTGGCCCCATGGGTTCCTCTGGGGCTAAGTGGCAGAGCTGGGAGATCTGGGTTTTAATCCCGGCCCTGCCACTGCACAGACCTTGGGCTAAGCCTCCTTCTGGCCAACTTGAAATGCATCGAGCAGTGAGAGTCCATTTAGGAAGGGGCTTATGAGGGTAAGGAGGATGGCCTGGAAATGCAACACCGATACCAAAACATGCCCTGCAGGTTGCAGCCACCGCATGGGATGGTTTTCAAACCTTTGATGGGAGGCCGAGTAGAACTGAACAAAATCTAAGACTGCCTTCTTTCCATGACAACAGAGTCTCTCATCTCAGCACATGACATTTTGTCATGGATCAAG +6 41766612 41767580 JUNC00002625 1 - 41766612 41767580 255,0,0 2 55,13 0,955 CTGGTCAAAGACGACATGGACCCCAATCTTCTCTCGATGATTACGGCCAGACCACCTAGAACATGGATATAAGCTTGTTACTAGTAAAAACTGGCATATATAATCGTAGGTAGAGGTAGGTATTAATAAAACAGGAAACTGAGAAAATTTAACCTTCTGGTTTTGGTTTTCCCATCAATAAAATGGTGGTAATCAATGTTCACACTTATTTATTTATTTATTTATTTATTTATTTATTTATTTTCAGACAGGGTCTCACTCTGTTGCCCAGGCTGCAGTGCAGTGGTGTGATTGATAAAGGCTCACTGCAGCTTCAACCTTCCCGAGCTCAGGTGATCCCCCCACCTCAGCCTCCCGAGTAGGTGGGACTATAGGTGCACACCACCATGCCTGGCTAATTTTTGCATTTTTTGTAGAGACAGGGTTTTAACCATGTTGCCGAGGCTGGTCTCAAACTCCTGAGCTCAAGCAATCTGCCTACCTTGGCCTCCCAAAGTTCTAGGGTTACAGTCATGAGCCACTGCGCCTGGCCCATTTATTTGCTGTTAATGTTCCATGAAGTACTCTAAACCTAGAGATCCCTAAACAAGCAATTCAAGCAACAGTCTTTTTGTAGTTTCCTGCCTGTGCTTCTTTTGCAGATTATTTTATACTTCATCCTGAAGATGAATCTACATAAGTCCAAGCCAGGCTAGAAAATCTTTAAGCACAGACACCCAGTGCAGGCTGGATAATTAATTACACAGGAGCATGAGAGCAAAAGCCCAGTCCTGCTGGGCTCCCACATCAATTTGCCTTAGTGGATGTGAACAGTTTAATCTCCATGGAAAACACGAGGTCTCAAGTGAACAGGCTTCCCCTCCCAACACCATAGAAGACATTTGCCCTCACCAAAGCCTTTGGAGCAGCTATTCCCACAGCTCTCACCCAGCTGGGACTCCCACAGCAAGCTCACCTGAATCTTTTAA +6 84856456 84862316 JUNC00002772 1 - 84856456 84862316 255,0,0 2 48,53 0,5807 TTTACTCTGCAGCTCATTAACTTGACTTTGGAAATGTCTTATAAGTACCTAGGAAATTGTTAGAAGTGCTGAGTTAATAAAGCTCCTCCTTGAACACTATCTGAAAAATATAAGAAAGCCAACTCACAAAAGGGAACTCTTATATACTGTTGGTGGGAATGTAAGTAAGCACAGCTACTATGGAAAAAAATTATAGAGATTTCTCAAAAAACTAAAAATAGAATTACCGTATTATCTAGCAATCCCACTACTGGGTATCTATCCAAAGGAAAAGAAATTTGTATATCAAACGGACACCTGGCACTTGCATGTTTATTGCAGCACTATTCACAATAGCAAAGATATGGAATCAATCTAAGTTCCCATCAACAGACGAAAAGAAATGTGATATATATACCCAATGGAATACAATTCAGCCTTAAATAAGAGGGAAATAATGTCATACGCAGTTAAATGGATGGAACTGAAGGTCATTATGTTAAGTGAAGTAAGCTAGGCACAGAAAAATACTGCATGTTCTCACTCATGTATGAGAGCTAAAAAACTTGATCTCATGAACACAGAGAATAGAACAATAGATACCAGAGATGGGAGTGGGGGGTTGGGAAGGATGAAGAGAGGTTGGTTATGGGTACAAACATACATACAGTTAGATAGATGAAATAAGTTCTAATGTTTAATAGCAGACCAGGGTGACTGTACTTGGCAACAATATTATGTGTATTTCACGGTAACTAGAAGAGAGAGTTTGAAGTGATACCAAGACATAAAAATGATGAACACTCAAGGTGACAGAAATCCCAATTACCCTGACTTGATTATTACACATTCTTTGCATGTAGCAAACACTCACATGTATCCTATAAATATCTAAATCATTATATATCAATAAAAGGGGGGGAAAATAAATATGAACATTTTTCATGACTAGAAAAAGTCCTTTTAAAATAAAGCTCTATAAAGTCTACACATTTGGAAAAAGGTAGACAATATTATCTTTATACTCTGTGGTTAAGTATGATAGTTTTTATTTAAGAACTATGGCTTAAATATAACATAAGCCTCATAATATCTGTGAAATAAATGTCAGTCTCATATCTTTAGGTAAACTGAAATATATGAACAGATACTTCCATCAAGAAACTAAGAGAAATTGGGAAAAATTTATCTTATCTCTAATGCTGTATTTTTAATGATCAGACTATGCTGATTTCTCTTTAACTCAATTTTACAATATTACATTTATAACGCTGGGAAAATAACATTTATGGCACTATGGCATCATGATACTACTTTTCCTACCTCTTCCAGACTATTCCTTAATACAAAATCACTTGATAGAATGCTTATAGTTGCCCTCTGACACCCTATTCCTATGGAGGATAAGGCCAGAGAGGTTCTCTGTAAGGCTGTATGGATTGTTACTTACGAAAGGCACTTTAAGTATTTTTCATTAACGTAGTCAAATAGTCCTTTTTTTTCTCAGCGATCTGAAGTTAAGTTACAGACTTCATTACATTTTGCCCCTAAATACTTCAAATGCCACTTTGGGAGGTCAAAGTGGGTGGATCACTTGAGGTCAGAAGTTCGAGACCAACCTGGCCAACATGGTGAAACCCCATCTCTACCAAAAATATAAAAATTAGTCAGGCATGGTGGCGTGCACCTGTAATCCCAGCTACTTGGAAGGCTGAGACAGAAGAATCGCTTGAACCTGGGAGGCAGAGGTTGCAGTGAGCTGAGACAGCCTGGGTGACAGAAGGAGACTGTCTCAAAGAAACAAAACCTTCAGATGCTTCAACATACATCTCTTTAGAATAAGGGCATTTCTCACAAAACCATAACATTTTATACCTAAGAAAATTTATAGTAATTTCTTAGTATCTAATAGCCAGTCCATATTCAAACTTCCCGAATTGAGGAATTTATTTTTAATTAAGATGGTATTAAATAATTACACTTGGTGAGATAATTCTTTCTTGATGTTTTTCTCCCTGACCTTCTTTCTTAAAATTGAAGCCTACAGACAAGATGAACACAATAAACATTCATATATTTTTCAACTTGATAACCCAAAAGGTGGCTGATTTGCCACCTTTTTCTTTTTCTTGGATATACCATTTGAAAGTAGATTGTAGATAAAAGTCCTATTTTAATGCTATTTTTATAAGTGACAAACAGTACAGAAGACAGGCCACATGGTAATGCTCATAGATATTCTAAGTAGATATACCTATTTTTCAAAAAGGAAGACACATTTATTTGAGGAAATGCTGATGAGTTGGCAGCATTCCCTATAATAATTCTTAGAAGACTATGCTCCATCAGAGCAAGAGCTGTGTATCTCTCTTTTCATTATTTTATTCCAGTGCCTGATTGTCTATTGACCTGAATCATTATTTGTTAATAAAAGGTCATGGGGCAGGAAGAATGCCTAAGGATAACTGTCAGCTGGGACTTGGGGACTCGGTTCAACATCCAAGTTCAAGTACTAGGGTCCTAATCATGTTAGACATATACCTCACCACCTGTGACAAATTTAGAGACTTATTTGAGAAGTATAGGGACAAACTGTAAAAATATTTATTTATATGGTGTAGTTTTAAAGTAGAAAAATTAAATATATAGTATTATCTTCTTTAAAAAAATTAAATATATAGTATTATCTTCTTTAAAAAAGAAATCAAGCCAATAATAGTTTCCCATTTCTATCACCAAGCTGAAAAAGTTAAAACATTTCCAAATAAAACAATGTTAACTTTAATCTCATTAAAGTCAAAACGAGTTTATTCAAGTCAAAAGATAAAACAGATTTGTCATCTACCTCTTGAGTTGCTATTTTACGATTTAGTTCAGCTACTTTGGAAGAAGAATTTTCTACTGCATTTTGGCAAAGGAGTTTCTCTATTTCTCTCTGATGAGATGCTTTGAGGGATGCAATGTGTGCTGCAGTATCTTCCTTGACCCTACAGAGCAAATGAAAGAAAACCACACATAAAAGTGGCTCTTCAACCTCTAATTCAGAGTTAGCACAAAAACAGTAGAAGAAACCTGGGAATTAAGCAAGGGGAAAATGGGCACCAATATTAACTGCACTGTGTTAGGTGCTTTATACATTACTACTGCCTATGAGTTCTTATGACCAATCCATGAGCTACTATTTACTGAATTTTATAAGGGCAAACTTGAGGAAATAGTGATATATCTCTACTTTACCCCAACCTGAGGCTACACTTTCTGAGGGAATCTAGGGTTCCAAAAGCAAAGCACTTGACATTTGATCTTTTAAGATAAACTCAGATTTCTTGCCCAGATGCTTCTTTGGAAAATCACAGCTGTTCCCTTTAGAATCTCCAATGCCTTGTAGGCAGCATTAACCAATTTTCAACCCATGCATGAGACTCCAGGGCTAGTTGCATATTGTGAGAATGAGACAATCTTTATGGTAATTCTCAGGACATATTATGCTTGTAATATGCAATTAATTAATTGGACTGTAAGTGTGAAAGACTATTCATTTGTCTCTGTGATAAAATTTATTCAGAGGCAGTTCAGTAAATGGTTAAAATTATAAATAAAGATTTGGTATATTTCCTACAGATAAACATTTTCTCTCTTATTGGAAAGACTTTGAAAATAACAGCATTTAACTTAAAGCATGTATTTTAGTAAACTGCTGTTGCTATATATTAATATGGTGCTTCCCAACTATAGACAAATTTATCCATTGCTTATTTCACCATATGAATATTTCTCCATGGAAACAGTGTGAGTAGGGTCCTAGAAAAAGCATAAAAACCTAAATTAAAGCACATATAACTGAATGATAGTAAAAAGTCTAAATATTGGACTAGAAAAGTTTATTCTTCACTTCTCAGGCCCAGACTTATCTAAGAAACAGATTGACAATATTAACAAATTTTAAATGATGTACAGGAAAATTAATTGGAAAACCTGAATTTTAAATTCTTTCCTGGAAAAAATGTGATTCTTCTTAAAGAGCTTCATGCTACTGTCATAATATTCAGGACTTAAGAAACAAGCTTTCCATTTGGTAGCAATATGACTGCCATATGGTAACAACATGACTCACTCACTGCCATGGCTGCCAGTCAATAAGCATTTGTTGGGTGTTTCCCATGCTCATTCTCCTCCTACAGCACATTGCTGTGGGGATGCAAAAGCAGCATAATCTTACACCAGTTCTCAAGGAACTACAAGGCTGTGTATGATTAATCCTTTCTTGGTACATATTAGTAAGTTATGTAGGAGTTGGGGGTAGGAGAAGGCAATTTCCTGGGAAGTTCCACGAAATGTTTTGTGAAGGTAGTACTTATAACAAGAATGAGGAAAAGGCATGAGGAAATGTGACAGATTAGGTAAGTAGAGTAGTATTAGAATAGTGATGAGACTTGAGATGCTAGGCAGAGAAATTTAGTCTTGATATTTAAGGAAGTTGAGAGTCACCGGAGAGCAGAGTGACCCAGGTATATGAGTATTTAAAGAATATTAATGATTTTTAAGCAAGATGACTGGGGAAGGTTAAAAAAAAGGAAGATTGGAGTGGGAAAATTGTATAATGGCAGCACTGGTTTATAAAAGTAAGGAAGACTGGAAAATGGGAAACATTTATACCCAAGGAATAACAGAGATAAAGCTGGGCTGGTAATTAAAATAAATATAGACAGGGAAGGGGAGCTTGGGAAGGGGAACTAATTTACAGGTAAAAAAAGGAAACTTGAGGTATTAGTGATTCCAAGAGAAGCAATCCAATAAATAGTTCTAGGTTTCACATCTGAACTCAGCAAAGAGAACAAGACTTTAAAAGTAAATCTAGAAGTCATGGTTCACAGAAACCGTAAGAATGGGTAAGTTTAGCAACAGAGAAAGAACAGAGAATAAACAGAAAACTAAGAGCTGAATGCTGGTAGAAAGCCAGATTTTTGGGGATGGAGAAGGAAGAAGGGTCAGAATTCAGTAAAAGAAACAGAAAAGGAGTTTTGGTCAGAGCAGCCAAGCCACAGAGGTTAAGGGGAGAAGCCATTTCAGGAAAAGTAGTGTGGTCAACATTGTTGAAGAAAATGAGGAAGGAGAAAATGCCCCTGGACTGGGAAAGCTCATTATTGAGACCAGGTTTAAAAAGGGTGATGATGAGAGCCATATTTAAGAGGAAGTGAGATAGAAATAAAAGAAAATACTATTTGTGAAAAACGGGTGCTAAAAGACAGGATTGGATTGGAACCTGATCAGGTGAAAAAAAGATCTGAGGATGTCAAAAGACTCAGAGGAAAGAGACAGTGGATAGCAAATGGAGATAGTGAGAAGGAAATGGATGAGTATAAAGTTCCTTAAGGCATGAATGGATGTGATCACAAGGTCAATTTGAAGAGTTATTTTTGAAAAGAGGACCTTTTCCTTTTAGACTGCCAACTAGCTTTAACAATAACTGTTTTCTGAAAGCAATCTGTCTATAGCAATAATAAACTATGACATTCCTTTGCTAGAAGAGCCCAAAGAGTCAGCTTATTAATACCATTTGGGTACATGCCAAAAGGCAAGACTGGCGAGCCCAAGGACACAATTTGGGTTGACACAGATGTTCAGGTACATCAGAATTTGGGCCACCAGCAGGTAGAAGACCCTATCCTGAAACAGACACATGCTCAGCAATATAATATAATTTATTAATATTCACTTATGAATTCAGTTCGGGGTCAAACTAATTATCTGGAGGAAATATATCGTATAATTTTTCTATTTTTATCTTTTTATTTTTACAAATATTTATAAATAATTTAACATTTTACCTTCTCATGGAGTTTTCAAATTGATTCATCACTGCTTCAGATTTCATCTTCAG
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/translated_bed_seqs.bed Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,14 @@ +track name="novel_junctioni_translations" description="test" +15 40902460 40907575 JUNC00019210_2 1 + 40902460 40907575 255,0,0 2 35,37 0,5078 SYENSEKVRNGRNKKLEDNYCEIT +15 40902461 40907549 JUNC00019210_3 1 + 40902461 40907549 255,0,0 2 34,11 0,5077 HMKIVRKSEMEGIRN +9 17406 18053 JUNC00000003_1 1 - 17406 18053 255,0,0 2 73,26 0,621 APWTSGPCRYKKYVFLDPLAGAVTKTHVMLGAE +9 17438 18052 JUNC00000003_2 1 - 17438 18052 255,0,0 2 41,25 0,589 PPGRAGPAGTRSMSSWTPWLVL +9 17404 18051 JUNC00000003_3 1 - 17404 18051 255,0,0 2 75,24 0,623 PLDERALQVQEVCLPGPPGWCCNKDPCDAGGRD +8 27369376 27370079 JUNC00000874_2 1 + 27369376 27370079 255,0,0 2 51,48 0,655 PTSCNPSDMSHGYVTVKGYHKAKATHRGPWLVA +7 99077597 99079392 JUNC00001897_1 1 + 99077597 99079392 255,0,0 2 64,35 0,1760 HLLSLFFQLGVHICIPNFRNNSVSHCTQLRVLL +7 99077598 99079360 JUNC00001897_2 1 + 99077598 99079360 255,0,0 2 63,3 0,1759 IYCLYFSNLVCISVYPTSEITV +7 148909514 148910831 JUNC00002152_1 1 + 148909514 148910831 255,0,0 2 60,39 0,1278 DQQDLADRDIPTDPNSGENKSLSSQHMTFCHGS +7 148909515 148910811 JUNC00002152_2 1 + 148909515 148910811 255,0,0 2 59,19 0,1277 ISRIWQTEIFPRIPIQVRTRVSHLST +6 41766614 41767580 JUNC00002625_1 1 - 41766614 41767580 255,0,0 2 53,13 0,953 LKDSGGLAVIIERRLGSMSSLT +6 41766612 41767578 JUNC00002625_3 1 - 41766612 41767578 255,0,0 2 55,11 0,955 KRFRWSGRNHREKIGVHVVFDQ +6 84856497 84862316 JUNC00002772_1 1 - 84856497 84862316 255,0,0 2 7,53 0,5766 LKMKSEAVMNQFENSMRRYL
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/translated_bed_seqs.fa Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,26 @@ +>JUNC00019210_2 pep:novel chromosome:GRCh37:15:40902460:40907575:1 depth:1 +SYENSEKVRNGRNKKLEDNYCEIT +>JUNC00019210_3 pep:novel chromosome:GRCh37:15:40902461:40907549:1 depth:1 +HMKIVRKSEMEGIRN +>JUNC00000003_1 pep:novel chromosome:GRCh37:9:17406:18053:-1 depth:1 +APWTSGPCRYKKYVFLDPLAGAVTKTHVMLGAE +>JUNC00000003_2 pep:novel chromosome:GRCh37:9:17438:18052:-1 depth:1 +PPGRAGPAGTRSMSSWTPWLVL +>JUNC00000003_3 pep:novel chromosome:GRCh37:9:17404:18051:-1 depth:1 +PLDERALQVQEVCLPGPPGWCCNKDPCDAGGRD +>JUNC00000874_2 pep:novel chromosome:GRCh37:8:27369376:27370079:1 depth:1 +PTSCNPSDMSHGYVTVKGYHKAKATHRGPWLVA +>JUNC00001897_1 pep:novel chromosome:GRCh37:7:99077597:99079392:1 depth:1 +HLLSLFFQLGVHICIPNFRNNSVSHCTQLRVLL +>JUNC00001897_2 pep:novel chromosome:GRCh37:7:99077598:99079360:1 depth:1 +IYCLYFSNLVCISVYPTSEITV +>JUNC00002152_1 pep:novel chromosome:GRCh37:7:148909514:148910831:1 depth:1 +DQQDLADRDIPTDPNSGENKSLSSQHMTFCHGS +>JUNC00002152_2 pep:novel chromosome:GRCh37:7:148909515:148910811:1 depth:1 +ISRIWQTEIFPRIPIQVRTRVSHLST +>JUNC00002625_1 pep:novel chromosome:GRCh37:6:41766614:41767580:-1 depth:1 +LKDSGGLAVIIERRLGSMSSLT +>JUNC00002625_3 pep:novel chromosome:GRCh37:6:41766612:41767578:-1 depth:1 +KRFRWSGRNHREKIGVHVVFDQ +>JUNC00002772_1 pep:novel chromosome:GRCh37:6:84856497:84862316:-1 depth:1 +LKMKSEAVMNQFENSMRRYL
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/translated_bed_sequences.bed Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,14 @@ +track name="novel_junctioni_translations" description="test" +15 40902460 40907575 JUNC00019210_2 1 + 40902460 40907575 255,0,0 2 35,37 0,5078 SYENSEKVRNGRNKKLEDNYCEIT +15 40902461 40907549 JUNC00019210_3 1 + 40902461 40907549 255,0,0 2 34,11 0,5077 HMKIVRKSEMEGIRN +9 17406 18053 JUNC00000003_1 1 - 17406 18053 255,0,0 2 73,26 0,621 APWTSGPCRYKKYVFLDPLAGAVTKTHVMLGAE +9 17438 18052 JUNC00000003_2 1 - 17438 18052 255,0,0 2 41,25 0,589 PPGRAGPAGTRSMSSWTPWLVL +9 17404 18051 JUNC00000003_3 1 - 17404 18051 255,0,0 2 75,24 0,623 PLDERALQVQEVCLPGPPGWCCNKDPCDAGGRD +8 27369376 27370079 JUNC00000874_2 1 + 27369376 27370079 255,0,0 2 51,48 0,655 PTSCNPSDMSHGYVTVKGYHKAKATHRGPWLVA +7 99077597 99079392 JUNC00001897_1 1 + 99077597 99079392 255,0,0 2 64,35 0,1760 HLLSLFFQLGVHICIPNFRNNSVSHCTQLRVLL +7 99077598 99079360 JUNC00001897_2 1 + 99077598 99079360 255,0,0 2 63,3 0,1759 IYCLYFSNLVCISVYPTSEITV +7 148909514 148910831 JUNC00002152_1 1 + 148909514 148910831 255,0,0 2 60,39 0,1278 DQQDLADRDIPTDPNSGENKSLSSQHMTFCHGS +7 148909515 148910811 JUNC00002152_2 1 + 148909515 148910811 255,0,0 2 59,19 0,1277 ISRIWQTEIFPRIPIQVRTRVSHLST +6 41766614 41767580 JUNC00002625_1 1 - 41766614 41767580 255,0,0 2 53,13 0,953 LKDSGGLAVIIERRLGSMSSLT +6 41766612 41767578 JUNC00002625_3 1 - 41766612 41767578 255,0,0 2 55,11 0,955 KRFRWSGRNHREKIGVHVVFDQ +6 84856497 84862316 JUNC00002772_1 1 - 84856497 84862316 255,0,0 2 7,53 0,5766 LKMKSEAVMNQFENSMRRYL
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/translated_bed_sequences.fa Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,26 @@ +>JUNC00019210_2 pep:novel chromosome:GRCh37:15:40902460:40907575:1 depth:1 +SYENSEKVRNGRNKKLEDNYCEIT +>JUNC00019210_3 pep:novel chromosome:GRCh37:15:40902461:40907549:1 depth:1 +HMKIVRKSEMEGIRN +>JUNC00000003_1 pep:novel chromosome:GRCh37:9:17406:18053:-1 depth:1 +APWTSGPCRYKKYVFLDPLAGAVTKTHVMLGAE +>JUNC00000003_2 pep:novel chromosome:GRCh37:9:17438:18052:-1 depth:1 +PPGRAGPAGTRSMSSWTPWLVL +>JUNC00000003_3 pep:novel chromosome:GRCh37:9:17404:18051:-1 depth:1 +PLDERALQVQEVCLPGPPGWCCNKDPCDAGGRD +>JUNC00000874_2 pep:novel chromosome:GRCh37:8:27369376:27370079:1 depth:1 +PTSCNPSDMSHGYVTVKGYHKAKATHRGPWLVA +>JUNC00001897_1 pep:novel chromosome:GRCh37:7:99077597:99079392:1 depth:1 +HLLSLFFQLGVHICIPNFRNNSVSHCTQLRVLL +>JUNC00001897_2 pep:novel chromosome:GRCh37:7:99077598:99079360:1 depth:1 +IYCLYFSNLVCISVYPTSEITV +>JUNC00002152_1 pep:novel chromosome:GRCh37:7:148909514:148910831:1 depth:1 +DQQDLADRDIPTDPNSGENKSLSSQHMTFCHGS +>JUNC00002152_2 pep:novel chromosome:GRCh37:7:148909515:148910811:1 depth:1 +ISRIWQTEIFPRIPIQVRTRVSHLST +>JUNC00002625_1 pep:novel chromosome:GRCh37:6:41766614:41767580:-1 depth:1 +LKDSGGLAVIIERRLGSMSSLT +>JUNC00002625_3 pep:novel chromosome:GRCh37:6:41766612:41767578:-1 depth:1 +KRFRWSGRNHREKIGVHVVFDQ +>JUNC00002772_1 pep:novel chromosome:GRCh37:6:84856497:84862316:-1 depth:1 +LKMKSEAVMNQFENSMRRYL
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_dependencies.xml Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,6 @@ +<?xml version="1.0"?> +<tool_dependency> + <package name="biopython" version="1.62"> + <repository changeset_revision="6b6a336db91a" name="package_biopython_1_62" owner="biopython" toolshed="https://testtoolshed.g2.bx.psu.edu" /> + </package> +</tool_dependency>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/translate_bed_sequences.py Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,375 @@ +#!/usr/bin/env python +""" +# +#------------------------------------------------------------------------------ +# University of Minnesota +# Copyright 2014, Regents of the University of Minnesota +#------------------------------------------------------------------------------ +# Author: +# +# James E Johnson +# +#------------------------------------------------------------------------------ +""" + +""" +Input: BED file (12 column) + 13th sequence column appended by extract_genomic_dna +Output: Fasta of 3-frame translations of the spliced sequence + +""" + +import sys,re,os.path +import tempfile +import optparse +from optparse import OptionParser +from Bio.Seq import reverse_complement, transcribe, back_transcribe, translate + +class BedEntry( object ): + def __init__(self, line): + self.line = line + try: + fields = line.rstrip('\r\n').split('\t') + (chrom,chromStart,chromEnd,name,score,strand,thickStart,thickEnd,itemRgb,blockCount,blockSizes,blockStarts) = fields[0:12] + seq = fields[12] if len(fields) > 12 else None + self.chrom = chrom + self.chromStart = int(chromStart) + self.chromEnd = int(chromEnd) + self.name = name + self.score = int(score) + self.strand = strand + self.thickStart = int(thickStart) + self.thickEnd = int(thickEnd) + self.itemRgb = itemRgb + self.blockCount = int(blockCount) + self.blockSizes = [int(x) for x in blockSizes.split(',')] + self.blockStarts = [int(x) for x in blockStarts.split(',')] + self.seq = seq + except Exception, e: + print >> sys.stderr, "Unable to read Bed entry" % e + exit(1) + def __str__(self): + return '%s\t%d\t%d\t%s\t%d\t%s\t%d\t%d\t%s\t%d\t%s\t%s%s' % ( + self.chrom, self.chromStart, self.chromEnd, self.name, self.score, self.strand, self.thickStart, self.thickEnd, self.itemRgb, self.blockCount, + ','.join([str(x) for x in self.blockSizes]), + ','.join([str(x) for x in self.blockStarts]), + '\t%s' % self.seq if self.seq else '') + def get_splice_junctions(self): + splice_juncs = [] + for i in range(self.blockCount - 1): + splice_junc = "%s:%d_%d" % (self.chrom, self.chromStart + self.blockSizes[i], self.chromStart + self.blockStarts[i+1]) + splice_juncs.append(splice_junc) + return splice_juncs + def get_exon_seqs(self): + exons = [] + for i in range(self.blockCount): + # splice_junc = "%s:%d_%d" % (self.chrom, self.chromStart + self.blockSizes[i], self.chromStart + self.blockStarts[i+1]) + exons.append(self.seq[self.blockStarts[i]:self.blockStarts[i] + self.blockSizes[i]]) + if self.strand == '-': #reverse complement + exons.reverse() + for i,s in enumerate(exons): + exons[i] = reverse_complement(s) + return exons + def get_spliced_seq(self): + seq = ''.join(self.get_exon_seqs()) + return seq + def get_translation(self,sequence=None): + translation = None + seq = sequence if sequence else self.get_spliced_seq() + if seq: + seqlen = len(seq) / 3 * 3; + if seqlen >= 3: + translation = translate(seq[:seqlen]) + return translation + def get_translations(self): + translations = [] + seq = self.get_spliced_seq() + if seq: + for i in range(3): + translation = self.get_translation(sequence=seq[i:]) + if translation: + translations.append(translation) + return translations + ## (start,end) + def get_subrange(self,tstart,tstop): + chromStart = self.chromStart + chromEnd = self.chromEnd + r = range(self.blockCount) + if self.strand == '-': + r.reverse() + bStart = 0 + for x in r: + bEnd = bStart + self.blockSizes[x] + if bStart <= tstart < bEnd: + if self.strand == '+': + chromStart = self.chromStart + self.blockStarts[x] + (tstart - bStart) + else: + chromEnd = self.chromStart + self.blockStarts[x] + (tstart - bStart) + if bStart <= tstop < bEnd: + if self.strand == '+': + chromEnd = self.chromStart + self.blockStarts[x] + (tstop - bStart) + else: + chromStart = self.chromStart + self.blockStarts[x] + self.blockSizes[x] - (tstop - bStart) + bStart += self.blockSizes[x] + return(chromStart,chromEnd) + #get the blocks for sub range + def get_blocks(self,chromStart,chromEnd): + tblockCount = 0 + tblockSizes = [] + tblockStarts = [] + for x in range(self.blockCount): + bStart = self.chromStart + self.blockStarts[x] + bEnd = bStart + self.blockSizes[x] + if bStart > chromEnd: + break + if bEnd < chromStart: + continue + cStart = max(chromStart,bStart) + tblockStarts.append(cStart - chromStart) + tblockSizes.append(min(chromEnd,bEnd) - cStart) + tblockCount += 1 + ## print >> sys.stderr, "tblockCount: %d tblockStarts: %s tblockSizes: %s" % (tblockCount,tblockStarts,tblockSizes) + return (tblockCount,tblockSizes,tblockStarts) + ## [(start,end,seq,blockCount,blockSizes,blockStarts),(start,end,seq,blockCount,blockSizes,blockStarts),(start,end,seq,blockCount,blockSizes,blockStarts)] + ## filter: ignore translation if stop codon in first exon after ignore_left_bp + def get_filterd_translations(self,untrimmed=False,filtering=True,ignore_left_bp=0,ignore_right_bp=0,debug=False): + translations = [None,None,None,None,None,None] + seq = self.get_spliced_seq() + ignore = (ignore_left_bp if self.strand == '+' else ignore_right_bp) / 3 + block_sum = sum(self.blockSizes) + exon_sizes = [x for x in self.blockSizes] + if self.strand == '-': + exon_sizes.reverse() + splice_sites = [sum(exon_sizes[:x]) / 3 for x in range(1,len(exon_sizes))] + if debug: + print >> sys.stderr, "splice_sites: %s" % splice_sites + junc = splice_sites[0] if len(splice_sites) > 0 else exon_sizes[0] + if seq: + for i in range(3): + translation = self.get_translation(sequence=seq[i:]) + if translation: + tstart = 0 + tstop = len(translation) + offset = (block_sum - i) % 3 + if debug: + print >> sys.stderr, "frame: %d\ttstart: %d tstop: %d offset: %d\t%s" % (i,tstart,tstop,offset,translation) + if not untrimmed: + tstart = translation.rfind('*',0,junc) + 1 + stop = translation.find('*',junc) + tstop = stop if stop >= 0 else len(translation) + offset = (block_sum - i) % 3 + trimmed = translation[tstart:tstop] + if debug: + print >> sys.stderr, "frame: %d\ttstart: %d tstop: %d offset: %d\t%s" % (i,tstart,tstop,offset,trimmed) + if filtering and tstart > ignore: + continue + #get genomic locations for start and end + if self.strand == '+': + chromStart = self.chromStart + i + (tstart * 3) + chromEnd = self.chromEnd - offset - (len(translation) - tstop) * 3 + else: + chromStart = self.chromStart + offset + (len(translation) - tstop) * 3 + chromEnd = self.chromEnd - i - (tstart * 3) + #get the blocks for this translation + (tblockCount,tblockSizes,tblockStarts) = self.get_blocks(chromStart,chromEnd) + translations[i] = (chromStart,chromEnd,trimmed,tblockCount,tblockSizes,tblockStarts) + if debug: + print >> sys.stderr, "tblockCount: %d tblockStarts: %s tblockSizes: %s" % (tblockCount,tblockStarts,tblockSizes) + # translations[i] = (chromStart,chromEnd,trimmed,tblockCount,tblockSizes,tblockStarts) + return translations + def get_seq_id(self,seqtype='unk:unk',reference='',frame=None): + ## Ensembl fasta ID format + # >ID SEQTYPE:STATUS LOCATION GENE TRANSCRIPT + # >ENSP00000328693 pep:splice chromosome:NCBI35:1:904515:910768:1 gene:ENSG00000158815:transcript:ENST00000328693 gene_biotype:protein_coding transcript_biotype:protein_coding + frame_name = '' + chromStart = self.chromStart + chromEnd = self.chromEnd + strand = 1 if self.strand == '+' else -1 + if frame != None: + block_sum = sum(self.blockSizes) + offset = (block_sum - frame) % 3 + frame_name = '_' + str(frame + 1) + if self.strand == '+': + chromStart += frame + chromEnd -= offset + else: + chromStart += offset + chromEnd -= frame + location = "chromosome:%s:%s:%s:%s:%s" % (reference,self.chrom,chromStart,chromEnd,strand) + seq_id = "%s%s %s %s" % (self.name,frame_name,seqtype,location) + return seq_id + def get_line(self, start_offset = 0, end_offset = 0): + if start_offset or end_offset: + s_offset = start_offset if start_offset else 0 + e_offset = end_offset if end_offset else 0 + if s_offset > self.chromStart: + s_offset = self.chromStart + chrStart = self.chromStart - s_offset + chrEnd = self.chromEnd + e_offset + blkSizes = self.blockSizes + blkSizes[0] += s_offset + blkSizes[-1] += e_offset + blkStarts = self.blockStarts + for i in range(1,self.blockCount): + blkStarts[i] += s_offset + items = [str(x) for x in [self.chrom,chrStart,chrEnd,self.name,self.score,self.strand,self.thickStart,self.thickEnd,self.itemRgb,self.blockCount,','.join([str(x) for x in blkSizes]),','.join([str(x) for x in blkStarts])]] + return '\t'.join(items) + '\n' + return self.line + +def __main__(): + #Parse Command Line + parser = optparse.OptionParser() + parser.add_option( '-i', '--input', dest='input', help='BED file (tophat junctions.bed) with sequence column added' ) + parser.add_option( '-o', '--output', dest='output', help='Translations of spliced sequence') + parser.add_option( '-b', '--bed_format', dest='bed_format', action='store_true', default=False, help='Append translations to bed file instead of fasta' ) + parser.add_option( '-D', '--fa_db', dest='fa_db', default=None, help='Prefix DB identifier for fasta ID line, e.g. generic' ) + parser.add_option( '-s', '--fa_sep', dest='fa_sep', default='|', help='fasta ID separator defaults to pipe char, e.g. generic|ProtID|description' ) + parser.add_option( '-B', '--bed', dest='bed', default=None, help='Output a bed file with added 13th column having translation' ) + parser.add_option( '-G', '--gff3', dest='gff', default=None, help='Output translations to a GFF3 file' ) + parser.add_option( '-S', '--seqtype', dest='seqtype', default='pep:splice', help='SEQTYPE:STATUS for fasta ID line' ) + parser.add_option( '-P', '--id_prefix', dest='id_prefix', default='', help='prefix for the sequence ID' ) + parser.add_option( '-R', '--reference', dest='reference', default=None, help='Genome Reference Name for fasta ID location ' ) + parser.add_option( '-r', '--refsource', dest='refsource', default=None, help='Source for Genome Reference, e.g. Ensembl, UCSC, or NCBI' ) + parser.add_option( '-Q', '--score_name', dest='score_name', default=None, help='include in the fasta ID line score_name:score ' ) + parser.add_option( '-l', '--leading_bp', dest='leading_bp', type='int', default=None, help='leading number of base pairs to ignore when filtering' ) + parser.add_option( '-t', '--trailing_bp', dest='trailing_bp', type='int', default=None, help='trailing number of base pairs to ignore when filtering' ) + parser.add_option( '-U', '--unfiltered', dest='filtering', action='store_false', default=True, help='Do NOT filterout translation with stop codon in the first exon' ) + parser.add_option( '-u', '--untrimmed', dest='untrimmed', action='store_true', default=False, help='Do NOT trim from splice site to stop codon' ) + parser.add_option( '-L', '--min_length', dest='min_length', type='int', default=None, help='Minimun length (to first stop codon)' ) + parser.add_option( '-M', '--max_stop_codons', dest='max_stop_codons', type='int', default=None, help='Filter out translations with more than max_stop_codons' ) + parser.add_option( '-d', '--debug', dest='debug', action='store_true', default=False, help='Turn on wrapper debugging to stdout' ) + (options, args) = parser.parse_args() + # Input files + if options.input != None: + try: + inputPath = os.path.abspath(options.input) + inputFile = open(inputPath, 'r') + except Exception, e: + print >> sys.stderr, "failed: %s" % e + exit(2) + else: + inputFile = sys.stdin + # Output files + bed_fh = None + gff_fh = None + gff_fa_file = None + gff_fa = None + outFile = None + if options.output == None: + #write to stdout + outFile = sys.stdout + if options.gff: + gff_fa_file = tempfile.NamedTemporaryFile(prefix='gff_fasta_',suffix=".fa",dir=os.getcwd()).name + gff_fa = open(gff_fa_file,'w') + else: + try: + outPath = os.path.abspath(options.output) + outFile = open(outPath, 'w') + except Exception, e: + print >> sys.stderr, "failed: %s" % e + exit(3) + if options.gff: + gff_fa_file = outPath + if options.bed: + bed_fh = open(options.bed,'w') + bed_fh.write('track name="%s" description="%s" \n' % ('novel_junctioni_translations','test')) + if options.gff: + gff_fh = open(options.gff,'w') + gff_fh.write("##gff-version 3.2.1\n") + if options.reference: + gff_fh.write("##genome-build %s %s\n" % (options.refsource if options.refsource else 'unknown', options.reference)) + leading_bp = 0 + trailing_bp = 0 + if options.leading_bp: + if options.leading_bp >= 0: + leading_bp = options.leading_bp + else: + print >> sys.stderr, "failed: leading_bp must be positive" + exit(5) + if options.trailing_bp: + if options.trailing_bp >= 0: + trailing_bp = options.trailing_bp + else: + print >> sys.stderr, "failed: trailing_bp must be positive" + exit(5) + # Scan bed file + try: + for i, line in enumerate( inputFile ): + if line.startswith('track'): + if outFile and options.bed_format: + outFile.write(line) + continue + entry = BedEntry(line) + strand = 1 if entry.strand == '+' else -1 + translations = entry.get_translations() + if options.debug: + exon_seqs = entry.get_exon_seqs() + exon_sizes = [len(seq) for seq in exon_seqs] + splice_sites = [sum(exon_sizes[:x]) / 3 for x in range(1,len(exon_sizes))] + print >> sys.stderr, entry.name + print >> sys.stderr, line.rstrip('\r\n') + print >> sys.stderr, "exons: %s" % exon_seqs + print >> sys.stderr, "%s" % splice_sites + for i,translation in enumerate(translations): + print >> sys.stderr, "frame %d: %s" % (i+1,translation) + print >> sys.stderr, "splice: %s" % (''.join(['^' if (((j*3)+i)/3) in splice_sites else '-' for j in range(len(translation))])) + print >> sys.stderr, "" + if options.bed_format: + tx_entry = "%s\t%s\n" % (line.rstrip('\r\n'),'\t'.join(translations)) + outFile.write(tx_entry) + else: + translations = entry.get_filterd_translations(untrimmed=options.untrimmed,filtering=options.filtering,ignore_left_bp=leading_bp,ignore_right_bp=trailing_bp,debug=options.debug) + for i,tx in enumerate(translations): + if tx: + (chromStart,chromEnd,translation,blockCount,blockSizes,blockStarts) = tx + if options.min_length != None and len(translation) < options.min_length: + continue + if options.max_stop_codons != None and translation.count('*') > options.max_stop_codons: + continue + frame_name = '_%s' % (i + 1) + pep_id = "%s%s%s" % (options.id_prefix,entry.name,frame_name) + if bed_fh: + bed_fh.write('%s\t%d\t%d\t%s\t%d\t%s\t%d\t%d\t%s\t%d\t%s\t%s\t%s\n' % (str(entry.chrom),chromStart,chromEnd,pep_id,entry.score,entry.strand,chromStart,chromEnd,entry.itemRgb,blockCount,','.join([str(x) for x in blockSizes]),','.join([str(x) for x in blockStarts]),translation)) + location = "chromosome:%s:%s:%s:%s:%s" % (options.reference,entry.chrom,chromStart,chromEnd,strand) + score = " %s:%s" % (options.score_name,entry.score) if options.score_name else '' + seq_description = "%s %s%s" % (options.seqtype, location, score) + seq_id = "%s " % pep_id + if options.fa_db: + seq_id = "%s%s%s%s" % (options.fa_db,options.fa_sep,pep_id,options.fa_sep) + fa_id = "%s%s" % (seq_id,seq_description) + fa_entry = ">%s\n%s\n" % (fa_id,translation) + outFile.write(fa_entry) + if gff_fh: + if gff_fa: + gff_fa.write(fa_entry) + gff_fh.write("##sequence-region %s %d %d\n" % (entry.chrom,chromStart + 1,chromEnd - 1)) + gff_fh.write("%s\t%s\t%s\t%s\t%s\t%s\t%s\t%d\tID=%s\n" % (entry.chrom,'splice_junc','gene',chromStart + 1,chromEnd - 1,entry.score,entry.strand,0,pep_id)) + for x in range(blockCount): + start = chromStart+blockStarts[x] + 1 + end = start + blockSizes[x] - 1 + phase = (3 - sum(blockSizes[:x]) % 3) % 3 + gff_fh.write("%s\t%s\t%s\t%s\t%s\t%s\t%s\t%d\tParent=%s;ID=%s_%d\n" % (entry.chrom,'splice_junc','CDS',start,end,entry.score,entry.strand,phase,pep_id,pep_id,x)) + """ + ##gff-version 3 + ##sequence-region 19 1 287484 + 19 MassSpec peptide 282299 287484 10.0 - 0 ID=TEARLSFYSGHSSFGMYCMVFLALYVQ + 19 MassSpec CDS 287474 287484 . - 0 Parent=TEARLSFYSGHSSFGMYCMVFLALYVQ;transcript_id=ENST00000269812 + 19 MassSpec CDS 282752 282809 . - 1 Parent=TEARLSFYSGHSSFGMYCMVFLALYVQ;transcript_id=ENST00000269812 + 19 MassSpec CDS 282299 282310 . - 0 Parent=TEARLSFYSGHSSFGMYCMVFLALYVQ;transcript_id=ENST00000269812 + """ + if bed_fh: + bed_fh.close() + if gff_fh: + if gff_fa: + gff_fa.close() + else: + outFile.close() + gff_fa = open(gff_fa_file,'r') + gff_fh.write("##FASTA\n") + for i, line in enumerate(gff_fa): + gff_fh.write(line) + gff_fh.close() + except Exception, e: + print >> sys.stderr, "failed: Error reading %s - %s" % (options.input if options.input else 'stdin',e) + +if __name__ == "__main__" : __main__() +
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/translate_bed_sequences.xml Wed Jan 13 07:38:29 2016 -0500 @@ -0,0 +1,128 @@ +<?xml version="1.0"?> +<tool id="translate_bed_sequences" name="Translate BED Sequences" version="0.1.0"> + <description>3 frame translation of BED augmented with a sequence column</description> + <requirements> + <requirement type="package" version="1.62">biopython</requirement> + <requirement type="python-module">Bio</requirement> + </requirements> + <command interpreter="python"> + translate_bed_sequences.py --input "$input" + #if $fa_db: + --fa_db='$fa_db' + #end if + #if $fa_sep: + --fa_sep='$fa_sep' + #end if + #if $id_prefix: + --id_prefix='$id_prefix' + #end if + #if $reference: + --reference $reference + #else: + --reference ${input.metadata.dbkey} + #end if + #if $refsource: + --refsource $refsource + #end if + #if $seqtype: + --seqtype $seqtype + #end if + #if $score_name: + --score_name $score_name + #end if + #if $filter.filterseqs == 'yes': + #if $filter.leading_bp: + --leading_bp $filter.leading_bp + #end if + #if $filter.trailing_bp: + --trailing_bp $filter.trailing_bp + #end if + #else: + --unfiltered + #end if + #if $trim.trimseqs == 'no': + --untrimmed + #if str($trim.max_stop_codons) != '': + --max_stop_codons $trim.max_stop_codons + #end if + #end if + #if str($min_length) != '': + --min_length $min_length + #end if + --bed $translated_bed + --output "$output" + </command> + <inputs> + <param name="input" type="data" format="bed" label="BED file with added sequence column" + help="Output from 'Extract Genomic DNA' run on tophat junctions.bed "/> + <param name="fa_db" type="text" value="" optional="true" label="fasta ID source, e.g. generic" + help="Any Compomics application such as PeptideShaker, requires a source"> + </param> + <param name="fa_sep" type="text" value="" optional="true" label="fasta ID source, e.g. generic" + help="Only used when a fasta ID source is given, default to the pipe character"> + </param> + <param name="id_prefix" type="text" value="" optional="true" label="ID prefix for generated IDs" + help="Can be used to distinguish samples"> + <validator type="regex" message="Allowed chars:a-z A-Z 0-9 _ - |">^[a-zA-Z0-9_-|]*$</validator> + </param> + <param name="refsource" type="text" value="Ensembl" optional="true" label="Genome reference source" + help=""/> + <param name="reference" type="text" value="" optional="true" label="Genome reference name" + help="By default, the database metadata will be used."/> + <param name="seqtype" type="text" value="" optional="true" label="The SEQTYPE:STATUS to include in the fasta ID lines" + help="For example: pep:splice"/> + <param name="score_name" type="text" value="" optional="true" label="Add the bed score field fasta ID line with this tag name" + help="For example: with the tag name 'depth' and bed score 12: depth:12"/> + <conditional name="filter"> + <param name="filterseqs" type="select" label="Filter out translations with stop codons before the splice site"> + <option value="yes" selected="true">Yes</option> + <option value="no">No</option> + </param> + <when value="yes"> + <param name="leading_bp" type="integer" value="" min="0" optional="true" label="Stop codon filtering start position base pairs" + help="Do not reject translation is stop_codons are within base pairs of the BED start position for positive strand"/> + <param name="trailing_bp" type="integer" value="" min="0" optional="true" label="Stop codon filtering end position base pairs" + help="Do not reject translation is stop_codons are within base pairs of the BED end position for negative strand"/> + </when> + <when value="no"/> + </conditional> + <conditional name="trim"> + <param name="trimseqs" type="select" label="Trim translations to stop codons"> + <option value="yes" selected="true">Yes</option> + <option value="no">No</option> + </param> + <when value="no"> + <param name="max_stop_codons" type="integer" value="" min="0" optional="true" label="Maximum number of stop codons allowed in a translation to be reported"/> + </when> + </conditional> + <param name="min_length" type="integer" value="" min="0" optional="true" label="Minimum length of a translation to be reported"/> + </inputs> + <stdio> + <exit_code range="1:" level="fatal" description="Error" /> + </stdio> + <outputs> + <data name="translated_bed" metadata_source="input" format="bed" label="${tool.name} on ${on_string} bed"> + </data> + <data name="output" metadata_source="input" format="fasta" label="${tool.name} on ${on_string} fasta"> + </data> + </outputs> + <tests> + <test> + <param name="input" value="Extract_Genomic_DNA.bed" ftype="bed" dbkey="hg19"/> + <param name="reference" value="GRCh37"/> + <param name="seqtype" value="pep:novel"/> + <param name="score_name" value="depth"/> + <output name="output" file="translated_bed_sequences.fa"/> + </test> + </tests> + <help> +**Translate BED Sequences** + +This tool takes a BED input file that has been processed +by the Galaxy tool "Extract Genomic DNA" to add a 13th column with the transcript sequence. + +It generates a peptide fasta file with the 3-frame translations of the spliced sequence +defined by each entry in the input BED file. + + </help> +</tool>