![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wiki is it working? | Rafael_RPS | Wiki Discussion | 1 | 07-01-2019 07:38 AM |
samtools view region not working properly | ianh | Bioinformatics | 1 | 11-04-2013 02:17 AM |
samtools mpileup -r not working | tracecakes | Bioinformatics | 3 | 08-14-2013 06:55 AM |
samtools: using -l in mpileup not working as expected | hmviit | Bioinformatics | 0 | 01-12-2012 07:08 AM |
Working Of Abyss | aarifjindani | Bioinformatics | 0 | 11-06-2011 10:05 PM |
![]() |
|
Thread Tools |
![]() |
#1 |
Junior Member
Location: Europe Join Date: Jul 2018
Posts: 6
|
![]()
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" |
![]() |
![]() |
![]() |
#2 |
Banned
Location: Milwaukie Join Date: Sep 2020
Posts: 1
|
![]()
So strange that it happened
|
![]() |
![]() |
![]() |
Thread Tools | |
|
|