New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1964564 =	17 (1.050)	3 (0.200)	3/238	NT	16.8%	coding (733/2526 nt)	yegE	phosphodiesterase domain‑containing diguanylate cyclase
	?	NC_002947	1964605 =	14 (0.930)					coding (774/2526 nt)	yegE	phosphodiesterase domain‑containing diguanylate cyclase

CCGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCC  >  NC_002947/1964459‑1964750                                                                                                                                                   |                                                                                                                                                  ccGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTCAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCttt                                                                                                                                                >  1:294728/1‑150 (MQ=255)   ggcggcCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGc                                                                                                                                              <  1:222991/150‑1 (MQ=255)               cTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCa                                                                                                                                  <  2:22649/150‑1 (MQ=255)                     tGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGgaaga                                                                                                                            >  2:360672/1‑150 (MQ=255)                         gCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCTCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGaaaaa                                                                                                                        <  2:405367/150‑1 (MQ=255)                            ttGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAagag                                                                                                                     >  1:188636/1‑150 (MQ=255)                            gttgtaatcaccgcatccccgatggcctgcagggtaatctgcgcacgctctttttcttcctgcagggcgctggcaaagGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCttt                                                                                                                                                >  1:327541‑M1/79‑123 (MQ=255)                            gttgtaatcaccgcatccccgatggcctgcagggtaatctgcgcacgctctttttcttcctgcagggcgctggcaaagGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCttt                                                                                                                                                <  2:327541‑M1/45‑1 (MQ=255)                                                 ttGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAg                                                                                                <  2:188636/150‑1 (MQ=255)                                                               gctgATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGgaa                                                                                                                              <  2:65385/106‑1 (MQ=255)                                                               gctgATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGgaa                                                                                                                              >  1:65385/1‑106 (MQ=255)                                                                    tCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGa                                                                             <  1:360672/150‑1 (MQ=255)                                                                         aaCCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACa                                                                        >  2:216476/1‑150 (MQ=255)                                                                          aCCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACaa                                                                       <  2:294729/150‑1 (MQ=255)                                                                                 gACCGCAATGTGCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGt                                                                <  1:442075/150‑1 (MQ=255)                                                                                               tGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCg                                                  >  1:45982/1‑150 (MQ=255)                                                                                                      aaaggcctggcgctgggccaacagcttgctggagcggcgccaggcCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAAc                                                                      <  2:327541‑M2/78‑1 (MQ=255)                                                                                                      aaaggcctggcgctgggccaacagcttgctggagcggcgccaggcCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAAc                                                                      >  1:327541‑M2/46‑123 (MQ=255)                                                                                                          gcctggcgctgggccaacagcttgctggagcggcgccaggcCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAAc                                       >  2:104364‑M2/42‑150 (MQ=255)                                                                                                          GCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAAc                                       >  1:155778‑M1/1‑150 (MQ=255)                                                                                                                                    ggCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTg             <  2:45982/150‑1 (MQ=255)                                                                                                                                       ccAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGCTGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCa          >  2:217446/1‑150 (MQ=255)                                                                                                                                        cAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAgg                                                           >  2:345812/1‑100 (MQ=255)                                                                                                                                        cAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGAGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAgg                                                           <  1:345812/100‑1 (MQ=255)                                                                                                                                               aggcCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGcc  >  2:131010‑M2/5‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CCGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCC  >  NC_002947/1964459‑1964750

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