I tried to use poretoos to convert fast5 to fasta file but the below command did not create any file
> poretools fasta --type best 15
> poretools fasta --type best 15/*
> poretools fasta --type best 15/*.fast5
What did I miss you?
Thank you in advance,
•
written
by
Ric
1 answer
It may be that your files don't have the type of reads that best extracts. Have you tried using some of the other types, for example 2D, or no type at all? If that outputs something, it's likely that your file doesn't have any reads that match the type best, which doesn't seem to make much sense -- the poretools documentation seems to imply that something should be output when using best:
A type of “best” will extract the 2D read, if it exists. If not, it will extract either the template or complement read, whichever is available and has a better average Phred score.
Also keep in mind that the fasta file will be output to stdout; if you'd like them in a file, you'll have to redirect the output.
Log in to answer this question.
-
Convert Nanopore Fast5 files to Fasta format
written by Sowmya Pulapet •
-
How to convert base called fast5 to fastq
written by YANG •
-
Getting 0 bytes or 1 line from galaxy
written by dreamfeathers08 •
-
Getting an empty .fastq file as output when ran poretools command on multiple .fast5 files
written by SV •
-
Only 1 linkage group after filtering MSTmap result
written by Ric
-
samtools sort: truncated file. Aborting
written by m.t.lorenc •
-
faSize caused expandFaFastBuf
written by Ric
-
vcf to venn diagram
written by Ric
-
Failed to open var.raw.bcf: unknown file type
written by Ric
-
new using pore tools. extract fast
written by lagartija
How did you obtain these fast5 files? Are those basecalled?
By the way, since the basecallers can output fastq directly
poretoolsis, well, obsolete now.We ask sequence company to sequence our genome. Is there any other software available?
Then you could ask that company what they send you: raw or base called fast5. Or you could share one fast5 file with me/us and I can have a look. In any case, you can do the basecalling yourself, using
albacore, which you can get from the nanopore community forum.