New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5835978	25 (0.970)	3 (0.120)	3/276	3.4	11.4%	coding (73/1302 nt)	PP_5115	zinc protease
	?	NC_002947	= 5836015	23 (0.940)					coding (110/1302 nt)	PP_5115	zinc protease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGCACCCCACCTGCGCCATGGTTGAAGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAG  >  NC_002947/5835845‑5836111                                                                                                                                      |                                                                                                                                      cGCACCCCACCTGGGCCTTGGTTGAAGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAAGCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGaa                                                                                                                         <  4:124930/148‑1 (MQ=255)   cACCCCACCTGCGCCATGGTTGAAGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGa                                                                                                                       <  2:125398/148‑1 (MQ=255)                 cATGGTTGAAGCTGCACCTTCATTTCCTGAAGGAGCAGGTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCAc                                                                                                         <  4:72009/148‑1 (MQ=255)                          aGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGc                                                                                                 <  7:63133/147‑1 (MQ=255)                             tGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACAACCGTGCGCCACTGGTCAGCGTGc                                                                                             <  4:98873/148‑1 (MQ=255)                              gCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCa                                                                                            >  6:16086/1‑148 (MQ=255)                                 ccTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGaa                                                                                                                         >  3:31767/1‑116 (MQ=255)                                  cTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAg                                                                                                                        <  4:31767/116‑1 (MQ=255)                                   ttCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGGTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCtgt                                                                                       <  3:102870/148‑1 (MQ=255)                                          ccTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTGCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGGCAGCGTGCAGCTGTGGTAcc                                                                                 <  8:33733/147‑1 (MQ=255)                                              aaGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTa                                                                            >  2:23532/1‑148 (MQ=255)                                                     aGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCa                                                                     <  5:113044/148‑1 (MQ=255)                                                                          cccTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCcaca                                                >  2:56137/1‑148 (MQ=255)                                                                           ccTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCacac                                               >  6:58644/1‑148 (MQ=255)                                                                                   ccgcTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGt                                       <  5:16086/148‑1 (MQ=255)                                                                                                ccTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGACTGTCCCATGCCTTGGa                          >  7:112723/1‑148 (MQ=255)                                                                                                 cTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGa                          >  1:37354/1‑147 (MQ=255)                                                                                                  tCGAAAGGTTCTTTTACACTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGAACACCGTGCTCCACTGGTAAGCGTGCAGCTGTGGTACCACGTAGGCTCAAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGc                        >  1:99994/1‑148 (MQ=255)                                                                                                     aaaGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCAt                      >  3:50274/1‑147 (MQ=255)                                                                                                                   ccTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTCCCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAgg         >  3:69878/1‑146 (MQ=255)                                                                                                                    cTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGgca      >  2:21047/1‑148 (MQ=255)                                                                                                                     tCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGa                                                              >  4:120604/1‑91 (MQ=255)                                                                                                                     tCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGgcc      >  8:32787/1‑146 (MQ=255)                                                                                                                      cGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGaa                                                             <  3:120604/91‑1 (MQ=255)                                                                                                                         cAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCGGTCCCATGCCGTGGAGCATCTACTGTTCTAGGgacgcag  >  5:78775/1‑147 (MQ=255)                                                                                                                                      |                                                                                                                                      CGCACCCCACCTGCGCCATGGTTGAAGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAG  >  NC_002947/5835845‑5836111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.