-
Notifications
You must be signed in to change notification settings - Fork 50
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
[bioparser::FastqParser] error: invalid file format! #99
Comments
Hello Huandna, Best regards, |
|
Aren't you missing |
I am sorry it is so long ,I am sure I have ‘>’ for each sequence |
If the above snippet is the result of the head command, then you are missing |
oh,Maybe the problem is that the results displayed on the web page are different from those in Linux. I will find a standard fasta file to test again Thank you very much |
The raw text has them, my bad :) Can you run |
|
That looks fine as well, hmm. Can you run the same for |
this is
and the
I think the reason maybe |
The FASTA file can be wrapped, i.e. the new lines are not the problem. Well, that depends on which version of racon are you using. Which is it? :) |
oh ,my version is v1.3.1 |
Well, I have no idea then. Would you mind sharing the data with me via email so I can investigate further? |
You can also try and check if |
I find the issue Many thanks for your offer |
I ran into this same error (and solved it, but I just wanted to leave this here in case it's useful for anybody else in the future) In my case it was because my fastq file contained reads of length 0. Performing a trimming step with |
I'm having this same issue. I haven't been able to figure it out. I'm using racon version 1.3.3, installed via git clone. I ran it as part of an iterative script on a cluster. Out of 6 datasets, 5 worked fine. One didnt. Command and error given are: Command: ~/Tools/racon/build/bin/racon ../sks125VsHuman_unmapped.fastq sks125_readsVsCleanAssem.sam blob/sks125_clean.contigs.fasta -t 32 > ../sks125_racon1.fasta I've even gone back to see if it was something wrong with one my processing steps and its still the same. |
Hello, Sorry for the late response! |
Hi Robert, I have the same problem. Could you please help me out? This is my racon command:
The error file:
So far I have troubleshooted by:
I am using an HPC. Racon was installed using this https://anaconda.org/bioconda/racon. Any help would be greatly appreciated. Cheers, Camilo |
Hello Camilo, Best regards, |
Hi Robert, Yes, this is the command from minimap2:
I realize I use the raw read data here instead of the trimmed reads. Do you think that is the problem? Best, Camilo |
Without option |
Hi Robert, I ran the following command:
I did not get the same error above, but I got a different one after it run for a few minutes:
Do you have any suggestion? Cheers, Camilo |
Try using the original read file |
Thanks Robert! It is running. So do you recommend to stick to Cheers, Camilo |
However you want :) If you want |
Apologies for being so late. Yes it turns out there was an issue with one the reads in the fasta file. Just had to format it and it worked fine! |
Hi Rvaser, i have been very curious comparing racon with hypo polish. Hypo work well but i keep getting this problem: [racon::Polisher::initialize] loaded target sequences 8.372787 s [racon::Polisher::initialize] loaded sequences 2478.189205 s terminate called after throwing an instance of 'std::invalid_argument' what(): [bioparser::SamParser] error: invalid file format Aborted (core dumped), troubleshooting with raised issues doesn't help me: please can you help me with some ideas this is by command: i am using illumina paired end reads minimap2 -t 64 consensus.fasta ${reads} > illumina.samracon -t 48 CsM2.fastq illumina.sam draft.fa > polished_Male3.racon.faAlternatively i tried overalapping the reads with generate and use a paf. file still it didn't work, minimap2 -x ava-pb -t 8 reads.fa reads.fq > raw_hifiam.paf Please i need your help. thanks, Raissa |
i used your script to generate a single file from the paired end reads as input for my sequence file but sill not working |
Hi rvaser, |
If you want a |
I had tried everything else and this was the fix for me as well! For context, I'm using Minimap2 version 2.24 at time of posting. I am wondering if something with printing to STDOUT (default of minimap2) is causing formatting issues for Racon. Either way, directing an output using -o was the solution, so thank you for saving me the headache! |
I tried use fasta file from nanopore to polish my genome of ecoli with racon v1.3.1 here is my code
and I get the error
[racon::Polisher::initialize] loaded target sequences [bioparser::FastaParser] error: invalid file format!
Then I used the same fasta file to simulate a fastq file,and get the same error:
[racon::Polisher::initialize] loaded target sequences [bioparser::FastqParser] error: invalid file format!
I used fastq as racon material without any problems before. this is my first time to use fasta file ,Is there anything different Between them? or is it because my genome is too small? Ecoli only has one sequence (4MB) and my fasta file is 1.5GB
thanks for your time
huandna
The text was updated successfully, but these errors were encountered: