I,m trying to annotate the genome of a newly assembled ant. I am following your tutorial here:http://bioinf.uni-greifswald.de/bioinf/ ... s.Augustus, and I'm trying to obtain hints for gene prediction in augustus. I'm having some problem in getting Exon-exon junction hints. After aligning read to the genome scaffolds (a test of 10M reads) I run filterBAM as described in the tutorial:
Code: Select all
filterBam --uniq --paired --in Aligned.out.bam --out Aligned.out.sf.bam
Code: Select all
processed line 300001
-------------------------------------
Summary of filtered alignments:
-------------------------------------
unmapped : 0
percent identity: 561
coverage : 5372
not paired : 358611
quantiles of unspliced insert lengths:
[Quantiles]:vector::_M_range_check; insertlen.size()=0
unique : 0
-------------------------------------
Cmd line:
filterBam --uniq --paired --in Aligned.out.bam --out Aligned.out.sf.bam
-------------------------------------
Elapsed time: 0.000 seconds.
Thank you a lot
Francesco