Minimap2 error: [ERROR] dual gap penalties violating E1>E2 and O1+E1<O2+E2
0
0
Entering edit mode
4 weeks ago

Hi,

I'm trying to run minimap2 to assemble Oxford Nanopore reads to a reference but I keep getting the following error: [ERROR] dual gap penalties violating E1>E2 and O1+E1<O2+E2.

Here is the command I'm running:

minimap2 -s 10 -t 4 -k 19 -B 4 -O 6 -R '"@RG\tID:"$samplename"\tSM:"$samplename"\tLB:"$samplename""' '/home/alexbyrne/mnt/VI6Storage/Raw_NGS_Data/MinION_Data/1353_Test_Run/FAO54735_SPAdes_Influenza_Database_26-10-21/top_matches.fa' '/home/alexbyrne/mnt/VI6Storage/Raw_NGS_Data/MinION_Data/1353_Test_Run/FAO54735_pass_c4ec9f31_merged.fastq' | samtools view -@ 4 -u - | samtools sort -@ 4 -o '/home/alexbyrne/mnt/VI6Storage/Raw_NGS_Data/MinION_Data/1353_Test_Run/FAO54735_pass_c4ec9f31_map_sorted.bam'

I have tried altering the values for the different flags (-s, -k, -B and -O) but I keep getting the same error.

Do you know what could be causing this issue?

assembly Oxford minimap2 Nanopore • 258 views
ADD COMMENT
1
Entering edit mode

Did you try with the preset -x map-ont?

ADD REPLY
0
Entering edit mode

Thanks!

I tried adding -ax map-ont to the above options and got the same error, but if I remove the other options then it work ok.

Is -ax map-ont sufficient to replace the other options?

ADD REPLY
0
Entering edit mode

Just an added query, I checked the output bam file from the above using picard ValidateSamFile and its missing the read groups. Is there a reason the -R flag is being ignored?

ADD REPLY
0
Entering edit mode

I don't know if there are other issues, but you enclosed the RG tags inside single quotes, so the variable $samplename within the single quotes will be taken as literal, and won't be interpolated to its value.

On a side note, it doesn't seem the RG fields you are adding are correct (though they may be). In general, several samples are prepared simultaneously with one library kit, so having $samplename as LB filed is odd.

ADD REPLY
0
Entering edit mode

Is -ax map-ont sufficient to replace the other options?

Should be: it is a preset set by the developer of minimap2, after extensive testing.

ADD REPLY

Login before adding your answer.

Traffic: 1724 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6