You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 2 Current »

I am trying to simulate a full-scale experiment of ~30M, 100-nt reads for Arabidopsis, but I keep getting the same error (also at 16M reads, about to try 8M reads):

[INFO] Loading default PCR distribution
    preparing transcript sequences *******Problems reading 3: 23459833, 74> 23459834 into 100: null
check for the right species/genome version!

I looked at the GTF and FASTA files.  The maximum position for chromosome 3 in the GTF file is 23,459,804; the actual chromosome sequence length is 23,459,831.  For some reason the simulator is trying to go beyond the end of the last gene in chromosome 3?  The full command I'm using is:

flux-simulator -t simulator -x -l -s -p a_thaliana_flux.par

I have attached the .par file I'm using for 8M reads in case this will help: a_thaliana_flux.par

Thanks!

  • No labels