Header Leaderboard Ad

Collapse

Solid-color space indexing for human fasta file

Collapse

Announcement

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

  • Solid-color space indexing for human fasta file

    Hello,
    I am using BFAST for indexing color space data using the command:

    bfast index -f folder/human_bfast.fa -m '11111' -w 5 -A 1 -n 8

    If I try to index the fasta file with a few chromosomes, say till 17, the tool successfully creates .bif file. But if I use a bigger fasta file (ie. from chr 1 till 22) for indexing then the tool throws error:

    Currently on [contig,pos]:
    [------22,------662000]
    Sorting by thread...
    100.000 percent complete************************************************************
    In function "RGIndexMergeHelperFromDiskContig_8": Fatal Error[ReadFileError]. Message: Could not read in tmp lower.
    The file stream error was:: Bad file descriptor
    ***** Exiting due to errors *****
    ************************************************************
    I want to create solid color space index for human fasta for all chromosomes from 1 through 22, X,Y

    Please note that the step just prior to indexing (ie. creating reference in color space using the same fasta) was always successful. It creates .fa.cs.brg file successfully.
    Can anybody please help me with the problem ?

    Thanks in advance,
    Aswathy

  • #2
    Try with a longer index, I think the recommended were between 12-14 long? 5 bases with no open positions does not make much sense for human. Masking strategies are described in the publications/documentations.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Improved Targeted Sequencing: A Comprehensive Guide to Amplicon Sequencing
      by seqadmin



      Amplicon sequencing is a targeted approach that allows researchers to investigate specific regions of the genome. This technique is routinely used in applications such as variant identification, clinical research, and infectious disease surveillance. The amplicon sequencing process begins by designing primers that flank the regions of interest. The DNA sequences are then amplified through PCR (typically multiplex PCR) to produce amplicons complementary to the targets. RNA targets...
      03-21-2023, 01:49 PM
    • seqadmin
      Targeted Sequencing: Choosing Between Hybridization Capture and Amplicon Sequencing
      by seqadmin




      Targeted sequencing is an effective way to sequence and analyze specific genomic regions of interest. This method enables researchers to focus their efforts on their desired targets, as opposed to other methods like whole genome sequencing that involve the sequencing of total DNA. Utilizing targeted sequencing is an attractive option for many researchers because it is often faster, more cost-effective, and only generates applicable data. While there are many approaches...
      03-10-2023, 05:31 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-24-2023, 02:45 PM
    0 responses
    16 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2023, 12:26 PM
    0 responses
    17 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-17-2023, 12:32 PM
    0 responses
    17 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-15-2023, 12:42 PM
    0 responses
    24 views
    0 likes
    Last Post seqadmin  
    Working...
    X