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

TCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGC  >  NC_000913/2026547‑2026830                                                                                                                                              |                                                                                                                                               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                                                                                                                                          <  1:322354‑M1/24‑1 (MQ=255)                                                                  gacgcaccaacatcctttttctgtcattcaggtcgatcttgaccattttaaagcgattaATGACCAAAGCGGTCATTAATCGc                                                                                                                                          >  2:322354‑M1/60‑83 (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)                                                                                                                      gcgattaatgaccgctttggtcatTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTc                                                                                      <  1:322354‑M2/59‑1 (MQ=255)                                                                                                                      gcgattaatgaccgctttggtcatTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTc                                                                                      >  2:322354‑M2/25‑83 (MQ=255)                                                                                                                        gccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGCTCTATAAACGCGt                        >  2:805316/1‑143 (MQ=255)                                                                                                                           taatgaccgctttggtcatTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGCACAGTGCGCCACGGTTATATAAACGCGTTaa                     >  1:94269‑M2/20‑143 (MQ=255)                                                                                                                                 CCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGCCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCCCGGTTCTATAAACGCgttaagctgcc               <  1:1010400‑M1/143‑6 (MQ=255)                                                                                                                                  cagaGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCg              <  1:783194/140‑1 (MQ=255)                                                                                                                                   aaaGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTAGATAAACGCGTAAAGGTGTCGt             >  2:509960/1‑143 (MQ=255)                                                                                                                                            atTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGCGATATAAACGCGTTAAGGTGTCGTGCCACGGCt    >  2:974725‑M2/3‑143 (MQ=255)                                                                                                                                              tAATCGTTTTAAAATGGGCATGATCTACCTGAATGAATGAAGAAGGATGTTGGTGCGTCTGACACAATTTGGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGc  <  2:354902/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               TCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGC  >  NC_000913/2026547‑2026830

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.