* pq flag replaced with rq
* ReadCountMemFail no longer reported
* Update test data
This commit is contained in:
Nigel Delaney 2015-10-14 11:48:10 -07:00
Родитель db7e581c28
Коммит c4c1640c9c
3 изменённых файлов: 8 добавлений и 13 удалений

Просмотреть файл

@ -53,27 +53,19 @@ namespace Bio.IO.PacBio
get { return statusCounts [1]; }
}
/// <summary>
/// Count of subreads not added to consensus for allocating too much memory.
/// </summary>
/// <value>The ReadCount mem fail.</value>
public int ReadCountMemFail {
get { return statusCounts [2]; }
}
/// <summary>
/// Count of subreads not added to consensus for having too low a Z-score.
/// </summary>
/// <value>The ReadCount bad zscore.</value>
public int ReadCountBadZscore {
get { return statusCounts [3]; }
get { return statusCounts [2]; }
}
/// <summary>
/// Gets the read count other.
/// </summary>
/// <value>The read count other.</value>
public int ReadCountOther {
get{ return statusCounts [4]; }
get{ return statusCounts [3]; }
}
/// <summary>
@ -141,8 +133,11 @@ namespace Bio.IO.PacBio
} else if (v.Tag == "zm") {
HoleNumber = (int)Convert.ToInt32 (v.Value);
} else if (v.Tag == "pq") {
ReadQuality = Convert.ToSingle(v.Value);
} else if (v.Tag == "za") {
// This tag is now deprecated by the rq tag
ReadQuality = Convert.ToSingle (v.Value);
} else if (v.Tag == "rq") {
ReadQuality = Convert.ToSingle (v.Value);
}else if (v.Tag == "za") {
AvgZscore = (float)Convert.ToSingle (v.Value);
} else if (v.Tag == "rs") {
statusCounts = v.Value.Split (',').Skip (1).Select (x => Convert.ToInt32 (x)).ToArray ();

Просмотреть файл

@ -16,7 +16,7 @@ namespace Bio.Tests.PacBio
{
string fname = System.IO.Path.Combine ("TestUtils", "PacBio", "ccs.bam");
var csp = new PacBioCCSBamReader ();
var seqs = csp.Parse (fname).ToList();
var seqs = csp.Parse (fname).Select(z => z as PacBioCCSRead).ToList();
var seq4 = seqs [4];
Assert.AreEqual (146331, seq4.HoleNumber);

Двоичные данные
src/Source/Tests/TestData/TestUtils/PacBio/ccs.bam Executable file → Normal file

Двоичный файл не отображается.