Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • LankyCyril
    Junior Member
    • Aug 2013
    • 1

    #1

    BEAGLE errors: Allele not in marker file / Number of fields in .gprobs isn't %3

    I'm trying to phase a specific genome against the 1000 genomes data.

    What I've done so far:
    1. fastq-dump --split-3: SRA -> FastQ
    2. FastQC
    3. trim_galore
    4. bowtie2: alignment against hg19
    5. samtools, bcftools: SAM/BAM -> VCF
    6. vcftools: VCF -> VCF for chr1 with minQ=50
    7. since the resulting file doesn't have rsIDs in it (the column contains dots), I downloaded snp138.txt.gz from http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database and used the data from it to fill in the missing rsIDs (wrote a script; there's a few rsIDs that are present in the VCF, but missing in snp138, so I ditched the corresponding records)
    8. vcf2beagle.jar: mydata.vcf -> mydata.bgl.gz, mydata.int, mydata.markers
    9. Downloaded the "reference" files for chr1 from this site: http://hgdownload.cse.ucsc.edu/gbdb/hg19/1000Genomes/ (the bigger one, ALL.chr1.integrated_phase1_v3.20101123.snps_indels_svs.genotypes.vcf.gz and its index).
    10. vcf2beagle.jar: reference.vcf -> reference.bgl.gz, reference.int, reference.markers.


    Now I'm trying to run this with
    Code:
    java -jar beagle.jar unphased=mydata.bgl.gz phased=reference.bgl.gz markers=reference.markers missing=? out=phased
    Notice the markers=reference.markers option.
    I also tried it with markers=mydata.markers, but I'm not entirely sure which markers file I should be passing to Beagle (<-- question #1).
    (Also experimented with like=mydata.bgl.gz, didn't affect anything.)

    But there's no difference, because it fails after about 10 seconds with 'ERROR: Allele in data file but not in marker file: allele "C" for marker rs6603854', or a different allele and marker, depending on which marker file I'm using. How do I fix this? (<-- question #2)

    Sometimes (I couldn't replicate this now), it fails with "number of fields in the .gprobs file is not divisible by 3", even though I don't have any .gprobs file to begin with, and it is only generated by Beagle itself. What does this even mean? (<-- question #3)

    Also, am I using the correct files (are they all for hg19?) and tools, even? (<-- question #4)

    Any help will be much appreciated!
    Last edited by LankyCyril; 11-20-2013, 07:54 AM. Reason: fixed formatting

Latest Articles

Collapse

  • SEQadmin2
    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by SEQadmin2




    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
    Today, 06:46 AM
  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-28-2026, 09:51 AM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
67 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
47 views
0 reactions
Last Post SEQadmin2  
Working...