Comments (14)
Hi,
Yes, ema preproc
should produce 500 (-n
) files of the form ema-bin-000
, ema-bin-001
etc. (the ???
is a bash idiom for matching files of this form). I wonder if there is some quirk with that FASTQ file -- could you post the first 10 or so lines of it?
from ema.
from ema.
That FASTQ looks good -- how do the *.ema-ncnt
files look? Are they produced by ema count
?
from ema.
No clue how the .ema-ncnt files really look. They are not text. I have only two files: counts_file.ema-ncnt and counts_file.ema-fcnt. The following is the first 10 lines of counts_file.ema-ncnt:
� *
9���
-���
��
����
�
����
���
���
��c�
The log for the count is:
:: Loading 10X took 2.4 s
:: Dumped block 1
:: Dumped block 2
:: Dumped block 3
:: Dumped block 4
:: Dumped block 5
:: Dumped block 6
:: Dumped block 7
:: Dumped block 8
:: Counting took 2472.2 s
:: Reads with OK barcode: �Ò637,148,268 out of Ð683,891,843
:: Ignored th OK barcode: �Ò637,1480 reads
:: Dumped block 9
:: Printing took 8.6 s
:: Processed oo683,891,843 reads ( 470,201 MB uncompressed) in 2,483 s
--Thanks for your help!
from ema.
Can you please try doing this instead:
cat ~/EMA/some_clam_S3_L004_interleaved.fastq | ema preproc -w ~/EMA/4M-with-alts-february-2016.txt -n 500 -t 40 -o test *.ema-ncnt 2>&1 | tee preproc.log
Istest
directory empty or not? I have a feeling that passing -o ~/EMA
confuses the script because of ~
.
from ema.
from ema.
So the log contains all these non-ASCII characters? Weird... any chance you can share the first 10k records from the interleaved file so that I can take a look?
from ema.
from ema.
OK, that is certainly unexpected.
Which OS are you using? You can email me the data link (inumanag at mit dot edu
--- first 10k reads should be fine) and I can take a look. Currently, we don't have any small sample for these purposes.
from ema.
from ema.
Just following up on this one: did you manage to resolve this issue with another compiler or pre-compiled binaries?
from ema.
from ema.
Yeah, I believe that there is something wrong with the compiler on your cluster because the generated binary is completely corrupted--- probably because some compiler library versions got mixed up (this happens sometimes on cluster environments).
As for Long Ranger vs. EMA: yeah, they are. We found out that EMA's results are more robust, esp. in repetitive regions, so we always map with EMA.
from ema.
I have now given it one more shot and installed using conda. I can see the help text when I type "ema help" but the command yields a completely empty log file and no count files at all and no errors shown.
cat ~/EMA/some_clam_S3_L004_interleaved.fastq | ema count -w ~/EMA/4M-with-alts-february-2016.txt -o counts_file 2>&1 | tee count.log
Any thoughts?
from ema.
Related Issues (20)
- Option for preprocess and map to write and read gzipped files? HOT 2
- Interpreting read counts HOT 1
- reference lenght constraint HOT 2
- Compilation error HOT 2
- ema align core dumped
- How to use "-b" flag? HOT 2
- Platform Flag HOT 1
- Input and workflow unclear HOT 6
- 8 bp whitelist HOT 1
- ema align output SAM parsing error HOT 4
- Why is read density optimization off by default?
- What's the meaning of density probs?
- ema align errors
- Can it used for read of length 500bp ?
- EMA split-read handling
- `ema align` silently fails with `[noeol]` FASTA files HOT 1
- Not all input reads are included in `ema align` output
- Barcodes not encoded/decoded correctly for haplotag reads.
- Incorrect template length output
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 ema.