Skip to content
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

How to use spades operate a fasta file, especially simulated reads in fasta format from badread? #1305

Open
1 task done
SoSongzhi opened this issue May 22, 2024 · 8 comments

Comments

@SoSongzhi
Copy link

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
@asl
Copy link
Member

asl commented May 22, 2024

Will you please attach the input file? Per error message – it is not a valid FASTQ file

@SoSongzhi
Copy link
Author

thank you for your reply. the input is like this:
image

@asl
Copy link
Member

asl commented May 22, 2024

Please attach the file itself

@SoSongzhi
Copy link
Author

ok, thank you so much.
hiv_50_2k.fastq.zip

@asl
Copy link
Member

asl commented May 23, 2024

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.

@SoSongzhi
Copy link
Author

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.

@asl
Copy link
Member

asl commented May 23, 2024

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.

@SoSongzhi
Copy link
Author

thanks so much.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants