def test_filter_on_TSS(self):
rna_seq = FeatureSet('Transcripts-ex1.txt')
lower,upper = 5000000,10000000
rna_tss = rna_seq.filterByTSS(upper,lower)
self.assertEqual(len(rna_tss),3,
"Wrong number of transcripts filtered on TSS")
for rna_data in rna_tss:
self.assertTrue((rna_data.getTSS() >= lower and
rna_data.getTSS() <= upper),
"Transcript outside range")
