Git Product home page Git Product logo

scripts's Introduction

Michigan Geomicrobiology Lab

Welcome to the GitHub Repo for some general purpose NGS, Data analysis and mining scripts used in the lab. Some scripts implement our short-read QC, assembly, binning, etc. pipeline and depend on the presence of a number of third-party software. The rest are Bash scripts or in core Perl or Python. This means, if you have Perl or Python 3 installed, you won't need anything else to work with these scripts.

Since these scripts are actively being used by the Lab, you can expect full support for any issues. Please do let us know if you find any bugs or easier/quicker/more elegant solutions.

Language and OS Dependencies

The scripts should work with variuos flavors of Linux and other unix-like environments. Here is a list of easy to install languages that you'll need:

  • Perl version 5.10 +
  • Python version 3.5 +
  • R version 3 +

Contact

Please send questions or comments to [email protected].

Principal Investigator

Gregory J. Dick, gdick [AT] umich [DOT] edu

License

Geo-omics-scripts is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

Disclaimer

Geo-omics scripts are distributed in the hope that they will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

scripts's People

Contributors

robert102 avatar sunitj avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

scripts's Issues

esomTrain.pl does not generate a wts file

Hello,

I ran your esomTrain.pl script while following the tetramerFreqs/Binning tutorial but the output is not a file with the extension .wts and hence, I'm not able to load the file in the ESOM gui as a previous work file. Any suggestions on how else I can train the data on the command line please?

Thanks

gff2tbl.pl problem

I want to get the NCBI tbl format using gff2tbl.pl, can you give me an example for JGI's GFF and the gene product file from JGI? Thank you.

gff2tbl not providing locus_tag for CRISPR

The gff file (and the IMG genome from which it is derived) does not have a locus tag for CRISPR:

Ga0115370_1019  CRT repeat_region   34874   36813   .   1   .   ID=Ga0115370_1019.29;Average_repeat_length=37;Number_of_repeats=27;rpt_type=CRISPR;rpt_unit=34948..34985;
...
Ga0115370_1111  Prodigal V2.6.3 February, 2016  CDS 64316   64669   .   1   0   ID=Ga0115370_1111.63;conf=98.59;gc_cont=0.486;locus_tag=Ga0115370_111160;
...
Ga0115370_1111  CRT repeat_region   80076   80257   .   1   .   ID=Ga0115370_1111.79;Average_repeat_length=35;Number_of_repeats=3;rpt_type=CRISPR;rpt_unit=80149..80184;

So gff2tbl does not provide a value for locus_tag for CRISPR hypothetical protein:

>Feature Ga0115370_1019
...
31990   34623   CDS
            locus_tag   Ga0115370_101928
            product  CRISPR-associated helicase, Cas3 family
...
>Feature Ga0115370_1111
80076   80257   CRISPR
            locus_tag   
            note    hypothetical protein
            note    locus=''
64316   64669   CDS
            locus_tag   Ga0115370_111160
            product  translation initiation factor 1 (eIF-1/SUI1)
...

This makes it problematic for tbl2asn:

[tbl2asn 24.9] ERROR: Qualifier 'locus_tag' has no value on misc_feature feature 
at lcl|Ga0115370_1111:80076-80257, relative line 5693
[tbl2asn 24.9] ERROR: Unknown feature CRISPR

IMG does not have locus labels for CRISPR regions. Perhaps user has to manually label CRISPR regions with respect to their locations in the scaffold, e.g. Ga0115370_101928.5 and Ga0115370_111175.5, in order to satisfy tbl2asn? Or they are not parsed in gff2tbl, or unique locus_tags have to be generated on the fly in gff2tbl.

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.