Hello,
We would like to try and use iRep with some long read data we have. Our reads were assembled and binned, yielding several bins which are represented by just 1 large circular contig. All of the reads were mapped to all of the bins (concatenated multifasta file) using minimap2 generating a samfile, which was then sorted prior to using iRep.
I now ran iRep as follows for each bin (bin1 as an example here):
iRep -f bins/bin1.fa -s align-bins.sorted.sam -o d3-irep-bin1
Unfortunately it looks like iRep did not work, and the average coverage is significantly lower than expected. The bin1 (represented by a single large contig) actually has a mean depth of ~130, and no t0.11 as indicated in the iRep output PDF (see attached files).
I also tested running iRep with bin*.fa for analysis of all bins with the same result.
Hoping to get some input on this, as we're excited to use iRep on our long read datasets. Thanks.
d3-irep-bin1.pdf
Hello,
We would like to try and use iRep with some long read data we have. Our reads were assembled and binned, yielding several bins which are represented by just 1 large circular contig. All of the reads were mapped to all of the bins (concatenated multifasta file) using minimap2 generating a samfile, which was then sorted prior to using iRep.
I now ran iRep as follows for each bin (bin1 as an example here):
iRep -f bins/bin1.fa -s align-bins.sorted.sam -o d3-irep-bin1Unfortunately it looks like iRep did not work, and the average coverage is significantly lower than expected. The bin1 (represented by a single large contig) actually has a mean depth of ~130, and no t0.11 as indicated in the iRep output PDF (see attached files).
I also tested running iRep with bin*.fa for analysis of all bins with the same result.
Hoping to get some input on this, as we're excited to use iRep on our long read datasets. Thanks.
d3-irep-bin1.pdf
