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

GCGATTGGAAGGCCCGCATTGTCGCCATCAACGAAGGCGTGACCCCGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAAC  >  NC_002947/1964415‑1964713                                                                                                                                                      |                                                                                                                                                      gCGATTGGAAGGCCCGCATTGTCGCCATCAACGAAGGCGTGACCCCGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTg                                                                                                                                                       >  2:155778/1‑150 (MQ=255)    ggccggattgccagtgggtcaactgctcggcggccgggttcatgtagccgatattgccctcgacatcggttgtaatcaccgcatccccgatggcctgcagggtaatctgcgcacgctctttttcttcctgcagggcgctggcaaagGCCt                                                                                                                                                    >  2:131010‑M1/147‑150 (MQ=255)                                     gCGTGACCCCGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGcc                                                                                                                   >  1:22649/1‑150 (MQ=255)                                         gttcatgtagccgatattgccctcgacatcggttgtaatcaccgcatccccgatggcctgcagggtaatctgcgcacgctctttttcttcctgcagggcgctggcaaagGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGc                                                                                                               >  2:104364‑M1/110‑150 (MQ=255)                                             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)                                                                       cTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTg                                                                                                                                <  2:386342/103‑1 (MQ=255)                                                                       cTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTg                                                                                                                                >  1:386342/1‑103 (MQ=255)                                                                        ttGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAagag                                                                                >  1:188636/1‑150 (MQ=255)                                                                        gttgtaatcaccgcatccccgatggcctgcagggtaatctgcgcacgctctttttcttcctgcagggcgctggcaaagGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCttt                                                                                                           <  2:327541‑M1/45‑1 (MQ=255)                                                                        gttgtaatcaccgcatccccgatggcctgcagggtaatctgcgcacgctctttttcttcctgcagggcgctggcaaagGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCttt                                                                                                           >  1:327541‑M1/79‑123 (MQ=255)                                                                                             ttGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAg                                                           <  2:188636/150‑1 (MQ=255)                                                                                                           gctgATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGgaa                                                                                         >  1:65385/1‑106 (MQ=255)                                                                                                           gctgATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGgaa                                                                                         <  2:65385/106‑1 (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)                                                                                                                                                      GCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAAc  >  1:155778‑M1/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCGATTGGAAGGCCCGCATTGTCGCCATCAACGAAGGCGTGACCCCGGCGGCCAAGGCCTTCAGTGATGCCTTGGGTGAAGGCTCGCGCATGTTGCTGCGGGTGCTGCTGATCACCAACCTGCTGACCGCAATGTTCCTGATCGCCATTGCCTGGCGCCGCTCCAGCAAGCTGTTGGCCCAGCGCCAGGCCTTTGCCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAAC  >  NC_002947/1964415‑1964713

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.