Determination of longest read-SOAPdenovo2

Question

[eischzj12]

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

[GenoMax]

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

[eischzj12]

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?

[GenoMax]

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.