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.

AGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCATGACCCGCTTGGGCG  >  NC_002947/5835870‑5836158                                                                                                                                                  |                                                                                                                                                aGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGc                                                                                                                                                <  7:63133/147‑1 (MQ=255)    tGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACAACCGTGCGCCACTGGTCAGCGTGc                                                                                                                                            <  4:98873/148‑1 (MQ=255)     gCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCa                                                                                                                                           >  6:16086/1‑148 (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)                                                                                                                              ccaccGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACGACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGTCGCCTGGCCAGTACTCAACCCtt                  <  8:58256/148‑2 (MQ=255)                                                                                                                               caccGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCa                 >  5:86666/1‑148 (MQ=255)                                                                                                                                     gcgcCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCATGAccc           >  4:16075/1‑148 (MQ=255)                                                                                                                                       gcCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCATGACCCGc         >  6:53682/1‑148 (MQ=255)                                                                                                                                         cACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCATGACCCGCtt       >  5:109775/1‑148 (MQ=255)                                                                                                                                              gTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCATGACCCGCTTGggcg  >  3:6614/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                AGCTGCACCTTCATTTCCTGAAGGAGCAGTTCAACCATGGCCGAACACCCCTCCCCGCCGCTAATGGCAGCCTCGAAAGGTTCTTTTACCCTCGCCAACGGCCTGCACGTCTGCCTGCGTGAAGACCACCGTGCGCCACTGGTCAGCGTGCAGCTGTGGTACCACGTAGGCTCCAGCTACGAACCTGAAGGCCACACCGGCCTGTCCCATGCCTTGGAGCATCTACTGTTCGAGGGCAGCAGCAAACTGGCGCCTGGCCAGTACTCAACCCTCATGACCCGCTTGGGCG  >  NC_002947/5835870‑5836158

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.