Thanks for looking at this Liubov. Are you actively using ATAC?
To be honest, I don't know how to read the chain file format. It's possible it's incorrectly representing our 1-based database data.
Hi Liguo, Thanks for getting back to me . The chain file I'm using is here: ftp://ftp.ensemblgenomes.org/pub/current/plants/assembly_chain/sorghum_bicolor/Sorghum_bicolor_v2_to_Sorghum_bicolor_NCBIv3.chain.gz Our internal (1-based) database has the following mapping: V2: 64231291 to 64603896 -> V3: 71417888 to 71790493 When I use our assembly converter tool (backed by crossmap) to 'see' this mapping using BED format, inputting this works: 1 64231291 64603896 but if I 'properly' 0-base the bed file:...
To be honest, I don't know how to read the chain file format. It's possible it's incorrectly representing our 1-based data.
Hi Liguo, Thanks for getting back to me . The chain file I'm using is here: ftp://ftp.ensemblgenomes.org/pub/current/plants/assembly_chain/sorghum_bicolor/Sorghum_bicolor_v2_to_Sorghum_bicolor_NCBIv3.chain.gz Our internal (1-based) database has the following mapping: V2: 64231291 to 64603896 -> V3: 71417888 to 71790493 When I use our assembly converter tool (backed by crossmap) to 'see' this mapping using BED format, inputtign this works: 1 64231291 64603896 but if I 'properly' 0-base the bed file:...
BED off by one projection error.
Mapping VCF over an inversion doesn't revcomp the alleles
matchExtender: Assertion `0' failed
A mere 18Tb free on that device, so I don't think that's the issue. Unfortunately,...
Possible. I'll check.
Uh?
What's weird is that I've seen the -crossspecies mode run in other cases.
What does seatac do?
Is this still an issue? Can you explain the output of leaff --partitionmap and how...
Error running ATAC in --crossspecies mode
Putting the sequence 'all on one line' in the input fasta fixed the issue we were...
Yeah, the 1-to-1 requirement is 'breaking' some of our alignments, so I'm trying...
Ah, I remember this now. I thought the pipeline did this itself, but obviously not....
Sorry, can't seem to edit to fix formatting...
KeyError in DNA.DNA.reversecomplement()
Thanks Brian, by space-based, do you mean zero offset? We aren't seeing this uniquely...
From which tool?
What should they do?
Can you add more details here?
Did this patch get merged?
Can you submit a patch to fix this?
Can I close this now?
Alignments sometimes off by one?