abs-tudelft / abs-tudelft/gen-mpi

Implementation using BWA-mem2

Đang mở
#1 1 bình luận 0 reaction 0 người được giao Xem trên GitHub
Ngôn ngữ chính
C
Star
3
Fork
1
Chỉ số merge pull request
Không có pull request nào được merge trong 30 ngày

Mô tả

Hi,
I am interested in using your software for variant calling.
I looked through the code, but it seems to be for long-reads (no paired-end files).`
time srun -N 2 -n 2 /home/tahmad/hawk/new/pc/minimap2-chrms/minimap2 -t 128 **-ax map-ont -R** '@RG\tID:None\tSM:sample\tPL:Pacbio\tLB:sample\tPU:lane' /scratch-shared/tahmad/bio_data/reference/GRCh38_no_alt_analysis_set.fasta /scratch-shared/tahmad/bio_data/long/HG002/ONT/HG002.fastq > /scratch-shared/tahmad/bio_data/long/HG002/ONT/out.sam`

Is it possible to use bwa-mem2 instead of bwa-mem for faster mapping?

Could you provide an example of how to run bwa-mem2, or bwa-mem, with paired-end reads, or
can I provide the coordinate-sorted bam file?

Thanks;

Hướng dẫn đóng góp

Chưa lập chỉ mục được hướng dẫn đóng góp cho kho mã nguồn này

Hướng nghiên cứu

The issue asks about integrating bwa-mem2 for paired-end reads. First, examine the existing mapping code to see how bwa-mem is called, likely in shell scripts or C wrappers. Check if the pipeline accepts coordinate-sorted BAM files as input. Look for configuration files or command-line arguments that specify the aligner. Determine what changes are needed to support bwa-mem2 and paired-end FASTQ inputs, then test with a small dataset.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.

Đánh giá

Công nghệ
c, shell
Lĩnh vực
bioinformatics
Loại issue
Tính năng
Độ khó
4/5
Thời gian dự kiến
3-5 ngày
Mức độ hoạt động
Đình trệ
Độ rõ ràng
Khá rõ ràng
Mức phù hợp với người mới
35/100

Nhận issue mới trong hộp thư của bạn

Bản tóm tắt ngắn những issue GitHub phù hợp với người mới.