Tag Archives: samtools

How to select genomic region using samtools

To select a genomic region using samtools, you can use the faidx command. This command is used to index a FASTA file and extract subsequences from it. Here is an example of how to use faidx to select a genomic … Continue reading

Posted in bioinformatics | Tagged | Leave a comment

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