converting fast5 to fasta failed

Question

[Ric]

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,

[WouterDeCoster]

How did you obtain these fast5 files? Are those basecalled?

By the way, since the basecallers can output fastq directly poretools is, well, obsolete now.

[Ric]

We ask sequence company to sequence our genome. Is there any other software available?

[WouterDeCoster]

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.

Answer 1

[Heather Ward]

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.

[Ric]

I tried all and 2d none of them have created an output.