Comments (10)
You can also run Gene2 through pxlssq
; if it fails there, it is almost certainly the line breaks:
pxlssq -s gene2_trimAl.txt
from phyx.
from phyx.
List.txt
[gene1_trimAl.txt](https://github.com/FePhyFoFum/phyx/files/7223294/ge
gene2_trimAl.txt
ne1_trimAl.txt)
from phyx.
I do not understand the problem. Why were the results incorrect? What were the results? What are the expected results? What was the command you used? Also, the upload of gene1_trimAl.txt failed.
from phyx.
gene1_trimAl.txt.txt
gene2_trimAl.txt.txt
name_List.txt
I'm sorry I wasn't clear enough.Gene1 can removes unwanted sequences normally, but other genes, such as Gene2, failed. I dont Know the reason?
from phyx.
I'm sorry but I still do not understand what went wrong.
What is wrong with the gene2 results? What were they? What were you expecting? What command did you use?
from phyx.
This is what I get when I examine your files:
Gene1 (original):
josephwb@Potiphar-Breen:issue$ pxlssq -s gene1_trimAl.txt
File type: fasta
Number of sequences: 84
Is aligned: true
Sequence length: 1086
--------- Nucl TABLE ----------
Nucl Total Proportion
A 19585 0.214691
C 24941 0.273404
G 22558 0.247281
T 21578 0.236539
- 2562 0.0280847
? 0 0
G+C 47499 0.520685
--------- Nucl TABLE ----------
Gene1 (processed):
josephwb@Potiphar-Breen:issue$ pxrms -s gene1_trimAl.txt -f List.txt | pxlssq
File type: fasta
Number of sequences: 12
Is aligned: true
Sequence length: 1086
--------- Nucl TABLE ----------
Nucl Total Proportion
A 2882 0.221148
C 3698 0.283763
G 3261 0.25023
T 3182 0.244168
- 9 0.000690608
? 0 0
G+C 6959 0.533993
--------- Nucl TABLE ----------
Gene2 (original):
josephwb@Potiphar-Breen:issue$ pxlssq -s gene2_trimAl.txt
File type: fasta
Number of sequences: 83
Is aligned: true
Sequence length: 945
--------- Nucl TABLE ----------
Nucl Total Proportion
A 15583 0.198674
C 19734 0.251597
G 18214 0.232218
T 15669 0.199771
- 9235 0.117741
? 0 0
G+C 37948 0.483815
--------- Nucl TABLE ----------
Gene2 (processed):
josephwb@Potiphar-Breen:issue$ pxrms -s gene2_trimAl.txt -f List.txt | pxlssq
File type: fasta
Number of sequences: 12
Is aligned: true
Sequence length: 945
--------- Nucl TABLE ----------
Nucl Total Proportion
A 2294 0.202293
C 2850 0.251323
G 2545 0.224427
T 2352 0.207407
- 1299 0.11455
? 0 0
G+C 5395 0.47575
--------- Nucl TABLE ----------
from phyx.
I note that gene2_trimAl.txt seems to have Windows line breaks:
josephwb@Potiphar-Breen:issue$ file gene1_trimAl.txt
gene1_trimAl.txt: ASCII text, with very long lines
josephwb@Potiphar-Breen:issue$ file gene2_trimAl.txt
gene2_trimAl.txt: ASCII text, with very long lines, with CRLF line terminators
Gene2 runs fine for me (I think? I still do not know why your results are considered wrong). However, you can convert the line breaks and retry things. You will need the program dos2unix
which on linuc you get with:
sudo apt-get install dos2unix
then do:
dos2unix gene2_trimAl.txt
You can check that things worked:
josephwb@Potiphar-Breen:issue$ file gene2_trimAl.txt
gene2_trimAl.txt: ASCII text, with very long lines
Finally, retry your pxrms
command.
from phyx.
Dear Joseph. Thanks a lot for your meticulous reply! I will try again.
from phyx.
Closing this since there has been no activity dum du jaroj .
from phyx.
Related Issues (20)
- phylip gets corrupted when piped HOT 3
- pxtlate detects AA when DNA is present HOT 2
- duplicated species in clade groupings with pxbp
- How to install phyx on linux without root access ? HOT 5
- pxrlt and pxrls only change the final name on the list HOT 9
- pxrlt and pxrmt doesn't output tip name with "-" correctly HOT 5
- pxaa2cdn -r not working as expected HOT 8
- strange behavior with pxbp on computing cluster HOT 7
- pxcomp outputs an enormous file, or rejects unaligned sequences. HOT 6
- phyx installed HOT 11
- pxaa2cdn -r (again) HOT 10
- infinite loop in superdouble implementation with gcc 13 HOT 3
- Unable to ./configure HOT 4
- Debian 12 install, pxseqgen HOT 3
- autoconf error in installation HOT 5
- pxssplit: command not found HOT 6
- pxrms .txt file using wildcard HOT 6
- Phyx errot: double free or corruption (!prev) HOT 4
- GHA failing for mac (armadillo) HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from phyx.