pdb-tools 2.5.0__py3-none-any.whl → 2.6.0__py3-none-any.whl

tests/test_pdb_selaltloc.py

@@ -1,943 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-#
-# Copyright 2018 João Pedro Rodrigues
-#
-# Licensed under the Apache License, Version 2.0 (the "License");
-# you may not use this file except in compliance with the License.
-# You may obtain a copy of the License at
-#
-#    http://www.apache.org/licenses/LICENSE-2.0
-#
-# Unless required by applicable law or agreed to in writing, software
-# distributed under the License is distributed on an "AS IS" BASIS,
-# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-# See the License for the specific language governing permissions and
-# limitations under the License.
-
-"""
-Unit Tests for `pdb_selaltloc`.
-"""
-import os
-import sys
-import unittest
-
-from config import data_dir
-from utils import OutputCapture
-
-
-class TestTool(unittest.TestCase):
-    """
-    Generic class for testing tools.
-    """
-
-    def setUp(self):
-        # Dynamically import the module
-        name = 'pdbtools.pdb_selaltloc'
-        self.module = __import__(name, fromlist=[''])
-
-    def exec_module(self):
-        """
-        Execs module.
-        """
-
-        with OutputCapture() as output:
-            try:
-                self.module.main()
-            except SystemExit as e:
-                self.retcode = e.code
-
-        self.stdout = output.stdout
-        self.stderr = output.stderr
-
-        return
-
-    def test_is_same_group_1(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we are in the SAME altloc group.
-
-        previous line: ""  # first line of the PDB file
-        current  line: " ALA    1"
-
-        The result of `is_another_altloc_group` should be False.
-        """
-        result = self.module.is_another_altloc_group(
-            ' ', '', '1', '', 'ALA', '', {}, {})
-
-        self.assertFalse(result)
-
-    def test_is_same_group_2(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we are in the SAME altloc group.
-
-        All `None` case.
-
-        The result of `is_another_altloc_group` should be False.
-        """
-        result = self.module.is_another_altloc_group(
-            ' ', None, None, None, None, None, {}, {})
-
-        self.assertFalse(result)
-
-    def test_is_same_group_3(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we are in the SAME altloc group.
-
-        Example: all parameters are the same as the previous line.
-
-        previous line: "BALA   12"
-        current  line: "BALA   12"
-
-        The result of `is_another_altloc_group` should be False.
-        """
-        result = self.module.is_another_altloc_group(
-            'B', 'B', '12', '12', 'ALA', 'ALA', {'B': None},
-            {'B': {('ALA', '12')}}
-            )
-        self.assertFalse(result)
-
-
-    def test_is_same_group_4(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we are in the SAME altloc group.
-
-        Example: Multiple residue altloc. This considers altloc spanning
-        several residues. See example dummy_altloc2.pdb.
-
-        previous line: "AGLU   25"
-        current  line: "ALEU   26"
-
-        It is important to note also the dictionary input.
-
-        The result of `is_another_altloc_group` should be False.
-        """
-        result = self.module.is_another_altloc_group(
-            'A', 'A', '26', '25', 'LEU', 'GLU', {'A': ['lines']},
-            {'A': {('GLU', '25')}}
-            )
-        self.assertFalse(result)
-
-    def test_is_same_group_5(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we are in the SAME altloc group.
-
-        Example: Multiple residue altloc. This considers altloc spanning
-        several residues. See example dummy_altloc2.pdb.
-
-        previous line: " GLU   25"
-        current  line: "AALA   25"
-
-        It is important to note also the dictionary input.
-
-        The result of `is_another_altloc_group` should be False.
-        """
-        result = self.module.is_another_altloc_group(
-            'A', ' ', '25', '25', 'ALA', 'GLU', {' ': ['lines']},
-            {' ': {('GLU', '25')}}
-            )
-        self.assertFalse(result)
-
-    def test_is_another_group_1(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we entered another altloc group.
-
-        previous line: "BPRO    1"
-        current  line: " ALA    2"
-
-        The result of `is_another_altloc_group` should be True.
-        """
-        result = self.module.is_another_altloc_group(
-            ' ', 'B', '2', '1', 'ALA', 'PRO', {'B': ['lines']},
-            {'B': {('PRO', '1')}}
-            )
-        self.assertTrue(result)
-
-    def test_is_another_group_2(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we entered another altloc group.
-
-        previous line: " ALA    1"
-        current  line: " ALA    2"
-
-        The result of `is_another_altloc_group` should be True.
-        """
-        result = self.module.is_another_altloc_group(
-            ' ', ' ', '2', '1', 'ALA', 'ALA', {' ': ['lines']},
-            {' ': {('ALA', '1')}}
-            )
-        self.assertTrue(result)
-
-    def test_is_another_group_3(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we entered another altloc group.
-
-        previous line: " GLU    1"
-        current  line: " ALA    1"
-
-        The result of `is_another_altloc_group` should be True.
-        """
-        result = self.module.is_another_altloc_group(
-            ' ', ' ', '1', '1', 'ALA', 'GLU', {' ': ['lines']},
-            {' ': {('GLU', '1')}}
-            )
-        self.assertTrue(result)
-
-    def test_is_another_group_4(self):
-        """
-        Test function to identify we entered another altloc group.
-
-        This indicates we entered another altloc group.
-
-        previous line: "AGLU   25"
-        current  line: "ALEU   26"
-
-        The result of `is_another_altloc_group` should be True.
-        """
-        result = self.module.is_another_altloc_group(
-            'A', 'A', '26', '25', 'LEU', 'GLU', {' ': ['lines'], 'A': ['lines']},
-            {' ': {('LEU', '25')}, 'A': {('GLU', '26')}}
-            )
-        self.assertTrue(result)
-
-    def test_all_same_residue(self):
-        """Test all same residue."""
-        inp = {
-            ' ': [
-                    "ATOM      3  N   ASN A   1      22.066  40.557   0.420  1.00  0.00           N  ",
-                    "ATOM      3  H   ASN A   1      21.629  41.305  -0.098  1.00  0.00           H  ",
-                    "ATOM      3  H2  ASN A   1      23.236  40.798   0.369  1.00  0.00           H  ",
-                    "ATOM      3  H3  ASN A   1      21.866  40.736   1.590  1.00  0.00           H  ",
-                    ],
-            'B': ["ATOM      3  CA BASN A   1      20.000  30.000   0.005  0.60  0.00           C  "],
-            'A': ["ATOM      3  CA AASN A   1      21.411  39.311   0.054  0.40  0.00           C  "],
-            }
-
-        result = self.module.all_same_residue(inp)
-        self.assertTrue(result)
-
-    def test_all_same_residue_false(self):
-        """Test all same residue."""
-        inp = {
-            'B': ["ATOM      3  CA BSER A   2      20.000  30.000   0.005  0.60  0.00           C  "], 'A': ["ATOM      3  CA AASN A   1      21.411  39.311   0.054  0.40  0.00           C  "], } 
-        result = self.module.all_same_residue(inp)
-        self.assertFalse(result)
-
-    def test_partial_altloc(self):
-        inp = {
-            'A': [
-                    "ATOM    333  CA AGLU A  26     -10.000  -3.000 -12.000  0.50  4.89           C  ",
-                    "ANISOU  333  CA AGLU A  26      576    620    663     31     42     45       C  ",
-                 ],
-            'B':[
-                "ATOM    333  CA CGLU A  26     -10.679  -3.437 -12.387  1.00  4.89           C  ",
-                "ANISOU  333  CA CGLU A  26      576    620    663     31     42     45       C  ",
-                ],
-            }
-
-        result = self.module.partial_altloc(inp)
-        self.assertFalse(result)
-
-
-    def test_default(self):
-        """$ pdb_selaltloc data/dummy_altloc.pdb"""
-
-        # Simulate input
-        # pdb_selaltloc dummy.pdb
-        sys.argv = ['', os.path.join(data_dir, 'dummy_altloc.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 36)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['SER', 'GLU', 'ALA', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 2
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-        # Make sure we picked the right GLU atom
-        atom_CA_GLU_26 = [
-            l[30:54] for l in self.stdout
-            if l[12:16] == " CA " and l[17:20] == "GLU" and l[22:26] == "  26"
-        ]
-        self.assertEqual(len(atom_CA_GLU_26), 2)
-        self.assertEqual(atom_CA_GLU_26[0], " -10.679  -3.437 -12.387")
-
-    def test_default2(self):
-        """$ pdb_selaltloc data/dummy_altloc2.pdb"""
-
-        # Simulate input
-        # pdb_selaltloc dummy.pdb
-        sys.argv = ['', os.path.join(data_dir, 'dummy_altloc2.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 37)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['SER', 'GLU', 'ALA', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 2
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-        # Make sure we picked the right GLU atom
-        atom_CA_GLU_26 = [
-            l[30:54] for l in self.stdout
-            if l[12:16] == " CA " and l[17:20] == "GLU" and l[22:26] == "  26"
-        ]
-        self.assertEqual(len(atom_CA_GLU_26), 2)
-        self.assertEqual(atom_CA_GLU_26[0], " -10.679  -3.437 -12.387")
-
-    def test_select_loc_A(self):
-        """$ pdb_selaltloc -A data/dummy_altloc.pdb"""
-
-        # Simulate input
-        sys.argv = ['', '-A', os.path.join(data_dir, 'dummy_altloc.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 36)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['SER', 'GLU', 'ALA', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 2
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-        # Make sure we picked the right GLU atom
-        atom_CA_GLU_26 = [
-            l[30:54] for l in self.stdout
-            if l[12:16] == " CA " and l[17:20] == "GLU" and l[22:26] == "  26"
-        ]
-        self.assertEqual(len(atom_CA_GLU_26), 2)
-        self.assertEqual(atom_CA_GLU_26[0], " -10.000  -3.000 -12.000")
-
-    def test_select_loc_A_2(self):
-        """$ pdb_selaltloc -A data/dummy_altloc2.pdb"""
-
-        # Simulate input
-        sys.argv = ['', '-A', os.path.join(data_dir, 'dummy_altloc2.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 37)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['SER', 'GLU', 'ALA', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 2
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-        # Make sure we picked the right GLU atom
-        atom_CA_GLU_26 = [
-            l[30:54] for l in self.stdout
-            if l[12:16] == " CA " and l[17:20] == "GLU" and l[22:26] == "  26"
-        ]
-        self.assertEqual(len(atom_CA_GLU_26), 2)
-        self.assertEqual(atom_CA_GLU_26[0], " -10.000  -3.000 -12.000")
-
-    def test_select_loc_B(self):
-        """$ pdb_selaltloc -B data/dummy_altloc.pdb"""
-
-        # Simulate input
-        sys.argv = ['', '-B', os.path.join(data_dir, 'dummy_altloc.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 46)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['PRO', 'GLU', 'ALA', 'ILE', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 4
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-    def test_select_loc_B_2(self):
-        """$ pdb_selaltloc -B data/dummy_altloc2.pdb"""
-
-        # Simulate input
-        sys.argv = ['', '-B', os.path.join(data_dir, 'dummy_altloc2.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 47)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['PRO', 'GLU', 'ALA', 'ILE', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 4
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-    def test_select_loc_C(self):
-        """$ pdb_selaltloc -C data/dummy_altloc.pdb"""
-
-        # Simulate input
-        sys.argv = ['', '-C', os.path.join(data_dir, 'dummy_altloc.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 52)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['SER', 'PRO', 'GLU', 'ALA', 'ILE', 'LEU'] for _ in range(8)
-        ] + ["GLU"] * 2
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-    def test_select_loc_C_2(self):
-        """$ pdb_selaltloc -C data/dummy_altloc2.pdb"""
-
-        # Simulate input
-        sys.argv = ['', '-C', os.path.join(data_dir, 'dummy_altloc2.pdb')]
-
-        # Execute the script
-        self.exec_module()
-
-        # Validate results
-        self.assertEqual(self.retcode, 0)  # ensure the program exited OK.
-        self.assertEqual(len(self.stdout), 61)
-        self.assertEqual(len(self.stderr), 0)  # no errors
-
-        # Make sure sequence is correct
-        expected = [
-            aa for aa in ['SER', 'GLU', 'PRO', 'GLU', 'ALA', 'ILE', 'LEU']
-            for _ in range(8)
-            ] + ["GLU"] * 2
-        observed = [
-            ln[17:20] for ln in self.stdout if ln.startswith(("ATOM", "ANISOU"))
-        ]
-
-        self.assertEqual(observed, expected)
-
-    def test_gives_same_dummy_A(self):
-        """Test dummy.pdb is not altered because there are not altlocs."""
-        sys.argv = ['', '-A', os.path.join(data_dir, 'dummy.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 203)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout[80],
-            "ATOM      3  CA  ASN A   1      21.411  39.311   0.054  0.40  0.00           C  ")
-
-    def test_gives_same_dummy_B(self):
-        """Test dummy.pdb is not altered because there are not altlocs."""
-        sys.argv = ['', '-B', os.path.join(data_dir, 'dummy.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 203)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout[80],
-            "ATOM      3  CA  ASN A   1      20.000  30.000   0.005  0.60  0.00           C  "
-            )
-
-    def test_gives_same_dummy_maxocc(self):
-        """Test dummy.pdb is not altered because there are not altlocs."""
-        sys.argv = ['', os.path.join(data_dir, 'dummy.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 203)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout[80],
-            "ATOM      3  CA  ASN A   1      20.000  30.000   0.005  0.60  0.00           C  "
-            )
-
-    def test_vu7_maxocc(self):
-        """Test vu7.pdb properly selects highest altloc."""
-        sys.argv = ['', os.path.join(data_dir, 'vu7.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 30)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout[15],
-            "HETATM 2910  C27 VU7 A 403     -20.472  24.444  21.209  0.70 25.90           C  ",
-            )
-        self.assertEqual(
-            self.stdout[23:30],
-            [
-                "HETATM 2918  C20 VU7 A 403     -25.101  22.166  21.562  1.00 19.10           C  ",
-                "HETATM 2919  C22 VU7 A 403     -25.005  23.859  23.317  1.00 20.89           C  ",
-                "HETATM 2921  C25 VU7 A 403     -21.980  23.567  22.880  0.70 25.16           C  ",
-                "HETATM 2923  C26 VU7 A 403     -20.710  24.002  22.511  0.70 29.65           C  ",
-                "HETATM 2925  C28 VU7 A 403     -21.503  24.451  20.265  0.70 24.53           C  ",
-                "HETATM 2927  C29 VU7 A 403     -22.775  24.020  20.615  0.70 24.22           C  ",
-                "HETATM 2929  F30 VU7 A 403     -23.749  24.031  19.675  0.70 26.05           F  ",
-                ]
-            )
-
-    def test_vu7_maxocc_B(self):
-        """Test vu7.pdb properly selects altloc B."""
-        sys.argv = ['', '-B', os.path.join(data_dir, 'vu7.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 30)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout[15],
-            "HETATM 2910  C27 VU7 A 403     -20.472  24.444  21.209  0.70 25.90           C  ",
-            )
-        self.assertEqual(
-            self.stdout[23:30],
-            [
-                "HETATM 2918  C20 VU7 A 403     -25.101  22.166  21.562  1.00 19.10           C  ",
-                "HETATM 2919  C22 VU7 A 403     -25.005  23.859  23.317  1.00 20.89           C  ",
-                "HETATM 2921  C25 VU7 A 403     -21.980  23.567  22.880  0.70 25.16           C  ",
-                "HETATM 2923  C26 VU7 A 403     -20.710  24.002  22.511  0.70 29.65           C  ",
-                "HETATM 2925  C28 VU7 A 403     -21.503  24.451  20.265  0.70 24.53           C  ",
-                "HETATM 2927  C29 VU7 A 403     -22.775  24.020  20.615  0.70 24.22           C  ",
-                "HETATM 2929  F30 VU7 A 403     -23.749  24.031  19.675  0.70 26.05           F  ",
-                ]
-            )
-
-    def test_vu7_maxocc_A(self):
-        """Test vu7.pdb properly selects altloc A."""
-        sys.argv = ['', '-A', os.path.join(data_dir, 'vu7.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 30)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout[15],
-            "HETATM 2909  C27 VU7 A 403     -20.473  24.388  21.073  0.30 23.87           C  ",
-            )
-        self.assertEqual(
-            self.stdout[23:30],
-            [
-                "HETATM 2918  C20 VU7 A 403     -25.101  22.166  21.562  1.00 19.10           C  ",
-                "HETATM 2919  C22 VU7 A 403     -25.005  23.859  23.317  1.00 20.89           C  ",
-                "HETATM 2920  C25 VU7 A 403     -22.814  24.071  20.642  0.30 22.86           C  ",
-                "HETATM 2922  C26 VU7 A 403     -21.562  24.484  20.224  0.30 22.96           C  ",
-                "HETATM 2924  C28 VU7 A 403     -20.656  23.873  22.343  0.30 25.28           C  ",
-                "HETATM 2926  C29 VU7 A 403     -21.907  23.454  22.776  0.30 23.59           C  ",
-                "HETATM 2928  F30 VU7 A 403     -22.039  22.963  24.024  0.30 24.64           F  ",
-                ]
-            )
-
-    def test_anisou_lines(self):
-        """
-        Test anisou.pdb is not altered because there are not altlocs.
-
-        anisou.pdb has ANISOU lines. Ensure bug reported in #130 is
-        corrected.
-        """
-        infile = os.path.join(data_dir, 'anisou.pdb')
-        sys.argv = ['', infile]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 24)
-        self.assertEqual(len(self.stderr), 0)
-        with open(infile, "r") as fin:
-            expected_lines = [l.strip(os.linesep) for l in fin.readlines()]
-        self.assertEqual(self.stdout, expected_lines)
-
-    def test_anisou_with_altloc_maxocc(self):
-        """
-        Test anisou_altloc.pdb properly selects the highest altloc.
-
-        anisou.pdb has ANISOU lines. Ensure bug reported in #130 is
-        corrected.
-        """
-        infile = os.path.join(data_dir, 'anisou_altloc.pdb')
-        sys.argv = ['', infile]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 24)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout,
-            [
-                "ATOM      1  N   ALA A  31     -12.806   6.423  23.735  0.60 70.81           N  ",
-                "ANISOU    1  N   ALA A  31     7836   7867  11203   2151   -675    -66       N  ",
-                "ATOM      2  CA  ALA A  31     -13.433   7.788  23.746  0.60 65.66           C  ",
-                "ANISOU    2  CA  ALA A  31     7296   7485  10167   1829   -642    -81       C  ",
-                "ATOM      3  C   ALA A  31     -12.448   8.891  24.124  1.00 63.44           C  ",
-                "ANISOU    3  C   ALA A  31     6818   7632   9654   1744   -656     20       C  ",
-                "ATOM      4  O   ALA A  31     -11.549   8.680  24.937  1.00 65.52           O  ",
-                "ANISOU    4  O   ALA A  31     6891   8010   9994   1853   -759    238       O  ",
-                "ATOM      5  CB  ALA A  31     -14.628   7.834  24.691  1.00 66.76           C  ",
-                "ANISOU    5  CB  ALA A  31     7636   7487  10242   1630   -748    143       C  ",
-                "ATOM      6  N   THR A  32     -12.659  10.075  23.550  1.00 58.59           N  ",
-                "ANISOU    6  N   THR A  32     6249   7244   8768   1532   -569   -122       N  ",
-                "ATOM      7  CA  THR A  32     -11.806  11.232  23.788  1.00 55.42           C  ",
-                "ANISOU    7  CA  THR A  32     5678   7217   8161   1398   -585    -53       C  ",
-                "ATOM      8  C   THR A  32     -12.119  11.857  25.148  1.00 52.59           C  ",
-                "ANISOU    8  C   THR A  32     5384   6947   7651   1206   -743    183       C  ",
-                "ATOM      9  O   THR A  32     -13.166  11.585  25.770  1.00 49.84           O  ",
-                "ANISOU    9  O   THR A  32     5232   6412   7293   1142   -803    278       O  ",
-                "ATOM     10  CB  THR A  32     -11.965  12.321  22.696  0.55 51.87           C  ",
-                "ANISOU   10  CB  THR A  32     5260   6948   7500   1220   -449   -254       C  ",
-                "ATOM     11  OG1 THR A  32     -13.294  12.844  22.720  1.00 48.11           O  ",
-                "ANISOU   11  OG1 THR A  32     5047   6344   6890   1015   -452   -282       O  ",
-                "ATOM     12  CG2 THR A  32     -11.660  11.781  21.303  1.00 54.23           C  ",
-                "ANISOU   12  CG2 THR A  32     5483   7248   7875   1398   -280   -508       C  ",
-                ]
-            )
-
-    def test_anisou_with_altloc_maxocc_A(self):
-        """
-        Test anisou_altloc.pdb selects altloc -A.
-
-        anisou.pdb has ANISOU lines. Ensure bug reported in #130 is
-        corrected.
-        """
-        infile = os.path.join(data_dir, 'anisou_altloc.pdb')
-        sys.argv = ['', '-A', infile]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 24)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout,
-            [
-                "ATOM      1  N   ALA A  31     -12.806   6.423  23.735  0.40 70.81           N  ",
-                "ANISOU    1  N   ALA A  31     7836   7867  11203   2151   -675    -66       N  ",
-                "ATOM      2  CA  ALA A  31     -13.433   7.788  23.746  0.60 65.66           C  ",
-                "ANISOU    2  CA  ALA A  31     7296   7485  10167   1829   -642    -81       C  ",
-                "ATOM      3  C   ALA A  31     -12.448   8.891  24.124  1.00 63.44           C  ",
-                "ANISOU    3  C   ALA A  31     6818   7632   9654   1744   -656     20       C  ",
-                "ATOM      4  O   ALA A  31     -11.549   8.680  24.937  1.00 65.52           O  ",
-                "ANISOU    4  O   ALA A  31     6891   8010   9994   1853   -759    238       O  ",
-                "ATOM      5  CB  ALA A  31     -14.628   7.834  24.691  1.00 66.76           C  ",
-                "ANISOU    5  CB  ALA A  31     7636   7487  10242   1630   -748    143       C  ",
-                "ATOM      6  N   THR A  32     -12.659  10.075  23.550  1.00 58.59           N  ",
-                "ANISOU    6  N   THR A  32     6249   7244   8768   1532   -569   -122       N  ",
-                "ATOM      7  CA  THR A  32     -11.806  11.232  23.788  1.00 55.42           C  ",
-                "ANISOU    7  CA  THR A  32     5678   7217   8161   1398   -585    -53       C  ",
-                "ATOM      8  C   THR A  32     -12.119  11.857  25.148  1.00 52.59           C  ",
-                "ANISOU    8  C   THR A  32     5384   6947   7651   1206   -743    183       C  ",
-                "ATOM      9  O   THR A  32     -13.166  11.585  25.770  1.00 49.84           O  ",
-                "ANISOU    9  O   THR A  32     5232   6412   7293   1142   -803    278       O  ",
-                "ATOM     10  CB  THR A  32     -11.965  12.321  22.696  0.45 51.87           C  ",
-                "ANISOU   10  CB  THR A  32     5260   6948   7500   1220   -449   -254       C  ",
-                "ATOM     11  OG1 THR A  32     -13.294  12.844  22.720  1.00 48.11           O  ",
-                "ANISOU   11  OG1 THR A  32     5047   6344   6890   1015   -452   -282       O  ",
-                "ATOM     12  CG2 THR A  32     -11.660  11.781  21.303  1.00 54.23           C  ",
-                "ANISOU   12  CG2 THR A  32     5483   7248   7875   1398   -280   -508       C  ",
-                ]
-            )
-
-    def test_anisou_with_altloc_maxocc_B(self):
-        """
-        Test anisou_altloc.pdb selects altloc -B.
-
-        anisou.pdb has ANISOU lines. Ensure bug reported in #130 is
-        corrected.
-        """
-        infile = os.path.join(data_dir, 'anisou_altloc.pdb')
-        sys.argv = ['', '-B', infile]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 24)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout,
-            [
-                "ATOM      1  N   ALA A  31     -12.806   6.423  23.735  0.60 70.81           N  ",
-                "ANISOU    1  N   ALA A  31     7836   7867  11203   2151   -675    -66       N  ",
-                "ATOM      2  CA  ALA A  31     -13.433   7.788  23.746  0.40 65.66           C  ",
-                "ANISOU    2  CA  ALA A  31     7296   7485  10167   1829   -642    -81       C  ",
-                "ATOM      3  C   ALA A  31     -12.448   8.891  24.124  1.00 63.44           C  ",
-                "ANISOU    3  C   ALA A  31     6818   7632   9654   1744   -656     20       C  ",
-                "ATOM      4  O   ALA A  31     -11.549   8.680  24.937  1.00 65.52           O  ",
-                "ANISOU    4  O   ALA A  31     6891   8010   9994   1853   -759    238       O  ",
-                "ATOM      5  CB  ALA A  31     -14.628   7.834  24.691  1.00 66.76           C  ",
-                "ANISOU    5  CB  ALA A  31     7636   7487  10242   1630   -748    143       C  ",
-                "ATOM      6  N   THR A  32     -12.659  10.075  23.550  1.00 58.59           N  ",
-                "ANISOU    6  N   THR A  32     6249   7244   8768   1532   -569   -122       N  ",
-                "ATOM      7  CA  THR A  32     -11.806  11.232  23.788  1.00 55.42           C  ",
-                "ANISOU    7  CA  THR A  32     5678   7217   8161   1398   -585    -53       C  ",
-                "ATOM      8  C   THR A  32     -12.119  11.857  25.148  1.00 52.59           C  ",
-                "ANISOU    8  C   THR A  32     5384   6947   7651   1206   -743    183       C  ",
-                "ATOM      9  O   THR A  32     -13.166  11.585  25.770  1.00 49.84           O  ",
-                "ANISOU    9  O   THR A  32     5232   6412   7293   1142   -803    278       O  ",
-                "ATOM     10  CB  THR A  32     -11.965  12.321  22.696  0.55 51.87           C  ",
-                "ANISOU   10  CB  THR A  32     5260   6948   7500   1220   -449   -254       C  ",
-                "ATOM     11  OG1 THR A  32     -13.294  12.844  22.720  1.00 48.11           O  ",
-                "ANISOU   11  OG1 THR A  32     5047   6344   6890   1015   -452   -282       O  ",
-                "ATOM     12  CG2 THR A  32     -11.660  11.781  21.303  1.00 54.23           C  ",
-                "ANISOU   12  CG2 THR A  32     5483   7248   7875   1398   -280   -508       C  ",
-                ]
-            )
-
-    def test_anisou_missing(self):
-        """Test raises error if there are missing anisou lines."""
-        infile = os.path.join(data_dir, 'anisou_missing.pdb')
-        sys.argv = ['', infile]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-
-    def test_captures_previous_residue_maxocc_A(self):
-        """
-        Test the previous residue is not lost.
-
-        In version 2.4.5 we realised that the residue previous to the
-        residues having alternative locations was missing. See dummy_altloc3.pdb.
-        """
-        sys.argv = ['', '-A', os.path.join(data_dir, 'dummy_altloc3.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 25)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout,
-            [
-                "ATOM    153  N   VAL A  18      -0.264 -17.574  22.788  0.00  0.00           N  ",
-                "ATOM    154  CA  VAL A  18      -0.201 -17.901  24.209  0.00  0.00           C  ",
-                "ATOM    155  C   VAL A  18      -1.047 -19.150  24.437  0.00  0.00           C  ",
-                "ATOM    156  O   VAL A  18      -2.260 -19.144  24.214  0.00  0.00           O  ",
-                "ATOM    157  CB  VAL A  18      -0.689 -16.724  25.084  0.00  0.00           C  ",
-                "ATOM    161  N   TRP A  19      -0.408 -20.230  24.882  0.00  0.00           N  ",
-                "ATOM    162  CA  TRP A  19      -1.091 -21.490  25.161  0.00  0.00           C  ",
-                "ATOM    163  C   TRP A  19      -1.303 -21.563  26.667  0.00  0.00           C  ",
-                "ATOM    164  O   TRP A  19      -0.357 -21.920  27.375  0.00  0.00           O  ",
-                "ATOM    165  CB  TRP A  19      -0.272 -22.670  24.635  0.00  0.00           C  ",
-                "ATOM    176  N   TYR A  20      -2.522 -21.226  27.083  0.00  0.00           N  ",
-                "ATOM    177  CA  TYR A  20      -2.898 -21.178  28.493  0.00  0.00           C  ",
-                "ATOM    178  C   TYR A  20      -3.718 -22.410  28.851  0.00  0.00           C  ",
-                "ATOM    179  O   TYR A  20      -4.629 -22.780  28.105  0.00  0.00           O  ",
-                "ATOM    180  CB  TYR A  20      -3.681 -19.898  28.795  0.00  0.00           C  ",
-                "ATOM    189  N   VAL A  21      -3.396 -23.034  29.982  0.40  0.00           N  ",
-                "ATOM    190  CA  VAL A  21      -4.121 -24.205  30.467  0.40  0.00           C  ",
-                "ATOM    191  C   VAL A  21      -4.530 -24.072  31.930  0.40  0.00           C  ",
-                "ATOM    192  O   VAL A  21      -3.835 -23.461  32.747  0.40  0.00           O  ",
-                "ATOM    193  CB  VAL A  21      -3.289 -25.497  30.298  0.40  0.00           C  ",
-                "ATOM    189  N   PRO A  22      -3.396 -23.034  29.982  0.40  0.00           N  ",
-                "ATOM    190  CA  PRO A  22      -4.121 -24.205  30.467  0.40  0.00           C  ",
-                "ATOM    191  C   PRO A  22      -4.530 -24.072  31.930  0.40  0.00           C  ",
-                "ATOM    192  O   PRO A  22      -3.835 -23.461  32.747  0.40  0.00           O  ",
-                "ATOM    193  CB  PRO A  22      -3.289 -25.497  30.298  0.40  0.00           C  ",
-                ]
-            )
-
-    def test_captures_previous_residue_maxocc(self):
-        """
-        Test the previous residue is not lost.
-
-        In version 2.4.5 we realised that the residue previous to the
-        residues having alternative locations was missing. See dummy_altloc3.pdb.
-        """
-        sys.argv = ['', os.path.join(data_dir, 'dummy_altloc3.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 25)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout,
-            [
-                "ATOM    153  N   VAL A  18      -0.264 -17.574  22.788  0.00  0.00           N  ",
-                "ATOM    154  CA  VAL A  18      -0.201 -17.901  24.209  0.00  0.00           C  ",
-                "ATOM    155  C   VAL A  18      -1.047 -19.150  24.437  0.00  0.00           C  ",
-                "ATOM    156  O   VAL A  18      -2.260 -19.144  24.214  0.00  0.00           O  ",
-                "ATOM    157  CB  VAL A  18      -0.689 -16.724  25.084  0.00  0.00           C  ",
-                "ATOM    161  N   TRP A  19      -0.408 -20.230  24.882  0.00  0.00           N  ",
-                "ATOM    162  CA  TRP A  19      -1.091 -21.490  25.161  0.00  0.00           C  ",
-                "ATOM    163  C   TRP A  19      -1.303 -21.563  26.667  0.00  0.00           C  ",
-                "ATOM    164  O   TRP A  19      -0.357 -21.920  27.375  0.00  0.00           O  ",
-                "ATOM    165  CB  TRP A  19      -0.272 -22.670  24.635  0.00  0.00           C  ",
-                "ATOM    176  N   TYR A  20      -2.522 -21.226  27.083  0.00  0.00           N  ",
-                "ATOM    177  CA  TYR A  20      -2.898 -21.178  28.493  0.00  0.00           C  ",
-                "ATOM    178  C   TYR A  20      -3.718 -22.410  28.851  0.00  0.00           C  ",
-                "ATOM    179  O   TYR A  20      -4.629 -22.780  28.105  0.00  0.00           O  ",
-                "ATOM    180  CB  TYR A  20      -3.681 -19.898  28.795  0.00  0.00           C  ",
-                "ATOM    197  N   ALA A  21      -5.676 -24.647  32.284  0.60  0.00           N  ",
-                "ATOM    198  CA  ALA A  21      -6.012 -24.814  33.696  0.60  0.00           C  ",
-                "ATOM    199  C   ALA A  21      -4.990 -25.677  34.426  0.60  0.00           C  ",
-                "ATOM    200  O   ALA A  21      -4.494 -26.675  33.897  0.60  0.00           O  ",
-                "ATOM    201  CB  ALA A  21      -7.405 -25.428  33.847  0.60  0.00           C  ",
-                "ATOM    197  N   GLY A  22      -5.676 -24.647  32.284  0.60  0.00           N  ",
-                "ATOM    198  CA  GLY A  22      -6.012 -24.814  33.696  0.60  0.00           C  ",
-                "ATOM    199  C   GLY A  22      -4.990 -25.677  34.426  0.60  0.00           C  ",
-                "ATOM    200  O   GLY A  22      -4.494 -26.675  33.897  0.60  0.00           O  ",
-                "ATOM    201  CB  GLY A  22      -7.405 -25.428  33.847  0.60  0.00           C  ",
-                ]
-            )
-
-    def test_captures_previous_residue_maxocc_B(self):
-        """
-        Test the previous residue is not lost.
-
-        In version 2.4.5 we realised that the residue previous to the
-        residues having alternative locations was missing. See dummy_altloc3.pdb.
-        """
-        sys.argv = ['', '-B', os.path.join(data_dir, 'dummy_altloc3.pdb')]
-        self.exec_module()
-        self.assertEqual(self.retcode, 0)
-        self.assertEqual(len(self.stdout), 25)
-        self.assertEqual(len(self.stderr), 0)
-        self.assertEqual(
-            self.stdout,
-            [
-                "ATOM    153  N   VAL A  18      -0.264 -17.574  22.788  0.00  0.00           N  ",
-                "ATOM    154  CA  VAL A  18      -0.201 -17.901  24.209  0.00  0.00           C  ",
-                "ATOM    155  C   VAL A  18      -1.047 -19.150  24.437  0.00  0.00           C  ",
-                "ATOM    156  O   VAL A  18      -2.260 -19.144  24.214  0.00  0.00           O  ",
-                "ATOM    157  CB  VAL A  18      -0.689 -16.724  25.084  0.00  0.00           C  ",
-                "ATOM    161  N   TRP A  19      -0.408 -20.230  24.882  0.00  0.00           N  ",
-                "ATOM    162  CA  TRP A  19      -1.091 -21.490  25.161  0.00  0.00           C  ",
-                "ATOM    163  C   TRP A  19      -1.303 -21.563  26.667  0.00  0.00           C  ",
-                "ATOM    164  O   TRP A  19      -0.357 -21.920  27.375  0.00  0.00           O  ",
-                "ATOM    165  CB  TRP A  19      -0.272 -22.670  24.635  0.00  0.00           C  ",
-                "ATOM    176  N   TYR A  20      -2.522 -21.226  27.083  0.00  0.00           N  ",
-                "ATOM    177  CA  TYR A  20      -2.898 -21.178  28.493  0.00  0.00           C  ",
-                "ATOM    178  C   TYR A  20      -3.718 -22.410  28.851  0.00  0.00           C  ",
-                "ATOM    179  O   TYR A  20      -4.629 -22.780  28.105  0.00  0.00           O  ",
-                "ATOM    180  CB  TYR A  20      -3.681 -19.898  28.795  0.00  0.00           C  ",
-                "ATOM    197  N   ALA A  21      -5.676 -24.647  32.284  0.60  0.00           N  ",
-                "ATOM    198  CA  ALA A  21      -6.012 -24.814  33.696  0.60  0.00           C  ",
-                "ATOM    199  C   ALA A  21      -4.990 -25.677  34.426  0.60  0.00           C  ",
-                "ATOM    200  O   ALA A  21      -4.494 -26.675  33.897  0.60  0.00           O  ",
-                "ATOM    201  CB  ALA A  21      -7.405 -25.428  33.847  0.60  0.00           C  ",
-                "ATOM    197  N   GLY A  22      -5.676 -24.647  32.284  0.60  0.00           N  ",
-                "ATOM    198  CA  GLY A  22      -6.012 -24.814  33.696  0.60  0.00           C  ",
-                "ATOM    199  C   GLY A  22      -4.990 -25.677  34.426  0.60  0.00           C  ",
-                "ATOM    200  O   GLY A  22      -4.494 -26.675  33.897  0.60  0.00           O  ",
-                "ATOM    201  CB  GLY A  22      -7.405 -25.428  33.847  0.60  0.00           C  ",
-                ]
-            )
-
-
-    def test_file_not_found(self):
-        """$ pdb_selaltloc not_existing.pdb"""
-
-        afile = os.path.join(data_dir, 'not_existing.pdb')
-        sys.argv = ['', afile]
-
-        self.exec_module()
-
-        self.assertEqual(self.retcode, 1)  # exit code is 1 (error)
-        self.assertEqual(len(self.stdout), 0)  # nothing written to stdout
-        self.assertEqual(self.stderr[0][:22],
-                         "ERROR!! File not found")  # proper error message
-
-    @unittest.skipIf(os.getenv('SKIP_TTY_TESTS'), 'skip on GHA - no TTY')
-    def test_file_missing(self):
-        """$ pdb_selaltloc -A"""
-
-        sys.argv = ['', '-A']
-
-        self.exec_module()
-
-        self.assertEqual(self.retcode, 1)
-        self.assertEqual(len(self.stdout), 0)  # no output
-        self.assertEqual(self.stderr[0],
-                         "ERROR!! No data to process!")
-
-    @unittest.skipIf(os.getenv('SKIP_TTY_TESTS'), 'skip on GHA - no TTY')
-    def test_helptext(self):
-        """$ pdb_selaltloc"""
-
-        sys.argv = ['']
-
-        self.exec_module()
-
-        self.assertEqual(self.retcode, 1)  # ensure the program exited gracefully.
-        self.assertEqual(len(self.stdout), 0)  # no output
-        self.assertEqual(self.stderr, self.module.__doc__.split("\n")[:-1])
-
-    def test_invalid_option(self):
-        """$ pdb_selaltloc -AH data/dummy.pdb"""
-
-        sys.argv = ['', '-AH', os.path.join(data_dir, 'dummy.pdb')]
-
-        self.exec_module()
-
-        self.assertEqual(self.retcode, 1)
-        self.assertEqual(len(self.stdout), 0)
-        self.assertEqual(self.stderr[0][:47],
-                         "ERROR!! Alternate location identifiers must be ")
-
-    def test_not_an_option(self):
-        """$ pdb_selaltloc 20 data/dummy.pdb"""
-
-        sys.argv = ['', '20', os.path.join(data_dir, 'dummy.pdb')]
-
-        self.exec_module()
-
-        self.assertEqual(self.retcode, 1)
-        self.assertEqual(len(self.stdout), 0)
-        self.assertEqual(self.stderr[0],
-                         "ERROR! First argument is not an option: '20'")
-
-
-if __name__ == '__main__':
-    from config import test_dir
-
-    mpath = os.path.abspath(os.path.join(test_dir, '..'))
-    sys.path.insert(0, mpath)  # so we load dev files before  any installation
-
-    unittest.main()