Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Kaskere
    Junior Member
    • Jul 2018
    • 7

    #1

    Samtools is not working

    Hello,
    I wrote a script for obtaining bam files, however it looks like it is not understanding flags given. I wrote a flag for no header - bS but still keep getting errors. I already feel desperated

    #!/bin/bash -x

    module load bio/samtools/1.10

    export INPATH=/dir/subdir
    export INPATH1=/dir/subdir1
    export INPATH2=/dir/subdir2
    export INPATH3=/dir/subdir3
    export INPATH4=/dir/subdir4
    export INPATH5=/dir/subdir5
    export INPATH6=/dir/subdir6
    export INPATH7=/dir/subdir7
    export reference=/dir/subdir/refseq.fasta

    samtools faidx $reference

    for file in $INPATH/*.sam ;
    do
    bname=$(basename $file ".sam")
    echo $bname
    file=$INPATH/*.sam
    echo file: $file
    bamfile=$INPATH1/$bname".bam"
    namesort=$INPATH2/$bname"_sorted.bam"
    fixmate=$INPATH3/$bname"_sorted_fixmate.bam"
    positionsort=$INPATH4/$bname"_sorted_fixmate_position.bam"
    markdup=$INPATH5/$bname"_sorted_fixmate_position_markdup.bam"
    bai=$INPATH6/$bname".bai"
    stats=$INPATH7/$bname".txt"
    echo $bamfile $namesort $fixmate $positionsort $markdup

    samtools view -bT $reference $file > $bamfile
    samtools sort -@ 12 -n $bamfile > $namesort
    samtools fixmate -m $namesort > $fixmate
    samtools sort -@ 12 $fixmate > $positionsort
    samtools markdup $positionsort > $markdup
    samtools index > $markdup
    samtools flagstat $markdup > $stats
    done


    this is the error output

    ./samtools.sh: line 25: 348179 Segmentation fault samtools view -bT $reference $file > $bamfile
    + samtools sort -@ 12 -n file.bam
    samtools sort: failed to read header from "file.bam"
  • LucyaSmit
    Banned
    • Sep 2020
    • 1

    #2
    So strange that it happened

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
      by SEQadmin2



      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
      ...
      10-14-2026, 09:55 AM
    • SEQadmin2
      New Genomics Technologies Take Aim at Long-Standing Limits
      by SEQadmin2


      Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

      We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
      ...
      09-28-2026, 10:25 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 10-22-2026, 11:13 AM
    0 responses
    29 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-22-2026, 10:48 AM
    0 responses
    21 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    38 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    44 views
    0 reactions
    Last Post SEQadmin2  
    Working...