the pos in sam file generated by minimap2 can not be correctly parsed. · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

the pos in sam file generated by minimap2 can not be correctly parsed.

Hello my friends, I encoutered a problem that when I used "minimap2 -a reference.fa reads.fa > mapping.sam" to generate a sam file, I found if the "pos" in the sam was used together with "CIGAR" to deduce the mapping, the result was wrongly mapped. When I ignored the "pos", the result is correct. anyone encountered similar problems? Thanks.

cigar minimap2 pos

• 685 views • link

written 4.1 years ago by Aaron •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

how to get a sequence from read fragments with ref sequence?
written 4.1 years ago by Aaron •

Hello my friends, I currently get some read fragments and a reference genome sequence. How can I get the whole sequence from these read fragments …
Nanopore data error
written 5.4 years ago by mailard •

Hello! I use ONT data I mapped reads with minimap2: `-a -x splice -k 14` and etc - next use samtools to sort But when …
Question regarding minimap2 presets
written 5.9 years ago by Rogerio Ribeiro •

Hi biostars. I have a question regarding minimap2. For my biological problem, I am doing differential expression analysis and differential transcript usage for six (3 …
calctruequality in bbmap
written 8.2 years ago by coyk •

I'm trying to recalibrate Q scores of a NextSeq run using MiSeq contigs assembled with Tadpole. The commands I used to map the reads to …
Fastq and sam number of reads
written 8.4 years ago by ceboral •

Hi! I'm having some little problems trying to figure out if the mapping of my datasets were correct. I have downloaded the datasets in SRA …
Difficulty in python while loop
written 8.7 years ago by KVC_bioinfo

Hello All, I wrote following python code: pos = 100 cigar = [(0, 150), (3, 50), (0, 130), (2, 270), (0, 750),(2,50),(3, 70), (0,500)] result=list() …
My VCF files are of zero length?
written 9.0 years ago by mahnoornaseer97 •

i have sorted bam files , i used these command to generate VCF file but my file only has headers, it does not show any …
Forum: Stranger Things: unexpected limitations of popular tools
written 9.1 years ago by Istvan Albert

Twice in just as many days, we were left scratching our heads as mainstream tools broke unexpectedly and spectacularly on what seemed to be normal …
Calculate percentage mapped reads from Bowtie2/RSEM
written 9.6 years ago by ando.kelli •

Hey all, I used the following script to align my reads to a reference transcriptome (bowtie2), and calculate transcript abundance (RSEM). $TRINITY_HOME/util/align_and_estimate_abundance.pl --transcripts location_of_data/assembly.fasta --seqType …
BWA MEM mapping issue
written 11.5 years ago by marina.v.yurieva

I'm mapping two samples independently (PE reads) using bwa-mem: ``` bwa mem -t 3 -M ref.fasta sample1_R1.fastq sample1_R2.fastq > sample1.sam # (2121026 mapped reads) bwa …