New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	964991 =	27 (1.670)	3 (0.190)	3/250	NT	10.3%	coding (794/828 nt)	phnE	phosphonate ABC transporter permease
	?	NC_002947	965023 =	26 (1.650)					coding (826/828 nt) coding (2/768 nt)	phnE ptxC	phosphonate ABC transporter permease phosphonate transport system permease PtxC

TTTGTCGGTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGG  >  NC_002947/964876‑965134                                                                                                                    |                                                                                                                                                tttGTCGGTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGa                                                                                                               >  2:433516/1‑150 (MQ=255)   tGTCGGTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAAc                                                                                                             <  2:286725/150‑1 (MQ=255)    gTCGGTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAGcc                                                                                                            >  2:147210/1‑150 (MQ=255)      cagcccgccattgccgaccagtgcgtgcaggtccagctgcagatacgcgaacgaggccaccgccgctgccacgatggcgcccagcaacacccagttgatcacccggttcaTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCgg                                                                                                          >  2:300986‑M1/111‑150 (MQ=255)       ggTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCggg                                                                                                         >  2:392942/1‑150 (MQ=255)        gTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCgggg                                                                                                        >  2:388815/1‑149 (MQ=255)              gcgggcggGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGGGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCaa                                                                                                  >  2:351920/1‑150 (MQ=255)                gggCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGTGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGCGGATGGCATCAACCGGGTGCACAACt                                                                                                >  1:425430/1‑150 (MQ=255)                  gCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTgg                                                                                              <  1:290548/150‑1 (MQ=255)                   cTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTggg                                                                                             <  1:220628/150‑1 (MQ=255)                              gATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGATCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGgc                                                                                  <  1:106793/150‑1 (MQ=255)                                gcaggtccagctgcagatacgcgaacgaggccaccgccgctgccacgatggcgcccagcaacacccagttgatcacccggttcaTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGt                                                                                                        >  2:409831‑M1/85‑126 (MQ=255)                                gcaggtccagctgcagatacgcgaacgaggccaccgccgctgccacgatggcgcccagcaacacccagttgatcacccggttcaTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGt                                                                                                        <  1:409831‑M1/42‑1 (MQ=255)                                        tCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTgg                                                                        >  1:287269/1‑150 (MQ=255)                                                  tGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGt                                                              <  1:325843/150‑1 (MQ=255)                                                  tGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGt                                                              >  2:315871/1‑150 (MQ=255)                                                  tGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTATTGTAGTGGGGTATCCTTGCCTACCTGCTCAGCCCTCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGt                                                              <  1:433516/150‑1 (MQ=255)                                                               gcAAGTGGCGAGCATGCTGTTGACGTTCTTGGTGCTAGGGATCCTGGCCGACCTTCTCAGCCGCCTGCTGCGCTTGATGCTGGCATGAACCGGGTGATCAACTGGGTGTTGGTGGGCGGTATCGTGGCAGTGGCGGTGGCGTCGTTcgcg                                                 <  2:425430/150‑1 (MQ=255)                                                                  aGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAGCCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTAt                                              <  1:147210/150‑1 (MQ=255)                                                                        gAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATCGGCAATGTGATCAACTGGGTGTTGCTGGGCGACATCATAGATCTGGCGTCGGGTCCGTTCGCGTATCTGCAg                                        >  1:342718/1‑150 (MQ=255)                                                                             tGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGa                                   <  1:392942/150‑1 (MQ=255)                                                                                       gTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGgtgt                                                                                          <  1:351389/85‑1 (MQ=255)                                                                                       gTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGgtgt                                                                                          >  2:351389/1‑85 (MQ=255)                                                                                                  tGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTc                   <  2:91418/145‑1 (MQ=255)                                                                                                  tGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTc                   >  1:91418/1‑145 (MQ=255)                                                                                                          cTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGc      >  1:116362/1‑150 (MQ=255)                                                                                                             gCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTgg   <  1:387645/150‑1 (MQ=255)                                                                                                              ccGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTc                                                        <  2:172051/96‑1 (MQ=255)                                                                                                              ccGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTc                                                        >  1:172051/1‑96 (MQ=255)                                                                                                              ccGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTggg  <  1:294335/150‑1 (MQ=255)                                                                                                                    |                                                                                                                                                TTTGTCGGTGCTGGCGGGCTGGGTCAGCAGATCGACCTGTCGATGCGCATGTTTGCGGGTGCGCAAGTGGCGAGCATGCTGTTGACGTTCTTGTTGCTGGTGATCCTGGCCGACCTGCTCAGCCGCCTGCTGCGCTGGAGGCTGGCATGAACCGGGTGATCAACTGGGTGTTGCTGGGCGCCATCGTGGCAGCGGCGGTGGCCTCGTTCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGG  >  NC_002947/964876‑965134

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ 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.