This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Determination of longest read-SOAPdenovo2

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

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

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?

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.

0 answers

No answers yet.

Log in to answer this question.