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

Questions on loading fastq.gz, and output.sam #111

Open
rl4940 opened this issue Jul 30, 2024 · 0 comments
Open

Questions on loading fastq.gz, and output.sam #111

rl4940 opened this issue Jul 30, 2024 · 0 comments

Comments

@rl4940
Copy link

rl4940 commented Jul 30, 2024

Hi Fritz team!

Thank you so much for developing this aligner. Currently I am using NGMLR along with minimap2, and I have 2 main questions regarding this:
below is my current code:

#### running ngmlr ####
ngmlr -t 20 -r $ref -q test_pipe1_standard/hg002head1.fastq.gz -o hg002.sam -x ont

### sorting and indexing the bam file ####
samtools view -b --threads 10 hg002.sam > hg002_out.bam
samtools sort hg002_out.bam -o hg002_out.sorted.bam
samtools index -M --threads 10 hg002_out.sorted.bam

Question 1: how to load multiple fastq.gz files into NGMLR? e.g. HG002 on GIAB has 3 fastq.gz files, so how to load them in?
Question 2: the output from NGMLR is SAM format, to eliminate unnecessary files, is there a way to re-set it into BAM output? if so, is there a way to auto-index the BAM file in NGMLR just like --write-index in minimap2?

Thanks!
R.L.

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