for i in *4p10.sam
do
    echo $i
    cat $i |samtools view -bS -t /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/WholeGenomeFasta/genome.fa.fai -  | samtools sort -o - - > $i.bam
    samtools index $i.bam
    /data/results/tools/align/rnaseqmut/rnaseqmut-master/bin/rnaseqmut.linux.x64 -r /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/WholeGenomeFasta/genome.fa -m 0 $i.bam  | awk -f /data/results/tools/align/rnaseqmut/count-conversions1.awk > $i.bam.conversion.stats
    cat !$
done
