Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

The x-axis chromosome ID is missing in the result of pafCoordsDotPlotly.R. #21

Open
lvqiang0120 opened this issue May 20, 2023 · 0 comments

Comments

@lvqiang0120
Copy link

lvqiang0120 commented May 20, 2023

Hi,
I have 11 chromosomes in both my reference and query sequences. However, after using Minimap2 to generate a PAF file and then using pafCoordsDotPlotly.R to plot it, the y-axis displays 11 chromosome IDs while the x-axis only displays 9 chromosome IDs. I am sure that my reference sequence contains all 11 chromosomes. I do not know why two chromosome ID numbers are missing on the x-axis.

There are my command line:
/minimap2 -x asm5 $ref $query -t $thread > $out_dir/"$prefix".minimap2.paf
pafCoordsDotPlotly.R -i "$prefix".minimap2.paf -o "$prefix".drplot -k 11 -q 10000 -m 200000 -r chr1,chr2,chr3,chr4,chr5,chr6,chr7,chr8,chr9,chr10,chr11 -x

Log file prompt:
PARAMETERS:
input (-i): cpxp.minimap2.paf
output (-o): cpxp-x.drplot
minimum query aggregate alignment length (-q): 10000
minimum alignment length (-m): 2e+05
plot size (-p): 15
show horizontal lines (-l): FALSE
number of reference chromosomes to keep (-k): 11
show % identity (-s): FALSE
show % identity for on-target alignments only (-t): FALSE
produce interactive plot (-x): FALSE
reference IDs to keep (-r): chr1,chr2,chr3,chr4,chr5,chr6,chr7,chr8,chr9,chr10,chr11

Number of alignments: 107168
Number of query sequences: 11

After filtering... Number of alignments: 1596
After filtering... Number of query sequences: 11

Warning: ‘timedatectl’ indicates the non-existent timezone name ‘n/a’
Warning: Your system is mis-configured: ‘/etc/localtime’ is not a symlink
Warning: It is strongly recommended to set envionment variable TZ to ‘Asia/Chongqing’ (or equivalent)

Looking forward to your reply
Thanks.
Lyu Qiang

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant