New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2026671 =	39 (0.720)	3 (0.060)	3/250	NT	7.2%	coding (1347/1695 nt)	dgcQ	putative diguanylate cyclase DgcQ
	?	NC_000913	2026688 =	40 (0.780)					coding (1330/1695 nt)	dgcQ	putative diguanylate cyclase DgcQ

CGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGT  >  NC_000913/2026529‑2026808                                                                                                                                               |                                                                                                                                          cGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCtga                                                                                                                                           <  2:724229/143‑1 (MQ=255) cGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCtga                                                                                                                                           >  1:51418/1‑143 (MQ=255)  agccgtggcacgacaccttaacgcgtttatatcgccgtggcgcactgttcgaaaaagcccgtccgctcgctaaattgtgtcagacgcaccaacatcctttttctgtcattcaggtcgatcttgaccattttaaagcgattaAt                                                                                                                                          >  2:974725‑M1/142‑143 (MQ=255)     cTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCCCCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTCCTAATTAATACGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGAc                                                                                                                                       >  1:594227/1‑143 (MQ=255)                   ttaacgcgtttatataaccgtggcgcactgtgcgaaaaagcccgtccgctcgctaaattgtgtcagacgcaccaacatcctttttctgtcattcaggtcgatcttgaccattttaaagcgattaATGACCAAAGCGGTCATTa                                                                                                                         >  1:94269‑M1/125‑143 (MQ=255)                          cgcgCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCtt                                                                                                                  >  1:103444/1‑143 (MQ=255)                                 ggCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTCAAATg                                                                                                           >  1:967187/1‑143 (MQ=255)                                 ggCAGAATCACACAAAACTCCGCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATg                                                                                                           <  1:259988/143‑1 (MQ=255)                                        tCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGa                                                                                                    <  1:594417/143‑1 (MQ=255)                                        tCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCACCATGGGTCATCACACCAGCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGa                                                                                                    >  1:486762/1‑143 (MQ=255)                                        tCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACCGCCAACCAATCCCGCAGCACGAGAAAGAACACGGTCGCCCGCCAGATGACCAAACCGGGCAGTAACCGCTGTCAAATGATCAAgg                                                                                                    >  2:943562/1‑142 (MQ=255)                                           cccaacgctcctcACCACCGACCCGCCCGGCAACGTCCTGCGCCCGCAAGGAACTGCTAATTAACCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCg                                                                                                 <  1:983078/136‑1 (MQ=255)                                            acaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGa                                                                                                >  1:769852/1‑143 (MQ=255)                                            acaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCACCAGCAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGCTCATTAATCGCTCTAAAATGGTCCACATcaa                                                                                                >  2:927744/1‑141 (MQ=255)                                                  ctcctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCAAATAACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGaa                                                                                          <  1:152773/143‑1 (MQ=255)                                                   tcctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGCTCGCCCTGGAt                                                                                         >  1:928833/1‑143 (MQ=255)                                                   tcctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACGGCGACTGACTCCGGCCCAGCGACATTATGCACGGGCGCACGCGTTCTGACCACAGCGGCCATGAATCGCTTTCCAAAGGTCACGAGCGACCtccat                                                                                         >  2:854093/1‑139 (MQ=255)                                                      tcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGAc                                                                                      <  2:483653/143‑1 (MQ=255)                                                       caccacCGACCCTCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACa                                                                                     <  2:657629/143‑1 (MQ=255)                                                        accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAAAGACAg                                                                                    >  2:869205/1‑143 (MQ=255)                                                          caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGGACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAGGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTCAATGACAGaa                                                                                  >  1:683647/1‑143 (MQ=255)                                                             cgcgccgcccgGCAACGTCCCGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaaa                                                                               <  1:913065/139‑1 (MQ=255)                                                                 ccgcccgGCACCGTCCTGCCCACGCAAGAAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGGCTGATGACCAAAGCGTTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                           <  2:128972/143‑1 (MQ=255)                                                                 ccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGCGAAAGAATACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                           <  2:121775/143‑1 (MQ=255)                                                                 ccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                           <  2:508145/143‑1 (MQ=255)                                                                    cccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGtt                                                                        >  2:274339/1‑143 (MQ=255)                                                                       ggCACCGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGt                                                                     <  2:193674/143‑1 (MQ=255)                                                                        gCAACGTACGGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTg                                                                    <  2:490760/143‑1 (MQ=255)                                                                         cAACGTGCTGCGACCGCGAGTAACTGCTGATTAACCCGGCAGCACGAGAAAGAACACGGTCGCCCGCCTGAGGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGc                                                                   <  1:992894/143‑1 (MQ=255)                                                                                    gacgcaccaacatcctttttctgtcattcaggtcgatcttgaccattttaaagcgattaATGACCAAAGCGGTCATTAATCGc                                                                                                                    >  2:322354‑M1/60‑83 (MQ=255)                                                                                    gacgcaccaacatcctttttctgtcattcaggtcgatcttgaccattttaaagcgattaATGACCAAAGCGGTCATTAATCGc                                                                                                                    <  1:322354‑M1/24‑1 (MQ=255)                                                                                          cggggaCTGCTAATTAATCCGGCAGCATGAGATAGAACACGGTCGCGCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagc                                                  <  2:904724/138‑1 (MQ=255)                                                                                          aaGGAACTGCTAAGTGATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagc                                                  <  2:826751/143‑1 (MQ=255)                                                                                                   cTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCGAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACGGAAAAAGGATGTTGGTGCGTGTGCCAAATTTTAGCTAGCAGagg                                         >  1:727654/1‑141 (MQ=255)                                                                                                   cTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACg                                         >  1:1030114/1‑143 (MQ=255)                                                                                                      cttgagccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGc                                      <  2:993449/137‑1 (MQ=255)                                                                                                           gccGGCAGCATGAGAAAGAACACGGTCGCCCGGCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGGCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttttt                                 <  2:855561/142‑1 (MQ=255)                                                                                                                       agaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTTCGTCTGACACAATTTAGCGAGCTGACGGTGTTTTTCGAACAGTgcgc                     >  1:55541/1‑143 (MQ=255)                                                                                                                        gaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGTCCTGCTTGTCAGAAATAGGATTTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTCTCGAACAGTGCGcc                    >  1:19586/1‑143 (MQ=255)                                                                                                                            gaaCACGGTCGCCCGCCTGATGACCGAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACgg                <  2:535581/143‑1 (MQ=255)                                                                                                                               cacGGTGGCCCGCCTGATGACCAAGGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTta             <  2:967187/143‑1 (MQ=255)                                                                                                                               cacGGTCGCCCGCCTGATGACCAACGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGTAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTta             <  1:61105/143‑1 (MQ=255)                                                                                                                                 cGGGCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTAAGCGAGCGGAAGGGCTTTTTCGAACAGTGCGCCACGGTtata           <  1:213080/143‑1 (MQ=255)                                                                                                                                   gTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGAcatata         >  2:132649/1‑143 (MQ=255)                                                                                                                                      gcccgccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAc      <  2:486762/143‑1 (MQ=255)                                                                                                                                      gcccgccTGATAGCCGCAGCGTCTATTAGTCGCTTTAAAATGGTGAAGATCGCTGTGAAGGGCGGAAAAAGGATGTTGGTGCCTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAc      <  2:401964/143‑1 (MQ=255)                                                                                                                                          gccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGCTCTATAAACGCGt  >  2:805316/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                          CGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGT  >  NC_000913/2026529‑2026808

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