Soru:
Albacore çalıştırılıyor, ancak 0 okuma çıktı veriyor
d_kennetz
2019-01-23 02:26:15 UTC
view on stackexchange narkive permalink

MinION tarafından SQK-LSK109 kitini ve FLO-MIN106 akış hücresini kullanarak komut satırı aracılığıyla üretilen arama verilerini temel almaya çalışıyorum. Albacore sürümüm en son sürümdür (v2.3.4). Aşağıdaki komut satırını kullanarak çalıştırmayı denedim (dizileme çalışmasının ve fast5s'nin adı farklı ve makine çıktısından varsayılan):

  python3.6 /home/dkennetz/.local/bin /read_fast5_basecaller.py -i ./seqrun/fast5/ -t 1 -s ./good_reads/ -f FLO-MIN106 -k SQK-LSK109 -r -o fastq  

-t 1 sadece bu test içindir. Çalışıyor ve neredeyse anında 0 hata ile bitiyor, ancak hiçbir okuma çıktı vermiyor. Daha önce aynı komut satırını SQK-LSK108 kitini kullanarak çalıştırdım ve hiçbir sorun yaşamadım.

Daha sonra h5_validate aracını kullanarak fast5'leri doğrulamaya başladım:

  python3.6 /home/dkennetz/.local/bin/h5_validate /home/dkennetz/.local/lib/python3.6/site-packages/h5_validator/schemas/single_read_fast5.yaml seqrun / fast5 / 0 / seqrun.fast5 -v 

Bunun çıktısı:

  Matcher <'UniqueGlobalKey ', type = group match_type = exact count = {' minimum_count ': 1,' maximum_count ': 1} <HDF5 dosyası "seqrun.fast5" (mode r) >Error at /: Matcher <'PreviousReadInfo', type = group match_type = exact count = {'minimum_count': 1, 'maximum_count' ile eşleştikten sonra > tatmin olmadı 1} >, <HDF5 dosyası "seqrun.fast5" (mode r) >  
ile eşleştikten sonra tatmin olmadı

Topluluk sayfasında bunun, başlıkta bazı ilgili meta eksikliğinden kaynaklanabileceğini okudum veri bilgisi? Bunun nasıl düzeltileceğiyle ilgili herhangi bir yardım çok takdir edilecektir!

Bu, ONT için topluluk sayfasında çapraz olarak yayınlanmıştır:

https: //community.nanoporetech. com / posts / albacore-basecalling-runni

Bir cevap:
Wouter De Coster
2019-01-23 13:07:54 UTC
view on stackexchange narkive permalink

MinKNOW’un son sürümleri tarafından oluşturulmuş çok hızlı5 dosyalarınız olduğuna inanıyorum. Bunları daha eski, tek okumalı tek dosya biçimine dönüştürmek için bir dönüştürme komut dosyası mevcuttur. AFAIK Albacore kullanımdan kaldırılıyor ve Albacore Guppy ile değiştirildiği için çok kafalı fast5 dosyalarını destekleyen hiçbir sürüm yapılmayacak.

Çoklu okuma fast5s'den tek okuma fast5s'e dönüştürmek için komut satırı aşağıdaki gibidir:

  $ multi_to_single_fast5 -i / yol / yol / hızlı5 / -s / yol / / kaydet / çıktı /  
Merhaba Wouter, cevabınız için teşekkürler! Bu beklendiği gibi çalıştı. Herhangi bir belirsizliği gidermek için, gelecekteki izleyiciler için komut satırını ekleyerek düzenledim. Kabul edip etmemeyi seçebilirsiniz, sadece netlik kattığını düşündüm! Bağlantı için minnettarım.


Bu Soru-Cevap, otomatik olarak İngilizce dilinden çevrilmiştir.Orijinal içerik, dağıtıldığı cc by-sa 4.0 lisansı için teşekkür ettiğimiz stackexchange'ta mevcuttur.
Loading...