fast5 demultiplex by barcode

fast5 demultiplex by barcode

pip3 install ont-fast5-api

cat sequencing_summary.txt | cut -f 21 | grep barcode[0-9] | sort | uniq > barcodes.log

while read -r barcode
  do
    head -n 1 sequencing_summary.txt > $barcode.txt
    cat sequencing_summary.txt | grep $barcode >> $barcode.txt
  
    for f in *.txt
      do 
        fast5_subset -i folder/with/your/fast5s -s output/folder/for/demultiplexed_fast5/$barcode -l $f -t 14
      done
  done < barcodes.log

Hi, I don't know. Guppy 6 should be demultiplexing your barcodes with the mentioned command. If it is not happening you can try the code in the gist above. Have a look at the sequencing_summary.txt file created by guppy and check if there are any 'barcode0x' in field 21 or if they are all 'unclassified' cat sequencing_summary.txt | cut -f 21 If there are barcodes, you should be able to use the script in this gist, but also, guppy should have done it.... If that leads to nothing you will have to ask this in the nanopore community, as both guppy and this script work for my setup, sorry.

…

On Thu, 17 Feb 2022 at 22:38, Fatih Tiras ***@***.***> wrote: ***@***.**** commented on this gist. ------------------------------ Hi Johannes, Now it looks like working. many thanks. but after basecalled, generated fast5 files are still all together in one file (in the workspace) and named like FAP...fast5_skip.....fast5. is there any way to separate that fast5 files according to each barcodes because ı wanna use deepsignal and tombo packages? many thanks. :) [image: Slayt1] <https://user-images.githubusercontent.com/75265865/154503922-16cb4b6e-ed48-4ef4-af4e-76fec4ad244e.JPG> Sincerely. — Reply to this email directly, view it on GitHub <https://gist.github.com/4af1e37ac96af3d9ed2104658993d944#gistcomment-4069063>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABBHB2WSPLGAKVRUNIPBM6TU3UB7ZANCNFSM5OOPDVPA> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: <JWDebler/fast5 demultiplex by ***@***.***>

Hi Johannes,

I am still working to fix my problem but many thanks for your help. I have really appreciated your effort for help. have a great days.:)

Sincerely.

No worries. Let me know once you succeed. I'm interested to find out why my suggestions didn't work. Good luck.

…

On Mon, 21 Feb 2022, 20:14 Fatih Tiras, ***@***.***> wrote: ***@***.**** commented on this gist. ------------------------------ Hi Johannes, I am still working to fix my problem but many thanks for your help. I have really appreciated your effort for help. have a great days.:) Sincerely. — Reply to this email directly, view it on GitHub <https://gist.github.com/4af1e37ac96af3d9ed2104658993d944#gistcomment-4072902>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABBHB2VN7RUE22EZ2WAB3JDU4IUCXANCNFSM5OOPDVPA> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: <JWDebler/fast5 demultiplex by ***@***.***>

Hi Johannes,

I wanna let you know that your code has worked on my last run and ı have been able to separate fast5 for each barcode. I dont know why but ıt has not been working my previous run (still generate - folder and put all fast5 in there like here nanoporetech/ont_fast5_api#68) so ı decided to basecall all my previous raw data again and see how it come up. I wanna personally say many thanks for your effort. have a great days.:)

Best wishes.

Hi Fatih, that's great to hear! Glad this code helped. Not sure why it didn't work at first, but glad you made it work! All the best for your experiment!

…

On Mon, 28 Feb 2022 at 20:28, Fatih Tiras ***@***.***> wrote: ***@***.**** commented on this gist. ------------------------------ Hi Johannes, I wanna let you know that your code has worked on my last run and ı have been able to separate fast5 for each barcode. I dont know why but ıt has not been working my previous run (still generate - folder and put all fast5 in there like here nanoporetech/ont_fast5_api#68 <nanoporetech/ont_fast5_api#68>) so ı decided to basecall all my previous raw data again and see how it come up. I wanna personally say many thanks for your effort. have a great days.:) Best wishes. — Reply to this email directly, view it on GitHub <https://gist.github.com/4af1e37ac96af3d9ed2104658993d944#gistcomment-4080360>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABBHB2S2GFQPYBIFO6RVJ33U5NS75ANCNFSM5OOPDVPA> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: <JWDebler/fast5 demultiplex by ***@***.***>

hi @BeatrizFaustino

I used demux_fast5 --input fast5_skip/ --save_path ./demutiplexedfast5/ --summary_file sequencing_summary.txt command for demultiplex my fast_skip folder for barcoded. You should use Minknown sequencing summary file for input.

Hi @BeatrizFaustino ,

the sequence summary file that ı used is provided after sample runs (sequenced) from MiKnown generally sequence_summarry_FAPxxxxx.txt. it's not generated a file after base calling like a guppy.
Yes fast5_skip includes all fast5 that is generated after sequencing by MinKnown and they don't separate by MinKnown (or you stopped the basecalling after sequencing finished but basecalling was still continuing). I hope this help you. thanks

Hi
is demux present in guppy or other software?

hi @BeatrizFaustino ,

I used this page as a demultiplex ont_fast5_api

Hi @Fatihlrcfs
Thank you very much, using the commands you gave me the analyzes ended now and I was successful

Hi, yes, the latest version of minknow can also demultiplex fast5s, but it has to do it by basecalling. If you have a minknow setup with a GPU I would just recall the raw data with the super accurate model.

…

On Sat, 17 Dec 2022, 06:19 B13F, ***@***.***> wrote: ***@***.**** commented on this gist. ------------------------------ Hi @Fatihlrcfs <https://github.com/Fatihlrcfs> Thank you very much, using the commands you gave me the analyzes ended now and I was successful — Reply to this email directly, view it on GitHub <https://gist.github.com/4af1e37ac96af3d9ed2104658993d944#gistcomment-4405609> or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABBHB2TQ2K2BXD7IXPBWQS3WNTTIHBFKMF2HI4TJMJ2XIZLTSKBKK5TBNR2WLJDHNFZXJJDOMFWWLK3UNBZGKYLEL52HS4DFQKSXMYLMOVS2I5DSOVS2I3TBNVS3W5DIOJSWCZC7OBQXE5DJMNUXAYLOORPWCY3UNF3GS5DZVRZXKYTKMVRXIX3UPFYGLK2HNFZXIQ3PNVWWK3TUUZ2G64DJMNZZDAVEOR4XAZNEM5UXG5FFOZQWY5LFVEYTCMZZG42TSNRTU52HE2LHM5SXFJTDOJSWC5DF> . You are receiving this email because you authored the thread. Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub> .