[9b26b7]: / third_party / nucleus / io / fastq_writer_test.cc

Download this file

128 lines (113 with data), 5.1 kB

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
/*
* Copyright 2018 Google LLC.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
* are met:
*
* 1. Redistributions of source code must retain the above copyright notice,
* this list of conditions and the following disclaimer.
*
* 2. Redistributions in binary form must reproduce the above copyright
* notice, this list of conditions and the following disclaimer in the
* documentation and/or other materials provided with the distribution.
*
* 3. Neither the name of the copyright holder nor the names of its
* contributors may be used to endorse or promote products derived from this
* software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
* AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
* ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
* LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
* CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
* SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
* INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
* CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
* ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
* POSSIBILITY OF SUCH DAMAGE.
*
*/
#include "third_party/nucleus/io/fastq_writer.h"
#include <utility>
#include <vector>
#include <gmock/gmock-generated-matchers.h>
#include <gmock/gmock-matchers.h>
#include <gmock/gmock-more-matchers.h>
#include "tensorflow/core/platform/test.h"
#include "third_party/nucleus/testing/test_utils.h"
#include "third_party/nucleus/util/utils.h"
#include "third_party/nucleus/core/status_matchers.h"
#include "tensorflow/core/lib/core/status.h"
#include "tensorflow/core/platform/env.h"
namespace nucleus {
using std::vector;
void CreateRecord(const string& id, const string& description,
const string& sequence, const string& quality,
nucleus::genomics::v1::FastqRecord* record) {
record->set_id(id);
record->set_description(description);
record->set_sequence(sequence);
record->set_quality(quality);
}
class FastqWriterTest : public ::testing::Test {
protected:
void SetUp() override {
nucleus::genomics::v1::FastqRecord first, second, third;
CreateRecord("NODESC:header", "", "GATTACA", "BB>B@FA", &first);
CreateRecord(
"M01321:49:000000000-A6HWP:1:1101:17009:2216", "1:N:0:1",
"CGTTAGCGCAGGGGGCATCTTCACACTGGTGACAGGTAACCGCCGTAGTAAAGGTTCCGCCTTTCACT",
"AAAAABF@BBBDGGGG?FFGFGHBFBFBFABBBHGGGFHHCEFGGGGG?FGFFHEDG3EFGGGHEGHG",
&second);
CreateRecord("FASTQ", "contains multiple spaces in description",
"CGGCTGGTCAGGCTGACATCGCCGCCGGCCTGCAGCGAGCCGCTGC",
"FAFAF;F/9;.:/;999B/9A.DFFF;-->.AAB/FC;9-@-=;=.", &third);
golden_ = {first, second, third};
}
vector<nucleus::genomics::v1::FastqRecord> golden_;
};
TEST_F(FastqWriterTest, WritingWorks) {
const string output_filename = MakeTempFile("writes_fastq.fastq");
std::unique_ptr<FastqWriter> writer =
std::move(FastqWriter::ToFile(output_filename,
nucleus::genomics::v1::FastqWriterOptions())
.ValueOrDie());
for (const nucleus::genomics::v1::FastqRecord& record : golden_) {
ASSERT_THAT(writer->Write(record), IsOK());
}
ASSERT_THAT(writer->Close(), IsOK());
writer.reset();
string contents;
TF_CHECK_OK(tensorflow::ReadFileToString(tensorflow::Env::Default(),
output_filename, &contents));
const string kExpectedFastqContent =
"@NODESC:header\n"
"GATTACA\n"
"+\n"
"BB>B@FA\n"
"@M01321:49:000000000-A6HWP:1:1101:17009:2216 1:N:0:1\n"
"CGTTAGCGCAGGGGGCATCTTCACACTGGTGACAGGTAACCGCCGTAGTAAAGGTTCCGCCTTTCACT\n"
"+\n"
"AAAAABF@BBBDGGGG?FFGFGHBFBFBFABBBHGGGFHHCEFGGGGG?FGFFHEDG3EFGGGHEGHG\n"
"@FASTQ contains multiple spaces in description\n"
"CGGCTGGTCAGGCTGACATCGCCGCCGGCCTGCAGCGAGCCGCTGC\n"
"+\n"
"FAFAF;F/9;.:/;999B/9A.DFFF;-->.AAB/FC;9-@-=;=.\n";
EXPECT_EQ(kExpectedFastqContent, contents);
}
TEST_F(FastqWriterTest, WritesGzippedFiles) {
const string output_filename = MakeTempFile("writes_fastq.fastq.gz");
std::unique_ptr<FastqWriter> writer =
std::move(FastqWriter::ToFile(output_filename,
nucleus::genomics::v1::FastqWriterOptions())
.ValueOrDie());
ASSERT_THAT(writer->Close(), IsOK());
string contents;
TF_CHECK_OK(tensorflow::ReadFileToString(tensorflow::Env::Default(),
output_filename, &contents));
EXPECT_THAT(IsGzipped(contents),
"FASTQ writer should be able to writed gzipped output");
}
} // namespace nucleus