Entering edit mode
12 months ago
gkarere
•
0
I use the this command line to run mapper.pl:
/Users/sysadmin/desktop/mirdeep2-0.1.0/src/mapper.pl study/miRNAseq-config_trim2.txt -d -e -h -i -j -m -p genome/hg19 -q -s miRNA_reads_Trim2.fa -t miRNA_reads_vs_hg19.arf -v
At the end of running all files, I get these messages:
#desc total mapped unmapped %mapped %unmapped
Use of uninitialized value $count2 in subtraction (-) at /Users/sysadmin/desktop/mirdeep2-0.1.0/src/mapper.pl line 731.
total: 26314065 Use of uninitialized value $count2 in print at /Users/sysadmin/desktop/mirdeep2-0.1.0/src/mapper.pl line 731.
26314065 Use of uninitialized value $count2 in division (/) at /Users/sysadmin/desktop/mirdeep2-0.1.0/src/mapper.pl line 732.
Use of uninitialized value $count2 in division (/) at /Users/sysadmin/desktop/mirdeep2-0.1.0/src/mapper.pl line 732.
0.000 1.000
This is an example of message in mapper log:
parsing fastq to fasta format
fastq2fasta.pl study/sRNA_99_S99_R1_001_trim2.fastq > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads.fa
converting rna to dna alphabet
rna2dna.pl dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads.fa > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_dna.fa
discarding sequences with non-canonical letters
fastaparse.pl dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_dna.fa -b > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_letters.fa 2>dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_discarded.fa
discarding short reads
fastaparse.pl dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_letters.fa -a 18 > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_no_short.fa 2>dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_too_short
collapsing reads
collapse_reads_md.pl dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_no_short.fa 188 > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_nr.fa
mapping reads to genome index
bowtie -p 1 -f -n 1 -e 80 -l 18 -a -m 5 --best --strata genome/hg19 --al dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/sRNA_99_S99_R1_001_trim2.fastq_mapped --un dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/sRNA_99_S99_R1_001_trim2.fastq_not_mapped dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/reads_nr.fa dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/mappings.bwt 2>bowtie.log
convert_bowtie_output.pl dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/mappings.bwt > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/mappings.arf
trimming unmapped nts in the 3' ends
parse_mappings.pl dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/mappings.arf -j > dir_mapper_188_sRNA_99_S99_R1_001_trim2.fastq_1479928259_22_12_2023_t_11_21_48/mappings_trim.arf
The error is after running python3.12 cutadapt
You have been warned not to add answers when you're not actually answering and to use the forum properly. Invest some effort into learning.
I finally got the mirdeep2-0.1.3 installed successfully in Mac computer (Sonoma 14.2.1).
I run this mapping script:
First I get this message:
Then, this run:
Results: the "arf" folder has zero bytes and "fa" folder has 250MB of data. Statistics shows no mapping was done.
Does anyone know why arf folder is zero after the run and what I can do?