minimap to make a comparison between reference to de novo assembly
0
0
Entering edit mode
2.9 years ago

Hi

I am trying to compare between my reference fasta file versus some de novo assembly (that includes just a sequences from family of proteins in fasta format after run blast).

I am using the following commands:

minimap2 -cx asm5 --cs /shared2/salmosim/lissa/fasta/TriTrypDB-41_TcruziSylvioX10-1_Genome.fasta CG_1.fasta

when I try to get a stats, I have that results

Number of mapped sequences: 990
Number of primary alignments: 7773
Number of secondary alignments: 7527
Number of primary alignments with >65535 CIGAR operations: 0
Number of bases in mapped sequences: 357654
Number of mapped bases: 489493
Number of insertions in [0,50): 12706
Number of insertions in [50,100): 10
Number of insertions in [100,300): 0
Number of insertions in [300,400): 0
Number of insertions in [400,1000): 0
Number of insertions in [1000,inf): 0
Number of deletions in [0,50): 661
Number of deletions in [50,100): 7
Number of deletions in [100,300): 0
Number of deletions in [300,400): 0
Number of deletions in [400,1000): 0
Number of deletions in [1000,inf): 0

Bur if try to call the SV everyone is 0.

minimap2 -cx asm5 --cs /shared2/salmosim/lissa/fasta/TriTrypDB-41_TcruziSylvioX10-1_Genome.fasta CG_1.fasta CG_1.fasta | sort -k6,6 -k8,8n | paftools.js call -L20000 - > var.txt

I am a beggining in this software, so actually i don't have idea what is happening.

All you help could be a lot of useful to me.

Thank you so much.

MINIMAP alignment • 867 views
ADD COMMENT
0
Entering edit mode

Hi Lissa Cruz Saavedra, I added some markdown to your post to highlight code and data examples using the code option (10101) from the format bar. Also, please avoid capslock in the title as this is generally considered impolite because it implies shouting.

ADD REPLY
0
Entering edit mode

Thank you so so much!. I didn't know, your so kind!

ADD REPLY

Login before adding your answer.

Traffic: 1868 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