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
      Advancing Precision Medicine for Rare Diseases in Children
      by seqadmin




      Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
      12-16-2024, 07:57 AM
    • seqadmin
      Recent Advances in Sequencing Technologies
      by seqadmin



      Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

      Long-Read Sequencing
      Long-read sequencing has seen remarkable advancements,...
      12-02-2024, 01:49 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 12-17-2024, 10:28 AM
    0 responses
    25 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-13-2024, 08:24 AM
    0 responses
    42 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-12-2024, 07:41 AM
    0 responses
    28 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-11-2024, 07:45 AM
    0 responses
    42 views
    0 likes
    Last Post seqadmin  
    Working...
    X