]0;/data/results/reference/mmu/mm9/mm9-ncRNAreczko@max:/data/results/reference/mmu/mm9/mm9-ncRNA$ export SHRIMP_FOLDER=/data/results/tools/align/SHRiMP_2_2_3
]0;/data/results/reference/mmu/mm9/mm9-ncRNAreczko@max:/data/results/reference/mmu/mm9/mm9-ncRNA$ python /data/results/tools/align/SHRiMP_2_2_3/SHRiMP_2_2_3/utils/project-db.py --seed 010111111111111,101011111111111,110101111111111,111010111111111,111101011111111,111110101111111,111111010111111,111111101011111,111111110101111,111111111010111,111111111101011,111111111110101,111111111111010,111111111111101,111111111111111  --h-flag --shrimp-mode cs mm9-ncRNA.fa

#max
./shrimp-map-ncRNA1.sh1  &> shrimp-map-ncRNA1.sh1.log &
#fix (avoid tmp bam file clashes)
cd b
./shrimp-map-ncRNA1.sh2  &> shrimp-map-ncRNA1.sh2.log &
mv b/* .


for i in *Mm.bam	 
do
    echo $i
    # add MD flags for bowtie format
    /data/results/tools/samtools/samtools-1.3/samtools calmd -b $i /data/results/reference/mmu/mm9/mRNA-genomic/mm9-genomic-mRNA.fa  >& /dev/null > foo.bam2
    samtools view -h foo.bam2  | awk -f /data/images/proton/DKlab/mr/parclip/paralyzer/sam2bowtie.awk  > $i.md.bt
    mv foo.bam2 $i.md.bam
done
cd /data/results/reference/mmu/mm9/ncRNA/
/home/reczko/bin/faToTwoBit mm9-ncRNA.fa mm9-ncRNA.2bit
#ini file:
/data/results/reference/mmu/mm9/ncRNA/mm9-ncRNA.2bit