Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error with samtools index during TALON run #153

Open
jadedavis5 opened this issue Aug 28, 2024 · 0 comments
Open

Error with samtools index during TALON run #153

jadedavis5 opened this issue Aug 28, 2024 · 0 comments

Comments

@jadedavis5
Copy link

Hi,
I am getting the following error when running TALON using a large plant reference genome:

Command output:
  [ 2024-08-28 02:47:36 ] Started talon_label_reads run.
  [ 2024-08-28 02:47:36 ] Splitting SAM by chromosome...
  [ 2024-08-28 02:47:36 ] -----Sorting and indexing...

Command error:
  [E::hts_idx_check_range] Region 548287349..548289692 cannot be stored in a bai index. Try using a csi index with min_shift = 14, n_lvls >= 6
  [E::sam_index] Read '5edb22b3-d236-4b8e-b51f-d110eca41670' with ref_name='chr2H', ref_length=675310294, flags=16, pos=548287350 cannot be indexed
  Traceback (most recent call last):
    File "/usr/local/bin/talon_label_reads", line 10, in <module>
      sys.exit(main())
    File "/usr/local/lib/python3.7/site-packages/talon/talon_label_reads.py", line 287, in main
      read_files = split_reads_by_chrom(options.sam_file, tmp_dir=options.tmp_dir, n_threads=options.threads)
    File "/usr/local/lib/python3.7/site-packages/talon/talon_label_reads.py", line 177, in split_reads_by_chrom
      pysam.index(bam_file)
    File "/usr/local/lib/python3.7/site-packages/pysam/utils.py", line 75, in __call__
      stderr))
  pysam.utils.SamtoolsError: 'samtools returned with error 1: stdout=, stderr=samtools index: failed to create index for "tmp_label_reads/raw/all_reads.sorted.bam": Numerical result out of range\n'

I am wondering if there is a way to store as a csi index instead?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant