New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4474083 =	38 (1.180)	4 (0.140)	4/226	3.4	11.5%	intergenic (‑30/‑41)	bdcA/bdcR	c‑di‑GMP‑binding biofilm dispersal mediator protein/transcriptional repressor for divergent bdcA
	?	NC_000913	4474127 =	27 (0.920)					coding (4/594 nt)	bdcR	transcriptional repressor for divergent bdcA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAC  >  NC_000913/4473997‑4474270                                                                                                                                   |                                                                                                                                                aCCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtctc                                                                                                                               <  3:36004/149‑1 (MQ=255)         cACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtctc                                                                                                                               <  8:885/141‑1 (MQ=255)         cACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtctc                                                                                                                               >  7:885/1‑141 (MQ=255)               tGCCACCGTGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTcc                                                                                                                        >  5:118667/1‑142 (MQ=255)               tGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTcc                                                                                                                        <  6:118667/142‑1 (MQ=255)                  cACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACACCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTAGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTcc                                                                                                               >  3:125672/1‑148 (MQ=255)                            gagaACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCAAGACGTTTCg                                                                                                    >  7:70863/1‑149 (MQ=255)                             agaACTGTCATACCTGTAAATGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTTCATAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGc                                                                                                   <  4:125672/149‑1 (MQ=255)                              gaACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCt                                                                                                  >  8:27217/1‑149 (MQ=255)                                      ttACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                            >  4:142780/1‑147 (MQ=255)                                                      ccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGgtcgtc                                                                                                                <  6:29965/111‑1 (MQ=255)                                                      ccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGgtcgtc                                                                                                                >  5:29965/1‑111 (MQ=255)                                                        aTAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                         >  2:128690/1‑148 (MQ=255)                                                         tAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaaaa                                                                       <  7:27217/149‑1 (MQ=255)                                                              aCTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCt                                                                  >  1:182868/1‑149 (MQ=255)                                                               cTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                                 <  2:182868/149‑1 (MQ=255)                                                                     gATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTc                                                             <  3:142780/147‑1 (MQ=255)                                                                       tttATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTc                                                                                                     >  8:123014/1‑105 (MQ=255)                                                                       tttATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTc                                                                                                     <  7:123014/105‑1 (MQ=255)                                                                              tGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGt                                                  >  6:162325/1‑149 (MQ=255)                                                                                        aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTg                                        >  6:139471/1‑149 (MQ=255)                                                                                              aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAgtg                                  <  1:128690/149‑1 (MQ=255)                                                                                                    ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGt                                                  <  6:132553/127‑1 (MQ=255)                                                                                                    ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGt                                                  >  5:132553/1‑127 (MQ=255)                                                                                                              cTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTg                   <  4:107929/148‑1 (MQ=255)                                                                                                                gTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa                  >  2:55303/1‑147 (MQ=255)                                                                                                                gTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa                  <  1:55303/147‑1 (MQ=255)                                                                                                                         ttttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTa       <  5:139471/149‑1 (MQ=255)                                                                                                                              tAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAc  >  6:153149/1‑149 (MQ=255)                                                                                                                              tAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAc  <  5:162325/149‑1 (MQ=255)                                                                                                                               aaTGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa                  >  2:73488/1‑132 (MQ=255)                                                                                                                               aaTGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa                  <  1:73488/132‑1 (MQ=255)                                                                                                                                   |                                                                                                                                                ACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAC  >  NC_000913/4473997‑4474270

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