How to use spades operate a fasta file, especially simulated reads in fasta format from badread?
SoSongzhi opened this issue · 8 comments
Description of bug
(bio) zhisong@zhisongs-MacBook-Air spades % bin/spades.py -s ../accuvirnew/bad_nano/hiv_50_2k.fastq -o testfasta --isolate
== Error == /Users/zhisong/Desktop/accuvirnew/bad_nano/hiv_50_2k.fastq is not a valid FASTQ file
Traceback (most recent call last):
File "/Users/zhisong/Desktop/spades/share/spades/spades_pipeline/support.py", line 146, in check_file_not_empty
if next(reads_iterator, None) is None:
File "/Users/zhisong/Desktop/spades/share/spades/spades_pipeline/common/SeqIO.py", line 123, in parse_fastq
raise Exception("/Users/zhisong/Desktop/accuvirnew/bad_nano/hiv_50_2k.fastq is not a valid FASTQ file")
Exception: /Users/zhisong/Desktop/accuvirnew/bad_nano/hiv_50_2k.fastq is not a valid FASTQ file
In case you have troubles running SPAdes, you can report an issue on our GitHub repository github.com/ablab/spades
Please provide us with params.txt and spades.log files from the output directory.
spades.log
no log
params.txt
bin/spades.py -s ../accuvirnew/bad_nano/hiv_50_2k.fastq -o testfasta --isolate
SPAdes version
Spades 3.16.0-dev
Operating System
macOS Ventura version 13.6
Python Version
python 3.8.4
Method of SPAdes installation
manual
No errors reported in spades.log
- Yes
Will you please attach the input file? Per error message – it is not a valid FASTQ file
thank you for your reply. the input is like this:
Please attach the file itself
ok, thank you so much.
hiv_50_2k.fastq.zip
I cannot reproduce the problem, everything is ok for me. So, likely there is some local problem. In any case, these reads are not valid SPAdes input as Illumina reads cannot be that long.
OK，thanks for your reply. And I want to know what if I want to assembly a short sequence like virus, how should I use these nanopore reads, thanks.
Per SPAdes manual (https://ablab.github.io/spades/input.html):
To run SPAdes you need at least one library of the following types:
- Illumina paired-end/high-quality mate-pairs/unpaired reads
- IonTorrent paired-end/high-quality mate-pairs/unpaired reads
- PacBio CCS reads
Illumina and IonTorrent libraries should not be assembled together. All other types of input data are compatible. SPAdes should not be used if only PacBio CLR, Oxford Nanopore, Sanger reads or additional contigs are available.
thanks so much.