Announcement

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

  • BFAST bfast.submit.pl configuration

    Could anyone please provide a working config file for bfast.submit.pl to use as an example? We tried to create a config file, but it's sometimes hard to guess which parameter name in the schema stands for which in the BFAST programs. Also, Eclipse reports that the provided xml schema is incorrect at some places. Despite our efforts, the bfast.submit.pl exited without producing output. (We're using the data from the latest BFAST version.) bfast.submit.pl looks like a valuable tool and it would be very useful to get it to run for our cluster.

    In that context, I'd like to know what the most efficient way of running BFAST is. I can use a node with 16 CPUs, up to 128 GB RAM. The 10 indexes for the human genome are 12 GB each so it's probably impossible to load them all into memory and keep enough space for the rest, especially when using pipes. As I noted, reading the indexes (one at a time as done by default) is the most time-consuming part in our case. Instead of splitting up the reads much and call multiple instances of bfast match with all indexes, I think it would be better to process all reads with one of the indexes in parallel.

    Thanks in advance for the help

    Barbara

  • #2
    Dears,
    this is my config file.
    I have not a sam file in the output, but only some sh script: match.sh, localign.sh, postprocess.sh, sam.sh.
    I don't know what is the next step because, as i write you, it's sometimes hard to guess which parameter name in the schema stands for which in the BFAST programs.
    A suggestion:
    if you have any validation errors with eclipse don't worry, they're not preclude the running of bfast.
    I am a beginner so i don't know if this can help you. I hope..
    Anyway, this is the BASIC xml version to run it:
    <?xml version="1.0" encoding="UTF-8"?>
    <bfastConfig xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="BfastConfig_5.xsd">
    <globalOptions>
    <bfastBin>/cbbl/users/vincenzo/bin/</bfastBin>
    <samtoolsBin>/cbbl/users/vincenzo/bin/</samtoolsBin>
    <qsubBin>/cbbl/users/vincenzo/bin/run_parallel/</qsubBin>
    <fastaFileName>/cbbl/users/vincenzo/bin/run_parallel/DH10B.fa</fastaFileName>
    <runDirectory>/cbbl/users/vincenzo/bin/run_parallel/</runDirectory>
    <readsDirectory>/cbbl/users/vincenzo/bin/run_parallel/</readsDirectory>
    <outputDirectory>/cbbl/users/vincenzo/bin/run_parallel/</outputDirectory>
    <tmpDirectory>/cbbl/users/vincenzo/bin/run_parallel/</tmpDirectory>
    <outputID>result</outputID>
    <cleanUsedIntermediateFiles>0</cleanUsedIntermediateFiles>
    <numReadsPerFASTQ localalignSplit="20" matchSplit="20">20</numReadsPerFASTQ>
    <timing>ON</timing>
    <queueType>SGE</queueType>
    <space>NT</space>
    </globalOptions>
    <matchOptions>
    <readCompression>gz</readCompression>
    <threads>1</threads>
    </matchOptions>
    <localalignOptions>
    <threads>1</threads>
    </localalignOptions>
    <postprocessOptions>
    <threads>1</threads>
    </postprocessOptions>
    <samOptions>
    <samtools>1</samtools>
    </samOptions>
    </bfastConfig>

    To add any parameter, you can inspire from XSD and from bfast.submit.pl.
    If you solved it, please share your code.
    Thanks a lot.
    Vincenzo

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Advanced Tools Transforming the Field of Cytogenomics
      by seqadmin


      At the intersection of cytogenetics and genomics lies the exciting field of cytogenomics. It focuses on studying chromosomes at a molecular scale, involving techniques that analyze either the whole genome or particular DNA sequences to examine variations in structure and behavior at the chromosomal or subchromosomal level. By integrating cytogenetic techniques with genomic analysis, researchers can effectively investigate chromosomal abnormalities related to diseases, particularly...
      09-26-2023, 06:26 AM
    • seqadmin
      How RNA-Seq is Transforming Cancer Studies
      by seqadmin



      Cancer research has been transformed through numerous molecular techniques, with RNA sequencing (RNA-seq) playing a crucial role in understanding the complexity of the disease. Maša Ivin, Ph.D., Scientific Writer at Lexogen, and Yvonne Goepel Ph.D., Product Manager at Lexogen, remarked that “The high-throughput nature of RNA-seq allows for rapid profiling and deep exploration of the transcriptome.” They emphasized its indispensable role in cancer research, aiding in biomarker...
      09-07-2023, 11:15 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 09:38 AM
    0 responses
    9 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 09-27-2023, 06:57 AM
    0 responses
    11 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 09-26-2023, 07:53 AM
    1 response
    23 views
    0 likes
    Last Post seed_phrase_metal_storage  
    Started by seqadmin, 09-25-2023, 07:42 AM
    0 responses
    17 views
    0 likes
    Last Post seqadmin  
    Working...
    X