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

Incomplete prunning #175

Open
Alteroldis opened this issue Jan 4, 2024 · 1 comment
Open

Incomplete prunning #175

Alteroldis opened this issue Jan 4, 2024 · 1 comment

Comments

@Alteroldis
Copy link

Hi Dr. Zhang,
At first, you create elegant tool. But with my species I be faced with incomplete prunning. My allele.ctg.table include many contigs which incorporated in same group and sequence. The contigs not included at removedb_Allele.txt, prunning.pairs.txt and at the other side, included at removedb_nonBest.txt. I used the command:

ALLHiC_prune -i Allele.ctg.tsv -b $finalbam.bam -r $assembly
ALLHiC_partition -b prunning.bam -r $assembly -e GATC -k 46
ALLHiC_rescue -b $finalbam.bam -r $assembly -c prunning.clusters.txt -i prunning.counts_GATC.txt
allhic extract $finalbam.bam $assembly --RE GATC
allhic optimize group${K}.txt $finalbam.clm
ALLHiC_build $assembly

v.0.9.8 from your repo.
P.S. About half of genes which used for creating allele.table were placed on the same chromosome (group). And half of contig pairs (alleles, haplotigs) from allele.table placed in one group at prunning.clusters.txt file.
P.P.S. "allhic alleles" cmd does not create any tables, only stdout like this:

02:23:31 ParseRecords | NOTICE Parse paffile genome.paf
02:29:15 ReadCSVLines | NOTICE Parse csvfile ../bwa-curr/filt.counts_GATC.txt
02:29:15 Run | NOTICE Success

@13061994
Copy link

13061994 commented Jan 8, 2024

Could you explain how to construct allele.ctg.table?

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

2 participants