Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BFAST index, not understand warnings

    Hello,

    I want to created index files in color space from NCBI build 36 in color space.
    Here is my command
    `which time` bfast index -f NCBI36.fa -n 8 -w 14 -i 1 -A 1 -m 1111111111111111111111

    Then I got warning like "startContig was less than zero.", "endContig was greater than the number of contigs in the reference genome.". But I am sure I am using the same FASTA reference genome file. So is there anything I should do for these warnings? And is anyone have experience about how long it takes to build this index file?

    Part of the outputs:
    ------------------------------------------------------
    Checking input parameters supplied by the user ...
    Validating fastaFileName NCBI36.fa.
    Validating tmpDir path ./.
    Input arguments look good!
    ************************************************************
    ************************************************************
    Printing Program Parameters:
    programMode: [ExecuteProgram]
    fastaFileName: NCBI36.fa
    space: [Color Space]
    mask: 1111111111111111111111
    depth: 0
    hashWidth: 14
    indexNumber: 1
    repeatMasker: [Not Using]
    startContig: 0
    startPos: 0
    endContig: 2147483647
    endPos: 2147483647
    exonsFileName: [Not Using]
    numThreads: 8
    tmpDir: ./
    timing: [Not Using]
    ************************************************************
    ************************************************************
    Reading in reference genome from NCBI36.fa.cs.brg.
    In total read 23 contigs for a total of 3022646526 bases
    ************************************************************
    Creating the index...
    ************************************************************
    Warning: startContig was less than zero.
    Defaulting to contig=1 and position=1.
    ************************************************************
    ************************************************************
    Warning: endContig was greater than the number of contigs in the reference genome.
    Defaulting to reference genome's end contig=23 and position=154913754.
    ************************************************************
    Currently on [contig,pos]:
    [-------4,---27000000]

    ------------------------------------------------------

  • #2
    Originally posted by zhanxw View Post
    Hello,

    I want to created index files in color space from NCBI build 36 in color space.
    Here is my command
    `which time` bfast index -f NCBI36.fa -n 8 -w 14 -i 1 -A 1 -m 1111111111111111111111

    Then I got warning like "startContig was less than zero.", "endContig was greater than the number of contigs in the reference genome.". But I am sure I am using the same FASTA reference genome file. So is there anything I should do for these warnings? And is anyone have experience about how long it takes to build this index file?

    Part of the outputs:
    ------------------------------------------------------
    Checking input parameters supplied by the user ...
    Validating fastaFileName NCBI36.fa.
    Validating tmpDir path ./.
    Input arguments look good!
    ************************************************************
    ************************************************************
    Printing Program Parameters:
    programMode: [ExecuteProgram]
    fastaFileName: NCBI36.fa
    space: [Color Space]
    mask: 1111111111111111111111
    depth: 0
    hashWidth: 14
    indexNumber: 1
    repeatMasker: [Not Using]
    startContig: 0
    startPos: 0
    endContig: 2147483647
    endPos: 2147483647
    exonsFileName: [Not Using]
    numThreads: 8
    tmpDir: ./
    timing: [Not Using]
    ************************************************************
    ************************************************************
    Reading in reference genome from NCBI36.fa.cs.brg.
    In total read 23 contigs for a total of 3022646526 bases
    ************************************************************
    Creating the index...
    ************************************************************
    Warning: startContig was less than zero.
    Defaulting to contig=1 and position=1.
    ************************************************************
    ************************************************************
    Warning: endContig was greater than the number of contigs in the reference genome.
    Defaulting to reference genome's end contig=23 and position=154913754.
    ************************************************************
    Currently on [contig,pos]:
    [-------4,---27000000]

    ------------------------------------------------------
    You can ignore the warnings if you didn't specify the "-s" or "-e" options. Each index takes about 4-5 hours to build (8 threads).

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Non-Coding RNA Research and Technologies
      by seqadmin




      Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

      Nobel Prize for MicroRNA Discovery
      This week,...
      10-07-2024, 08:07 AM
    • seqadmin
      Recent Developments in Metagenomics
      by seqadmin





      Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
      09-23-2024, 06:35 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 10-02-2024, 04:51 AM
    0 responses
    103 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-01-2024, 07:10 AM
    0 responses
    112 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 09-30-2024, 08:33 AM
    1 response
    115 views
    0 likes
    Last Post EmiTom
    by EmiTom
     
    Started by seqadmin, 09-26-2024, 12:57 PM
    0 responses
    21 views
    0 likes
    Last Post seqadmin  
    Working...
    X