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.

AACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCG  >  NC_002947/3820738‑3821030                                                                                                                                                   |                                                                                                                                                   aaCCGCAACAGCCACGTCTATCCCCTGTCAGGCGGCTGCCCGGGCTACAACCCGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCaaa                                                                                                                                                   <  2:128920/148‑1 (MQ=255)    cGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGc                                                                                                                                                >  7:34798/1‑148 (MQ=255)       aaCAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTgc                                                                                                                                             >  1:17162/1‑148 (MQ=255)               cGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGcccc                                                                                                                                     >  3:53131/1‑148 (MQ=255)                  cTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCtgt                                                                                                                                  <  3:153870/148‑1 (MQ=255)                      cACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTgg                                                                                                                                <  7:90438/146‑1 (MQ=255)                      cACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGGTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTgg                                                                                                                                <  3:80407/146‑1 (MQ=255)                                     tGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAg                                                                                                               >  8:140958/1‑148 (MQ=255)                                                            aGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCtgtg                                                                                                                                 >  1:147079/1‑107 (MQ=255)                                                             gCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTgg                                                                                                                                <  2:147079/107‑1 (MQ=255)                                                                       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)                                                                                                                                                   |                                                                                                                                                   AACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGTTACAGCGAAGAACCAGGTTACCTTCGCAAACGAAGGTGAAGCGCAGGCAGCGGGGTATCGCAAGGCAGGCAACTGCCGCTGACCGTTTCAAAGGCTGCGCGGCCCCTGTGGGAGCGGGGGCGCGTCAGCCGTTGACCAGGTTTACCGGCCGCTCACCGGCCAGCGCGCTGAGCAGGTTGTCCACCGCGCAACGCGCCATGGCCTGGCGGGTTTCCTCGGTCGCCGAGCCGATATGCG  >  NC_002947/3820738‑3821030

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.