Question: How do I carry out small RNA assembly?
gravatar for jaqx008
13 months ago by
jaqx00870 wrote:

Hello everyone. I have a small RNA library that I am trying to see if the target any kind of viruses. I plan to make these short sequences into longer reads or contigs and then see if they map to any viral genomes. I was advised to use velvet to carry out the assembly. I installed velvet using conda install velvet. to do the assembly I used the command I found in the velvet help according to how I understand it.

velveth output.fastq 191 -fastq sample.fastq

the problem is, I am not sure this is right because the output file were still short reads when I looked and looked as below


>NS500519:44:HHHTLBGX2:3:11406:4633:13982   21398134    0
>NS500519:44:HHHTLBGX2:3:11406:15103:13982  21398135    0
>NS500519:44:HHHTLBGX2:3:11406:23757:13983  21398136    0
>NS500519:44:HHHTLBGX2:3:11406:9883:13983   21398137    0
ADD COMMENTlink modified 13 months ago • written 13 months ago by jaqx00870

smallRNA's are of a size that should not need assembly. If you need to see if they map to viral genomes, you could do that with the data you have now. If the data originally did not come from an entity that was long to begin with there is no point in trying to assemble.

Instead you could look into from BBMap suite that can be used for read extension/error correction as an alternative.

ADD REPLYlink written 13 months ago by genomax90k

I apologize for late response. The data I have are small RNAs mostly piRNA and siRNAs that can likely target viral elements (in invertebrates). making the reads longer can help me blast my data against viral databases to see if the small RNAs originate from any viral genomes as I dont have a specific virus in mind. However, I will see if tadpole can help with the read extension. Thanks again

ADD REPLYlink written 13 months ago by jaqx00870
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 887 users visited in the last hour