Hi all, I have some results retrieved from ProDom database and I need to analyze the output from 3 inputs to check if they share the same conserved domains, the results points to same description of 'ribonucleoprotein', then intuitively we can conclude that they share same domains, but I need to automatize this process in a script writing some algorithm.
The inputs are:
#QUERY:
>Contig4 60S ribosomal protein L6
CGGCACGAGGCCCAGACGACGACACCAAAACGAAACACCCACCCGCTCCCTCCGCTGAAA
ACTTGGACCTCGCACTCCTACTTTTCTTCTTTTCCACATCTTCGTCACAATAGCCAATAT
GTCGGCGACTTTGGACTCAAAGTCTTTTGGCCAGACGAAGAAGTTTGGAAAGGGCGAGAG
GACCATCCCCAGCCAAAAGGCTTCAAAGTGGTATCCCACTGAGGACGAGCCACAGCCAAA
GAAAGTCCGCAAGACTATCCACCCAGCGAAGCCCCGGGCCAGTCTTCAACCGGGCACCAT
CCTCATCCTCCTCGCTGGTCGCTTCCGAGGCAAACGTGTTGTGCTCCTTAAGCACCTCCA
CCAAGGTGTTCTCCTCGTTACCGGTCCTTTCAAACTCAATGGCGTCCCTCTTCGGAGAGT
AAACGCCAGATACGTCATTGCCACCAGCTCGAAGGTGGACCTCACGGGAATTGACGACAA
GGTTCTGGAGAAAGCCTCGGAATCCGAGTACTTCACTCGTGAGAAGAAGGCCGAGAAGAA
GGGAGAGGAAGCTTTCTTCAAACAAGGAGAGAAGCCAGAGAAGAAGAAGGTCGTCAGCGC
CCGTGCCAACGACCAAAAGGCCATAGATCGGCCATTGTTGGCCACCATTAAGAAGGAGCA
GTTCCTCGCCAGCTACCTCTCCACCAGCTTCAGCCTCCGGAAAGGCGACAAGCCTCATGA
AATGAAGTGGTAAAACCGGCGGGTTGTCTACGTTGATGGCATAGAAAGAATGGGGATGTT
TTTTCCTTATTTTCAATTGTTTTGTTTCTTGGACGGGTCCAACGGGGGAAATTTTCTTTT
CAGGAAATAAATGGAAAAAAACATTAAAATATACAAACACTTCCAAAGGTCCTTTTTCAA
TGAAATAGTGTACATGGATTGGCCGATATTTCCCCTCGACTGATAATTTCAACGAAAACC
CTTCTAAATTCCCCACCTTTTCCGCAGG
NR-fungi:
↧
