Examples:

1. dataset downloads

We use a real dataset from human genome(hg18) with 427786 reads with the length of 75bps. Half of them cover all known ref genes' junction sites, the other half are directly extracted from coding region.

The dataset could be download here (fasta format or fastqformat).

2. run commands

perl runABMapper.pl -ref reflist.in -input hg18ens.sp.fa -overlap 3 -seed 10 -multi 100 -output hg18ens.sp.abmapper

3. output explanation