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.

GGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCGGTATCGTGTTCCA  >  NC_000913_3_pae_pgi/1563429‑1563713                                                                                                                                                     |                                                                                                                                         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)                                                                                                                                        ataaaCTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCGGTATCGTGTTCca  <  1:15926/148‑1 (MQ=255)                                                                                                                                           aaCTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCGGTATCgtgt      >  3:61127/1‑143 (MQ=255)                                                                                                                                           aaCTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCGGTATCgtgt      <  4:61127/143‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         GGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATTTCCGTTACCGGAAGACTGACTAAGCGGGATAATCCGGAAAAGCCCGTACCAAAATCATCTACCGATAAGCCCACGCCCATATCACGCAGGATCTGAATGCGCTTAAAGATTTCGGTATCGTGTTCCA  >  NC_000913_3_pae_pgi/1563429‑1563713

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