Samtools view error: header include @SQ line missing LN tag.

After running minimap2 aligning .fastq to index .mmi creating .bam and pipeing the result into samtools view the above error ensues. What can I do abou it?

minimap2

• 1,362 views • link

updated 3.4 years ago by Ram

• written 5.3 years ago by Digimon •

What is the exact command you used?

• 0 views • link

written 5.3 years ago by GenoMax

minimap2 -t 40 -ax splice -k14 -u f --secondary=no input.mmi othr_input.fastq.gz | samtools view -Sb > othr_input.bam

• 0 views • link

written 5.3 years ago by Digimon •

0 answers

No answers yet.

More posts like this

HLA alignment, approximation and frequencies
written 1.7 years ago by oumo •

Hello, I have tried aligning this fastq file (long reads from ONT sequencing) with minimap2. I came up with a sam file, converted to bam …
How to use minimap2?
written 4.1 years ago by MobiusT •

Hi, I am new to bioinformatics and I am tasked to do a read alignment. I have fastq.gz files, and I should run minimap2 over …
IGV Sashimi plot from ONT long read sequencing of cDNA amplicons
written 4.2 years ago by WRUZIC1 •

Hi - I'm trying to generate a sashimi plot using data from Oxford Nanopore sequencing of cDNA amplicons. I've aligned my fastq files to the …
find function with samtools - convert format to all subfolders
written 4.3 years ago by Korsocius

I would like to find all .bam files in subfolders and exec it with samtools. Concretely conversion of bam to cram to same folder and …
Samtools SORT or bedtools bamToFastq giving error when converting reads to fastq.
written 5.5 years ago by DNAngel

I am trying to convert mapped reads and unmapped reads into fastq format but I am having trouble with the last bit of my pipeline. …
Error when trying to convert BAM file to FASTQ.
written 5.5 years ago by DNAngel

I thought I was finally in the stretch but lo and behold a new issue arises. Here is my whole pipeline that I used: bwa …
bedtools questions ! can't sort/index by samtools after bedToBam command
written 5.8 years ago by zhaizhaoyu7 •

Hi all ! After I use **STAR** to alignment , I got the **BAM** file , this **BAM** file could index successfully by samtools , …
Why Samtools view don't extract sequence out of the given region?
written 6.3 years ago by kamanovae •

Hi! I have nanopore sequencing results in fastq format. I need to know the coverage of specific genes. To do this, I aligned the sequence …
getting error when indexing bam files
written 7.1 years ago by Sara

I have aligned my fastq files to hg19 using this command: bwa mem hg19.fa /DATA/myfile.fastq.gz then I made bam files from sam files using: samtools …
samtools sorting and indexing
written 9.1 years ago by ggman

Hi friends, I am attempting to sort my bam files that I obtained from my bowtie sam files. I am not indexing them appropriate according …