I have some data of Arabidopsis for small RNA analysis. Now i want to do a locall alignment with Rfam then exclude the sequences which match to Rfam from futher analysis. But the question is:
Do i need to extract the relevant Arabidopsis members of the Rfam before alignning? If it's best to do that, how should i deal with members with no species name? Do they also need to extract for matching? These members are as follows:
>RF00005;tRNA;AAVX01183092.1/735-806 :
GTCCCATAGCATTGTGGTTTGCGTATTCGCCTTGCAAGCGAGAGGATCTG
GGTTCGATTCCCAGGTGGGGCA
>RF00005;tRNA;AAVX01300616.1/338-420 :
GTAGTTGTGGCTGAGTGGTTAAGACAATGGACTAGAAATCCACTAGAAAA
CTACCCGCGCAGGTTCGAATCCTGCAGACTACG
>RF00005;tRNA;AAVX01004056.1/1665-1747 :
GAAGTATTGGCCGAGTGGCGAAGGCGATGGACTACAAATCCATTGAGAAA
CTTCCCGCGCTGTTTCGAATCCTGCCGACTACG
Looking forward to your reply~~
Thank you in advance!
Do i need to extract the relevant Arabidopsis members of the Rfam before alignning? If it's best to do that, how should i deal with members with no species name? Do they also need to extract for matching? These members are as follows:
>RF00005;tRNA;AAVX01183092.1/735-806 :
GTCCCATAGCATTGTGGTTTGCGTATTCGCCTTGCAAGCGAGAGGATCTG
GGTTCGATTCCCAGGTGGGGCA
>RF00005;tRNA;AAVX01300616.1/338-420 :
GTAGTTGTGGCTGAGTGGTTAAGACAATGGACTAGAAATCCACTAGAAAA
CTACCCGCGCAGGTTCGAATCCTGCAGACTACG
>RF00005;tRNA;AAVX01004056.1/1665-1747 :
GAAGTATTGGCCGAGTGGCGAAGGCGATGGACTACAAATCCATTGAGAAA
CTTCCCGCGCTGTTTCGAATCCTGCCGACTACG
Looking forward to your reply~~
Thank you in advance!
Comment