f = FastaFile('HSP100.fa', indexed=True, interface='container')
print "Iteration"
i = 0
for h, s in f:
print h, s
i += 1
if i == 4: break
print '---'
for h, s in f:
print h, s
rule()
f.reset()
print "Slicing"
print " Indexed"
print f[0]
print
print f[2]
print
print f.readOne()
rule()
print " Slicing"
print f[5:8]
print
print f[0:2]
rule()
print "As mapping"
f = FastaFile('HSP100.fa', indexed=True, interface='container')
print "Iteration"
i = 0
for h,s in f:
print h,s
i += 1
if i==4: break
print '---'
for h,s in f:
print h,s
rule()
f.reset()
print "Slicing"
print " Indexed"
print f[0]
print
print f[2]
print
print f.readOne()
rule()
print " Slicing"
print f[5:8]
print
print f[0:2]
rule()
print "As mapping"
