Comments (15)
Hi Robert,
Sorry for my super late reply, finally it did work, I reinstalled Racon again and it run without problems. Thanks!
from racon.
Hello,
could you please provide the full command which caused this error?
Best regards,
Robert
from racon.
from racon.
I am experiencing the same problem. I first run Graphmap:
/graphmap/bin/Linux-x64/graphmap align -r scaffolds.fasta -d long_reads.fastq -o long_reads_scaffolds.sam
And then I run Racon:
racon/bin/racon --sam long_reads.fastq long_reads_scaffolds.sam scaffolds.fasta nanopore_racon.fasta
And it gives me this error:
ERROR: Reads are not specified in a format which contains quality information. Exiting
However, I have checked my fastq file with the long reads and indeed, is in fastq format, any suggestion?
Thanks!
from racon.
Hello Elena,
I executed your commands on a small dataset containing lambda reads in FASTQ format and everything worked fine. Are you sure that your long_reads.fastq file is in a valid FASTQ format? If you are allowed to, please send me the output of the following commands: "head -n 16 long_reads.fastq" and "tail -n 16 long_reads.fastq".
Sorry for the delayed answer and best regards,
Robert
from racon.
Hello Robert,
Here I send you the the head and tail of my long reads. It looks fine I think, there are as many quality characters as nucleotides in the sequences in all the cases. So maybe the problem would be that there are some sequences in the middle without quality values? It his would be the case, I dont know why did it happen...
output_head_tail_long_reads_fastq.txt
from racon.
Hello Elena,
the sample file you provided looks fine. How did you obtain the whole long_reads.fastq file?
Best regards,
Robert
from racon.
Hello Robert,
My long reads fastq file is the result of 5 MinION flowcells. Every run output 1 fastq file, and I concatenated the 5 fastq files to obtain the long_reads.fastq file. I also think the file looks OK, so I guess there are some "intruder" sequences in my file without quality values or with uncorred number of quality values...its the obly explanation. Do you know how can I chekc this? I read somewhere about fastQValidator..
Best regards,
Elena
from racon.
Hello Elena,
I wrote you a simple python script for FASTQ validity (https://pastebin.com/6anD3d2q). Copy the code in a file and run it with 'python <file> long_reads.fastq'.
Best regards,
Robert
from racon.
Thanks! I will let you know that if that was the issue and if it works :) Thanks again!
Elena
from racon.
Hello Robert,
I run it with long_reads.fastq, it runs, it finishes and it does not say anything. I tried at well with the head and tail I gave you, but with one quality line trimmed on purpose, and it points it out. So, I assume the fastq file is ok...
Does the sam file need to have quality values as well or is enough with the read file having them? Anyway, when I used Graphmap I mapped the scaffolds against that same long_reads.fastq file, so I assume it should have them...
from racon.
Hello Elena,
I am not sure what the problem is. Have you the latest graphmap version? Did graphmap finish without errors? Please paste here the output of the command 'head -n 15 long_reads_scaffolds.sam' as well as 'tail -n 15 long_reads_scaffolds.sam'.
Best regards,
Robert
from racon.
Hello Robert,
No, graphmap did not give any error, I think I used the last version but I have just installed again and I am repeating the alignment. It should not be something related to the memory, right? The long_reads.fastq is 14.64 Gb, long_reads_scaffolds.sam is 19.04Gb and the scaffolds.fasta is 1.646 Gb.
I have just checked the head and tail of the sam file, they look different, is this normal?
/* rvaser removed the snip */
Thanks!!!
Elena
from racon.
Hello Elena,
the sam file looks okay to me. If you fail to run racon with the repeated alignments, please try minimap instead of graphmap.
Best regards,
Robert
from racon.
Hello Elena,
were you able to run racon one way or another?
Best regards,
Robert
from racon.
Related Issues (20)
- The number contig decrease after polishing HOT 1
- Compiling RACON for CUDA 11.4 & <$ make package> FAILS on Ubuntu 18.04 - SOLVED
- Empty Output file when using Racon HOT 7
- error: empty overlap set! and multiple consusnses sequnces HOT 6
- Issue with some files resulting in "overlap is not transmuted" HOT 3
- unequal length error HOT 1
- Improper HOT 1
- terminate called after throwing an instance of 'std::invalid_argument' what(): [bioparser::SamParser] error: invalid file format HOT 1
- Segmentation Fault (core dumped) HOT 1
- Output record of fixes HOT 1
- VERY slow performance with fastq.gz files HOT 3
- Racon error: empty overlap set!
- polishing contig with Nanopore Q20 read? HOT 1
- Use cDNA to polish exons
- Racon cannot work with duplucates in paf file HOT 1
- How to generate OVERLAP file? HOT 1
- installation on CentOS 7: 'Could NOT find ZLIB: Found unsuitable version "1.2.7", but required is at least "1.2.8"' HOT 1
- Empty Output (non-memory-related issue) HOT 3
- error: overlap is not transmuted!
- racon: invalid option -- 'L' HOT 3
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 racon.