New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4438181	31 (0.600)	4 (0.080)	4/248	NT	11.0%	coding (1341/1425 nt)	PP_3932	diguanylate cyclase
	?	NC_002947	= 4438272	36 (0.740)					coding (757/762 nt)	PP_3933	membrane protein

ACTTCACCGCCGAACGCCTGCGCCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAG  >  NC_002947/4438039‑4438329                                                                                                                                               |                                                                                                                                                     aCTTCACCGCCGAACGCCTGCGCCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGcc                                                                                                                                                >  5:14075/1‑149 (MQ=255)        cgccgAACGCCTGCGCCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCggctagtt                                                                                                                                              >  6:287231‑M1/1‑136 (MQ=255)        cgccgAACGCCTGCGCCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAgg                                                                                                                                         >  7:251708/1‑149 (MQ=255)               cGCCTGCGCCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTa                                                                                                                                  >  5:36781/1‑149 (MQ=255)                    gcgccgcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCggg                                                                                                                             <  1:25567/149‑1 (MQ=255)                         gcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAg                                                                                                                        <  5:109959/149‑1 (MQ=255)                         gcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAg                                                                                                                        >  2:156817/1‑149 (MQ=255)                          cgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAgg                                                                                                                       <  3:277143/149‑1 (MQ=255)                           gccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGc                                                                                                                      >  4:145063/1‑149 (MQ=255)                                  aGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCggctagttcagatcaccaggcgctgtgcttcggcctgtt                                                                                                               >  2:89793‑M1/1‑110 (MQ=255)                                          cACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGgc                                                                                                                                        <  5:160532/116‑1 (MQ=255)                                          cACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGgc                                                                                                                                        >  6:160532/1‑116 (MQ=255)                                          cACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAAc                                                                                                         <  5:167468/147‑1 (MQ=255)                                           aCGGGACACCGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACcgc                                                                                                      <  6:261029/149‑1 (MQ=255)                                                          aCCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGaa                                                                                       <  7:99846/149‑1 (MQ=255)                                                              gcgcCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCACCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGca                                                                                    >  2:291067/1‑148 (MQ=255)                                                                     cAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCt                                                                            <  2:234342/149‑1 (MQ=255)                                                                       gATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTgg                                                                          >  4:187106/1‑149 (MQ=255)                                                                            ccGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATc                                                                     >  7:181249/1‑149 (MQ=255)                                                                                gAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGaa                                                                 <  8:300026/149‑1 (MQ=255)                                                                                gAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGaa                                                                 <  2:113603/149‑1 (MQ=255)                                                                                gAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGCGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGaa                                                                 >  5:48926/1‑149 (MQ=255)                                                                                 aGTATCGGTTTGTCAGGCTGCGACGCAGACATTCCGGCCACCATCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGaa                                                                 <  6:36781/148‑1 (MQ=255)                                                                                 aGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGaa                                                                 <  4:246577/148‑1 (MQ=255)                                                                                      cGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCGCGACGAACTGCTGGCCggctagttcagatcaccaggcgctgtgcttcggcctgttcgacgcggttgcggccccgggccttggcccggtacagcgcctggtcggcac                                                            <  5:287231‑M1/148‑91 (MQ=255)                                                                                       ggTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCg                                                          >  1:156817‑M2/1‑148 (MQ=255)                                                                                        gtgtGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCgg                                                         <  7:289626‑M2/149‑3 (MQ=255)                                                                                              cAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCggctagttcagatcaccaggcgctgtgcttcggcctgttcgactcggttgcggccccgggccttggcccggtacagcgcctggtcggcactggccagca                                                   <  1:89793‑M1/149‑100 (MQ=255)                                                                                              cAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCggccagca                                                   >  6:14075‑M2/1‑141 (MQ=255)                                                                                               aGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTCGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCggccagcag                                                  <  1:231386‑M2/149‑10 (MQ=255)                                                                                               aGGCGGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACCGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCggccagcag                                                  >  7:189400‑M2/1‑140 (MQ=255)                                                                                                    gCGACGCATACATGCCGGCCCCCAGCCACGACGAACTGCTGGCCAGTGCCGACCAGCCGCTGTAACGGGCCAAGGCCCGGGGCCGCAACCCCGTCGAACAGGCCGAAGCACCGCGCCTGGTCATCTGAACTAGCCGGACAGCAGGTccc                                             >  6:44997/1‑149 (MQ=255)                                                                                                     cGACGCAGACATGCCGGCCCCCACCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCa                                            <  3:145063/149‑1 (MQ=255)                                                                                                         gCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGAc                                        <  8:251708/149‑1 (MQ=255)                                                                                                              cATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCt                                    <  2:303686/148‑1 (MQ=255)                                                                                                               aTGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                  <  3:187106/149‑1 (MQ=255)                                                                                                                tGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTgg                                 <  2:231386/149‑1 (MQ=255)                                                                                                                      cccccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCa                           <  3:227265/149‑1 (MQ=255)                                                                                                                           aGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcaca                      <  1:280853/149‑1 (MQ=255)                                                                                                                                    gaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGc             <  7:48858/149‑1 (MQ=255)                                                                                                                                       ctgctgGCCggctagttcagatcaccaggcgctgtgcttcggcctgttcgacgcggttgcggccccgggccttggcccggtacagcgcctggtcggcactggcgagcagttcgtcgaggctgggggccggcatgtctgcgtcgcagcct          <  1:231386‑M1/149‑141 (MQ=255)                                                                                                                                       ctgctgGCCggctagttcagatcaccaggcgctgtgcttcggccggttcgacgcggttgcggccccgggccttggcccggtacagcgcctggtcggcactggccagcagttcgtcgaggctgggggccggcatgtctgcgtcgccgcct          >  7:189400‑M1/1‑9 (MQ=255)                                                                                                                                        tgctgGCCggctagttcagatcaccaggcgctgtgcttcggcctgttcgacgcggttgcggccccgggccttggcccggtacagcgcctggtcggcactggccagcagttcgtcgaggctgggggccggcatgtctgcgtcgcagcctg         >  6:14075‑M1/1‑8 (MQ=255)                                                                                                                                        tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGAGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCggccagcagttcgtcgaggctgggggccggcatgtctgcgtcgcagcctg         <  1:89793‑M2/149‑51 (MQ=255)                                                                                                                                           ggtacAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAgg       >  1:94989/5‑148 (MQ=255)                                                                                                                                              ccggctagttcagatcaccaggcgctgtgcttcggcctgttcgacgcggttgcggccccgggccttggcccggtacagcgcctggtcggcactggccagcagttcgtcgaggctgggggccggcatgtctgcgtcgcagcctgccacac   <  7:289626‑M1/149‑148 (MQ=255)                                                                                                                                              ccAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACa   <  8:181249/149‑1 (MQ=255)                                                                                                                                               cggctagttcagatcaccaggcgctgtgcttcggcctgttcgacgcggttgcggccccgggccttggcccggtacagcgcctggtcggcactggccagcagttcgtcgaggctgggggccggcatgtctgcgtcgcagcctgccacacc  >  1:156817‑M1/1‑1 (MQ=255)                                                                                                                                               cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAg                          >  4:95152/1‑125 (MQ=255)                                                                                                                                               cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAg  >  5:232403/1‑149 (MQ=255)                                                                                                                                               cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAg  >  2:14088/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     ACTTCACCGCCGAACGCCTGCGCCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAG  >  NC_002947/4438039‑4438329

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