Biostar Beta. Not for public use.
SOAPdenovo trans min contig issue
0
Entering edit mode
14 months ago
Biogeek • 350

Hey guys,

I'm running soapdenovo TRANS and I am trying to give a cut-off contig length of minimum 200, rather than the 100 default using the -L flag.

./SOAPdenovo-Trans-127mer all -s config_file -o K91 -K 91 -L 200 -p 30

My output file .scafStatistics, still reports (sequences/scaffolds above 100bp). Why is the -L flag not working, should I be using some other flag?

Additionally, is it normal to keep gaps put out in the scafSeq file or do people close them in transcriptome assemblies?

Thanks.

ADD COMMENTlink
0
Entering edit mode

I have faced with the same problem about the min contig length. After running, you can remove all contigs less than 200 bp from the final assembly. Regarding gap, I used a software (sorry, don't remember its name), but its results weren't fantastic.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1