Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • colindaven
    replied
    From memory I think readseq.jar might be able to help, it is simple too.

    Leave a comment:


  • pmiguel
    replied
    The "bp_" prefix stands for "BioPerl". You might want to ask the BioPerl programmers what the problem might be. They have an email list. Possibly there are other ways to query them as well.

    --
    Phillip

    Leave a comment:


  • Having trouble converting genbank format to gff format

    I have a fairly large number of sequin annotation files generated by GenBank's automated PGAP pipeline that I want to convert into bed format. I've been able to use the asn2gb program from the ncbi toolkit to generate GenBank format files from those sequin files, but so far I haven't had much luck converting from genbank into gff. After getting into gff, conversion into bed won't be a problem. But so far I haven't had any luck going from genbank to gff.

    I am trying to use bp_genbank2gff3.pl, but am getting this error message:

    --------------------- WARNING ---------------------
    MSG: Bad LOCUS name? Changing [Staphylococcus_spHMPREF3292-1.0_Cont0>11789] to 'unknown' and length to Staphylococcus_spHMPREF3292-1.0_Cont0>11789
    ---------------------------------------------------

    I thought maybe it was unhappy because of the greater-than symbol in the LOCUS line. So I tried removing it, and then I get a new error message:

    ------------- EXCEPTION: Bio::Root::Exception -------------
    MSG: asking for tag value that does not exist date
    STACK: Error::throw
    STACK: Bio::Root::Root::throw /gsc/scripts/opt/genome/current/user/lib/perl/Bio/Root/Root.pm:357
    STACK: Bio::SeqFeature::Generic::get_tag_values /gsc/scripts/opt/genome/current/user/lib/perl/Bio/SeqFeature/Generic.pm:498
    STACK: main::gff_header /gsc/bin/bp_genbank2gff3.pl:895
    STACK: /gsc/bin/bp_genbank2gff3.pl:406
    -----------------------------------------------------------


    All my input files were generated by GenBank itself, and I used a GenBank script to generate the .gbk files I now have. So I would think my .gbk files are correctly formatted. Can anyone either advise me on what I'm doing wrong with the 'bp_genbank2gff3.pl' script? Or if there is a better way to convert from .gbk into .gff that would work too.


    Thanks,
    John Martin

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 06:54 AM
0 responses
10 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-14-2024, 07:24 AM
0 responses
17 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-13-2024, 08:58 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-12-2024, 02:20 PM
0 responses
17 views
0 likes
Last Post seqadmin  
Working...
X