Unconfigured Ad

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

    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

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 06-09-2026, 11:58 AM
    0 responses
    28 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    35 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    40 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-02-2026, 12:03 PM
    0 responses
    62 views
    0 reactions
    Last Post SEQadmin2  
    Working...