0
0
Entering edit mode
5.0 years ago
eischzj12 • 0

I'm working on using SOAPdenovo for short-read assembly. How do you determine the length of your longest read? I want to set this equal to the rd_len_cutoff.

Thanks

soapdenovo2 read length Assembly • 914 views
1
Entering edit mode

You can run stats.sh from BBMap on your dataset to get that information. What kind of data is this?

0
Entering edit mode

Paired end Illumina sequence data. Sorry for being vague, I'm not very familiar with the relevant information to provide. Still learning the ropes. BBMap is something I need to download separately then I take it?

0
Entering edit mode

I have linked BBMap in post above.

If your illumina sequencing was done for say 150 cycles (you would have 150 bp reads in that case) then the length of your longest (and perhaps all reads, if they have not been trimmed) would be 150.