spundhir / rna-seq Goto Github PK

Hello,

I am trying to map small RNA sequencing dataset to noncode v6 human fasta file. I have paired end reads for each samples. When I am running following command it is giving me this error. Please help!!

segemehl.x -i NONCODEv6_human.idx -d NONCODEv6_human.fa -q /data/user/mbansal/AG_DKC1_rawdata/small_rna_seq/Output_Teton-1_R1_001.trimmed.fastq -p /data/user/mbansal/AG_DKC1_rawdata/small_rna_seq/Output_Teton-1_R2_001.trimmed.fastq > Teton-1.sam

[SEGEMEHL] Wed Jan 25 10:40:09 2023: reading queries in '/data/user/mbansal/AG_DKC1_rawdata/small_rna_seq/Output_Teton-1_R1_001.trimmed.fastq'.
[SEGEMEHL] Wed Jan 25 10:40:16 2023: 5060435 query sequences found.
[SEGEMEHL] Wed Jan 25 10:40:16 2023: reading mates in '/data/user/mbansal/AG_DKC1_rawdata/small_rna_seq/Output_Teton-1_R2_001.trimmed.fastq'.
[SEGEMEHL] Wed Jan 25 10:40:22 2023: 1: Reading mates failed: mate and query files differ in size!
[SEGEMEHL] Wed Jan 25 10:40:22 2023: set->nooffiles 1 = 1 nooffiles
[SEGEMEHL] Wed Jan 25 10:40:22 2023: set->noofseqs 5060435 = 5038367 noofseqs
[SEGEMEHL] Wed Jan 25 10:40:22 2023: set->chunkindex->size 506 = 503 chunkindex->size(segemehl)