1. Starting with 20 bp reads in fastq format and an ambiguity masked reference
genome execute the following commands on a linux system:
MosaikBuild -fr /mnt/lustre/home/jdegner/OUR_DNASE/AMB_REF/hg_18_ambiguity.gz
-oa /mnt/lustre/home/jdegner/OUR_DNASE/AMB_REF/hg_18_ambiguity.dat
MosaikBuild -q
/mnt/lustre/home/jdegner/OUR_DNASE/MAPPABILITY_TESTS/SIMULATED_SNPS_CHR22/chr22_
SNPs.fastq -out
/mnt/lustre/home/jdegner/OUR_DNASE/MAPPABILITY_TESTS/SIMULATED_SNPS_CHR22/chr22_
SNPs.20bp.dat -st illumina
MosaikJump -ia /mnt/lustre/home/jdegner/OUR_DNASE/AMB_REF/hg_18_ambiguity.dat
-out
/mnt/lustre/home/jdegner/OUR_DNASE/AMB_REF/hg_18_ambiguity.jump.dat -hs 15
MosaikAligner -in
/mnt/lustre/home/jdegner/OUR_DNASE/MAPPABILITY_TESTS/SIMULATED_SNPS_CHR22/chr22_
SNPs.20bp.dat -out
/mnt/lustre/home/jdegner/OUR_DNASE/MAPPABILITY_TESTS/SIMULATED_SNPS_CHR22/chr22_
SNPs.20bp.alligned.dat -ia
/mnt/lustre/home/jdegner/OUR_DNASE/AMB_REF/hg_18_ambiguity.dat -j
/mnt/lustre/home/jdegner/OUR_DNASE/AMB_REF/hg_18_ambiguity.jump.dat -hs 15 -m
unique -mm 1 -mhp 100 -act 20 -p 5
While I expect this to produce an output with aligned reads, I receive this
error:
- loading reference sequence... ERROR: Unable to initialize the uncompressed
FastLZ buffer (Read).
Any suggestions for where to look for a solution would be most appreciated
Jack
Original issue reported on code.google.com by JFDeg...@gmail.com on 28 Apr 2010 at 6:30
Original issue reported on code.google.com by
JFDeg...@gmail.com
on 28 Apr 2010 at 6:30