artic-network / artic-ncov2019

ARTIC nanopore protocol for nCoV2019 novel coronavirus
Creative Commons Attribution 4.0 International
168 stars 166 forks source link

how to get my consensus sequence using artic-network MINion work flow using metagenomic data #47

Open mhelbom14 opened 3 years ago

mhelbom14 commented 3 years ago

Hello,

I used a metagenomic approach to sequence SARS-CoV2 genome and the data were very good with 100% coverage and 422x depth converge. now I want to assemble my reads and get the consensus sequence and the variant file. I want to use the workflow described here but I cant optimize my code to work with my data. I did the basecalling and the demultiplexing for the data and I have the pass fastq files. I want to start from the artic minion command. I tried to run the command with my data but got this error. " (artic-ncov2019) lab123@DESKTOP-2AI9TES:/mnt/d/Mohammad/Nanopore/Analysis/covid-19$ artic minion --normalise 200 --threads 4 --read-file /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01f475dc98*.fastq --fast5-directory /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fast5_pass/barcode01/ --sequencing-summary /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/sequencing_summary_FAL01681_f475dc98.txt usage: artic [-h] [-v] {extract,basecaller,demultiplex,minion,gather,guppyplex,filter,rampart,export,run} ... artic: error: unrecognized arguments: /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_100.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_101.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_102.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_103.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_104.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_105.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_106.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_107.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_108.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_109.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_11.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_110.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_12.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_13.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_14.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_15.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_16.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_17.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_18.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_19.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_2.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_20.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_21.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_22.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_23.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_24.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_25.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_26.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_27.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_28.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_29.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_3.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_30.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_31.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_32.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_33.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_34.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_35.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_36.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_37.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_38.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_39.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_4.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_40.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_41.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_42.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_43.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_44.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_45.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_46.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_47.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_48.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_49.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_5.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_50.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_51.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_52.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_53.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_54.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_55.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_56.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_57.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_58.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_59.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_6.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_60.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_61.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_62.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_63.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_64.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_65.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_66.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_67.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_68.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_69.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_7.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_70.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_71.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_72.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_73.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_74.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_75.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_76.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_77.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_78.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_79.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_8.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_80.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_81.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_82.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_83.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_84.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_85.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_86.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_87.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_88.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_89.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_9.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_90.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_91.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_92.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_93.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_94.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_95.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_96.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_97.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_98.fastq /mnt/c/data/covid-19/covid-19/1/20201102_1452_MN31812_FAL01681_d9c4ed35/fastq_pass/barcode01/FAL01681_pass_barcode01_f475dc98_99.fastq " how could I fix it? and for the "scheme-directory" I don't have scheme directory because I didn't use artic network primers. I used SISPA A/B method

will-rowe commented 3 years ago

Hi @mhelbom14. Thanks for getting in touch. Unfortunately this pipeline is designed to work with amplicon sequencing data, specifically the Network ARTIC schemes, and doesn't perform an assembly.

Maybe something like VGEA would be more suitable?

mhelbom14 commented 3 years ago

Thank you so much @will-rowe I mean the generation of the consensus sequence is it possible to edit the script so it would work with metagenomic data? I will try VGEA