Tag Archives: samtools

missing SAM header with minimap2 and samtools

When using minimap2 to map sequencing reads onto a reference, you can use this kind of command (be careful, this is wrong as you will see later): minimap2 -a -x map-pb test.fastq reference.fasta > minimap.samminimap2 -a -x map-pb test.fastq reference.fasta … Continue reading

Posted in bioinformatics | Tagged , | Leave a comment

failed to convert from cram to bam (parse error CIGAR character)

In order to convert a bioinformatic file from CRAM to BAM format, I naively used the samtools command available on a cluster but got this error: samtools view -T reference.fa -b -o output.bam input.cram [sam_header_read2] 3366 sequences loaded. [sam_read1] reference … Continue reading

Posted in bioinformatics | Tagged | Leave a comment