schneebergerlab / fixchr Goto Github PK
View Code? Open in Web Editor NEWLicense: MIT License
License: MIT License
Hi,
It says in the readme that a bioconda package is a TODO item. Are you happy to create one? Or I can do it if you like. My only hesitation is that if you want to create a pypi release this will change the bioconda recipe compared to if you don't want to create a pypi package (in which case you would just need to rely on github releases)
Hello again.
Although I've had success using fixchr on some scaffolds, I'm still struggling with the original ones that needed reversecomplementing despite generating pafs and sams with different flags in minimap2. Here's the paf generation with error:
$minimap2 -cx asm20 -L --eqx BAR3230_D13.fasta PSC355_D13.fasta >> PSC355xBAR3230_D13.paf
fixchr:
$fixchr -c PSC355xBAR3230_D13.paf -r BAR3230_D13.fasta -q PSC355_D13.fasta -F P --prefix D13_PSC355xBAR3230
#error
Traceback (most recent call last):
File "/home/test/miniconda3/bin/fixchr", line 4, in <module>
__import__('pkg_resources').run_script('fixchr==0.1', 'fixchr')
File "/home/test/miniconda3/lib/python3.8/site-packages/pkg_resources/__init__.py", line 662, in run_script
self.require(requires)[0].run_script(script_name, ns)
File "/home/test/miniconda3/lib/python3.8/site-packages/pkg_resources/__init__.py", line 1466, in run_script
exec(script_code, namespace, namespace)
File "/home/test/miniconda3/lib/python3.8/site-packages/fixchr-0.1-py3.8.egg/EGG-INFO/scripts/fixchr", line 6, in <module>
File "/home/test/miniconda3/lib/python3.8/site-packages/fixchr-0.1-py3.8.egg/fixchr/scripts/fixchr.py", line 105, in main
File "/home/test/miniconda3/lib/python3.8/site-packages/fixchr-0.1-py3.8.egg/fixchr/scripts/fixchr.py", line 63, in fixchr
File "/home/test/miniconda3/lib/python3.8/site-packages/fixchr-0.1-py3.8.egg/fixchr/scripts/func.py", line 632, in homchr
UnboundLocalError: local variable 'coords2' referenced before assignment
thanks!
Hi,
Thank you for making this tool. I've corrected the contig orientation and generated plots using FicChr; however my contigs are out of order.
This is a species1areference vs species1bquery
All chromosomes are arranged from ch1-ch24 in the original .fasta files, both before and after orientation fixing. Is there a way to correct this?
In the instructions it mentions an .agp file-- I do not have this for my assemblies (is there a tool out there which generates this?), or is there an easier way to correct this?
Hello I'm looking to determine variants between to homologous chromosomes by mapping a whole genome sequence to a reference chromosome via minimap2 and fixchr:
$ minimap2 -x asm20 --eqx Chrom24_CM040141_Aggsing.fa Agt_singleline.mfa >> Chr24_Chr10GST_AggxAgt_20.paf
fixchr -c Chr24_Chr10GST_AggxAgt_20.paf -r Chrom24_CM040141_Aggsing.fa -q Agt_singleline.mfa
but I keep getting an error with the paf, do the files need to be compressed?
Reading Coords - ERROR - astart is not int
Thank you
When I install the software running "python setup.py install", the error occurs:
Traceback`(most recent call last):
File "/data/01/user157/software/fixchr/setup.py", line 18, in <module>
long_description=open('README.rst').read(),
FileNotFoundError: [Errno 2] No such file or directory: 'README.rst'
I can not find the file "README.rst". Could you give me any suggestions?
Dear developers
I ran:
/path/to/apps/fixchr/bin/fixchr \
-c /path/to/some_comparative_genomics/syri_Hp1_to_pm1.delta.m.coords \
-r /path/to/some_comparative_genomics/assembly_pm1_v1_upper.fasta \
-q /path/to/some_comparative_genomics/Hp1.fasta \
--contig_size 5000 \
--dir /path/to/some_comparative_genomics/ \
--prefix Hp1_to_pm1
but the output had no prefix whatsoever.
ref.filtered.fa
qry.filtered.fa
homologous_strand_corrected.pdf
homologous_strand_corrected_alignments.txt
homologous_alignments.txt
homologous.pdf
input.pdf
input_alignments.txt
Which is harmless, unless you run more than one instance.
Kind regards,
AWebb
Hi Manish,
when I run 'fixchr -c S.P.filter.coords -r S.genome.fa -q P.genome.fa -F P',it showed 'Reading BAM/SAM file - ERROR - CIGAR string is not present in PAF at line 2364 8479 37257382 37251209 6116 6174 95.73 1 -1 chr01 chr22. Exiting.' .The S.P.filter.coords file is the output of 'show-coords -THrd' .
I'm wondering which file should be the input(alignment.paf) for fixchr and how to get it ?
Best,
Lily
./minimap2 -cx asm5 asm1.fa asm2.fa > aln.paf
fixchr -c aln.paf -r asm1.fa -q asm2.fa -F P
Hi! @mnshgl0110 .
How to fix it?
Thank you for your reply!
Thanks for the update! I believe I've reinstalled all of the dependencies and pulled the newest version 0.2, but now I'm running into a "hometools" no module issue. I've tried to install this using the standard conda install and python pip to no avail. I see this in the func.py script, line 45 but nowhere else...?
Originally posted by @cistarsa in #3 (comment)
when set "fixchr -F T", What does this Table file refer to?
my command is: "nucmer --mum --mincluster 500 -t 30 $geo1 $geo2 -p XD"
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.