New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	1563562 =	29 (1.670)	4 (0.240)	3/258	2.1	12.8%	coding (2172/2400 nt)	dosP	oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
	?	NC_000913_3_pae_pgi	1563577 =	27 (1.640)					coding (2157/2400 nt)	dosP	oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCG  >  NC_000913_3_pae_pgi/1563414‑1563700                                                                                                                                                     |                                                                                                                                           gAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGa                                                                                                                                           >  2:3655/1‑150 (MQ=255)  aaTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                                                                                                                                          >  2:173947/1‑150 (MQ=255)      aCGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACaaaa                                                                                                                                      >  1:30871/1‑150 (MQ=255)        gcgACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACt                                                                                                                                    >  2:127285/1‑150 (MQ=255)           aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt                                                                                                                                 <  3:52452/150‑1 (MQ=255)                ggATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa                                                                                                                            <  1:173947/150‑1 (MQ=255)                 gATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAt                                                                                                                           >  3:3542/1‑150 (MQ=255)                       gCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAt                                                                                                                     >  2:82552/1‑150 (MQ=255)                              aTCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGtt                                                                                                              >  3:35047/1‑150 (MQ=255)                                  cAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCTTTTTGACATCGTTCGTCAAACCTGTTTCCATTTTTGGTTGCCGtttgtt                                                                                                          >  4:42743/1‑146 (MQ=255)                                   aaaTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACGACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCgg                                                                                                         >  1:156186/1‑150 (MQ=255)                                    aaTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTAAATTTTGTTTTCCGTTTACGGa                                                                                                        <  4:27245/150‑1 (MQ=255)                                                gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTATTTCGGTCAGACAACGATCGACAAACCTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAg                                                                                            >  2:32320/1‑150 (MQ=255)                                                     gACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGa                                                                                       <  1:158973/150‑1 (MQ=255)                                                           ttCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATcc                                                                                 <  2:189035/150‑1 (MQ=255)                                                            tgcGCCACGATGGTTAAATTGAGTATTTGCCCAATGCTGGTAGTGGCATCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCCAAAAAACTTTTGTCAATTTGGATTTCCGTTACCGGAAGACTGACGAAGCGGGATAATCCg                                                                                <  4:112052/148‑1 (MQ=255)                                                            tCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCg                                                                                <  1:132441/150‑1 (MQ=255)                                                                 cACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGaaaa                                                                           <  4:3542/150‑1 (MQ=255)                                                                         ttAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTAc                                                                   >  3:25219/1‑150 (MQ=255)                                                                                   gCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAatcatc                                                         <  4:78209/150‑1 (MQ=255)                                                                                         cccAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGa                                                   <  3:112968/150‑1 (MQ=255)                                                                                                         gCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCata                                   >  2:127768/1‑150 (MQ=255)                                                                                                          cTTCACGTAAGGCAAGCATGCGTTTTTCGGTGAGACCCCGATCGACAAAACTTTTGTCATTTGTGATTTCCGTGACCGGAAGACGGATTATGCGAGATGCTCCGGAAAAGCCCGTACCAAAATCCTCTGCCGATAAGCCAACGCCCatat                                  <  4:35047/150‑1 (MQ=255)                                                                                                          cTTCAAGTAAGGCAAGGATGCGTTTTTCGGTTAGATAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGTAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCatct                                  >  2:16548/1‑148 (MQ=255)                                                                                                               aGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGc                             <  1:3655/150‑1 (MQ=255)                                                                                                                         ggATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGaa                   >  2:152257/1‑150 (MQ=255)                                                                                                                               gTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCt             <  1:82552/150‑1 (MQ=255)                                                                                                                                          accaaaCGACCGACAAACCTTTTGTAAAGTTTGATTTCCGTTACCGGAAGAATGACTAAGAGGGATAATCGGGAAAAGCGCGTACCAAAATCATCAACCGATAAGCCCACGCCCATCTCACGCAGGATCCGGTTGCGCTAAAAGATTTCg  <  2:156186/146‑1 (MQ=255)                                                                                                                                          aGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCg  <  1:127768/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           GAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCG  >  NC_000913_3_pae_pgi/1563414‑1563700

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