poretools icon indicating copy to clipboard operation
poretools copied to clipboard

poretools not finding fastq after subsequent albacore basecalling

Open biorover opened this issue 7 years ago • 3 comments

I ran a rapid library kit (RAD002) with basecalling, but stopped it before it could catch up. I then basecalled the skip folder using albacore. Poretools can find the sequences of the fast5 files basecalled by MinKNOW, but can't find the sequences basecalled by albacore. Upon inspection with h5ls, I found that Albacore had written the basecalls to Analyses/Basecall_1D_001/BaseCalled_template/ (Analyses/Basecall_1D_000/ only contains "Summary" group). However, from my cursory glance at Fast5File.py, it should be able to find basecalls no matter what the numbers on the end of "Basecall_1D_***" are. Any idea what's going on? Sample fast5 file attached. albacorecalledread.fast5.zip

biorover avatar Mar 23 '17 22:03 biorover

Use --group 1

nickloman avatar Mar 23 '17 22:03 nickloman

Works perfectly, thanks!

biorover avatar Mar 23 '17 23:03 biorover

Thanks, this was helpful

AhmedAbdelfattah avatar Apr 14 '17 09:04 AhmedAbdelfattah