def test__eq__(self):
# Check equality of FeatureSets
feature_set1 = FeatureSet()
feature_set2 = FeatureSet()
# Empty feature sets
self.assertEqual(feature_set1,feature_set2)
# Populate
feature_set1.addFeature(Feature('CG1000','chr1','1','2','+'))
feature_set2.addFeature(Feature('CG1000','chr1','1','2','+'))
self.assertEqual(feature_set1,feature_set2)
# Add second
feature_set1.addFeature(Feature('CG2000','chr1','1','2','+'))
self.assertNotEqual(feature_set1,feature_set2)
feature_set2.addFeature(Feature('CG2000','chr1','1','2','+'))
self.assertEqual(feature_set1,feature_set2)
# Add third
feature_set1.addFeature(Feature('CG2001','chr2',3,4,'-'))
feature_set2.addFeature(Feature('CG2002','chr2',3,5,'+'))
self.assertNotEqual(feature_set1,feature_set2)
