New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3820884	24 (0.830)	5 (0.180)	5/282	3.1	18.4%	intergenic (+9/+34)	endA/ptxD	endonuclease I/phosphonate dehydrogenase
	?	NC_002947	= 3820924	21 (0.750)					coding (957/963 nt)	ptxD	phosphonate dehydrogenase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTC  >  NC_002947/3820808‑3821070                                                                                                                     |                                                                                                                                                   aGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCa                                                                                                                     >  1:107036/1‑148 (MQ=255)         ttCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAgcgcgc                                                                                                               >  3:10349/1‑146 (MQ=255)            gCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGa                                                                                                            >  3:24633/1‑146 (MQ=255)               aaCGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAgg                                                                                                       <  7:137304/148‑1 (MQ=255)                           gcgcAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACcgc                                                                                            >  7:3536/1‑147 (MQ=255)                                        ggTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCCCCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATgg                                                                              >  1:148843/1‑148 (MQ=255)                                            tCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCACCGCGCCATGGCCTg                                                                          >  8:94860/1‑148 (MQ=255)                                               ctaggcaggcAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGc                                                                        <  5:91612/145‑1 (MQ=255)                                                   gcaggcaACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGGCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACCCGCCAGGCCATGGCGCGtt                                                                   <  2:16939/148‑1 (MQ=255)                                                       gcaACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCt                                                               <  8:34798/148‑1 (MQ=255)                                                        caACTGCCGCTGACCGTGTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTc                                                              <  2:9745/148‑1 (MQ=255)                                                                        tttCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGCGCCGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGata                                              >  5:54891/1‑148 (MQ=255)                                                                            aaaGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCg                                          >  5:36058/1‑148 (MQ=255)                                                                            aaaGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCg                                          >  7:115173/1‑148 (MQ=255)                                                                              aGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCggg                                        <  3:45582/148‑1 (MQ=255)                                                                                           ccTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACg                            >  5:49612/1‑147 (MQ=255)                                                                                            cTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGtt                          <  2:107036/148‑1 (MQ=255)                                                                                             tgtgGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTa                         <  2:17162/148‑1 (MQ=255)                                                                                                        gggCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCaa              <  1:80366/148‑1 (MQ=255)                                                                                                              gTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAgcgg        <  8:3536/148‑1 (MQ=255)                                                                                                                    cGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTc  <  6:143041/148‑1 (MQ=255)                                                                                                                     |                                                                                                                                                   AGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCGGGGTCGCTACCACGTTATCCAGCTGCAACAGCGGCGAGTC  >  NC_002947/3820808‑3821070

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.