Code Monkey home page Code Monkey logo

Comments (5)

ziyewang avatar ziyewang commented on July 30, 2024

Hi, Kika,

It seems that the kmer_profile file contains contigs longer than 500bp, but the coverage profile file contains all the contigs. Please remove the contigs no longer than 500bp for the coverage profile file and the contig file as README, and give it another try.

"Remove the contigs no longer than 500bp like this for the maxbin depth:

cat mb2_master_depth.txt | awk '{if ($2>500) print $0 }' | cut -f -1,4- > coverage_profile_f500.tsv

Remove the contigs no longer than 500bp like this for contig file:
cd path_to_MetaBinner
cd scripts

python Filter_tooshort.py final.contigs.fa 500
"

Please feel free to let us know if you have any further questions.

Best,
Ziye

from metabinner.

kikinocka avatar kikinocka commented on July 30, 2024

Hi Ziye,
I ran into the same problem even after replacing coverage and contig files. After examining number of lines in coverage and composite files, I noticed a different number of records. All contigs with length 500nt were missing in the composition file. After setting minimal length parameter in gen_kmer.py to 499, run_metabinner.sh command ran properly. So I think the minimal length parameter (length_threshold) reports contigs 1nt longer than the set value.

Kika

from metabinner.

ziyewang avatar ziyewang commented on July 30, 2024

Hi, Kika,

Good to know it is okay to run and thanks for your information. The program keeps the contigs longer than the "length_threshold". It is kind of confusing and we will fix the problem in the next release.

Best,
Ziye

from metabinner.

maggieZhangmj avatar maggieZhangmj commented on July 30, 2024

Hi, I wanted to run MetaBinner on my metagenome assembly. I have already produced coverage and composite profile files, but when I loaded them to run_metabinner.sh, the run ended with the error in component_binning.py script, see here.

Can you please take a look?

Thanks, Kika

Hi Kika!
I also ran into the same problems with you, I tried to change the "length_threshold", from length_threshold = int(sys.argv[2]) to length_threshold = int(sys.argv[2])-1, ending up with no change. I still cannot solve this problem, can you show me what specific details you changed in "gen_kmer.py"?

from metabinner.

kikinocka avatar kikinocka commented on July 30, 2024

Hi Ziye,
my gen_kmer.py command looked like: python gen_kmer.py contigs.fasta 499 4. The first number (499) is length_threshold parameter, the second (4) is the kmer_length parameter. So, I changed this parameter only in my command, not in the python script itself.

Hope this helps!
Kika

from metabinner.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.