minimap2 | versatile pairwise aligner | Genomics library
kandi X-RAY | minimap2 Summary
kandi X-RAY | minimap2 Summary
A versatile pairwise aligner for genomic and spliced nucleotide sequences
Support
Quality
Security
License
Reuse
Top functions reviewed by kandi - BETA
Currently covering the most popular Java, JavaScript and Python libraries. See a Sample of minimap2
minimap2 Key Features
minimap2 Examples and Code Snippets
for record in FASTQ.Reader(GzipDecompressorStream(open("data_file.fastq.gz")))
seq = string(sequence(record))
check = py"mappy"(seq,aligner)
end
close(reader)
Community Discussions
Trending Discussions on minimap2
QUESTION
I am attempting to merge x number of bam files produced via performing multiple alignments at once (on batches of y number of fastq files) into one single bam file in Nextflow.
So far I have the following when performing the alignment and sorting/indexing the resulting bam file:
...ANSWER
Answered 2021-Mar-04 at 15:53The simplest fix would probably to stop passing the static dirstring and runstring around via channels:
QUESTION
I am trying to process MinION cDNA amplicons using Porechop with Minimap2 and I am getting this error.
...ANSWER
Answered 2020-Sep-09 at 22:50First of all, you can always debug the problems like that specifying the flag --printshellcmds
. That would print all shell commands that Snakemake runs under the hood; you may try to run them manually and locate the problem.
As for why your rule doesn't produce any output, my guess is that samtools requires explicit filenames or -
to use stdin:
Samtools is designed to work on a stream. It regards an input file '-' as the standard input (stdin) and an output file '-' as the standard output (stdout). Several commands can thus be combined with Unix pipes. Samtools always output warning and error messages to the standard error output (stderr).
So try that:
QUESTION
I was trying to run a command which ideally looks like this,
minimap2 -a -x map-ont -t 20 /staging/reference.fasta fastq/sample01.fastq | samtools view -bS -F 4 - | samtools sort -o fastq_minon/sample01.bam
Similarly, I have multiple samples (referring to fastq/sample01.fastq) in the folder.
The snakemake file I wrote to automate this behaviour is, however, parsing all files at once in the command like,
minimap2 -a -x map-ont -t 1 /staging/reference.fasta fastq/sample02.fastq fastq/sample03.fastq fastq/sample01.fastq | samtools view -bS -F 4 - | samtools sort -o fastq_minon/sample02.bam fastq_minon/sample03.bam fastq_minon/sample01.bam
I have pasted the code and logs below. Please help me try to figure out this mistake.
Code
...ANSWER
Answered 2020-May-05 at 10:01The expand
function is used to create a list. Thus, in your rule minimap
, you're telling snakemake that you want all fastq files as input and that the rule will produce as many bam files. What you want is a rule that will be triggered for every sample using a wildcard:
Community Discussions, Code Snippets contain sources that include Stack Exchange Network
Vulnerabilities
No vulnerabilities reported
Install minimap2
Support
Reuse Trending Solutions
Find, review, and download reusable Libraries, Code Snippets, Cloud APIs from over 650 million Knowledge Items
Find more librariesStay Updated
Subscribe to our newsletter for trending solutions and developer bootcamps
Share this Page