New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3981803 =	14 (0.750)	3 (0.170)	3/238	NT	17.9%	intergenic (‑58/+39)	cobW/PP_mr43	adenosylcobalamin biosynthesis protein CobW/Cobalamin regulatory region
	?	NC_002947	3981841 =	14 (0.770)					intergenic (‑96/+1)	cobW/PP_mr43	adenosylcobalamin biosynthesis protein CobW/Cobalamin regulatory region

GCGCAGCAAGGTGGTCTTGCCCGAGCCGAGGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGT  >  NC_002947/3981671‑3981948                                                                                                                                     |                                                                                                                                                  gcgcAGCAAGGTGGTCTTGCCCGAGCCGAGGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGc                                                                                                                                  >  1:304906/1‑150 (MQ=255)            tggtCTTGCCCGAGCCGAGGCAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATCTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCCTGTTCGCCCCCGGATCAc                                                                                                                       >  2:317487/1‑150 (MQ=255)                   gCCCGAGCCGAGGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACcccgccc                                                                                                                <  2:304906/150‑1 (MQ=255)                     ccgagccgagGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCgg                                                                                                              >  2:434194/1‑150 (MQ=255)                      tcgacacagccactgcgccaacgcgcgggaaggcgtcatcccgccagcccgaacggggacatggaacggtgcaagccaggagacctgcctcgtcacgttttcgactttcaaCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGt                                                                                                             >  1:232838‑M1/112‑150 (MQ=255)                                 aGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACcccgcccg                                                                                                               >  2:124790/1‑137 (MQ=255)                                 aGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACcccgcccg                                                                                                               <  1:124790/137‑1 (MQ=255)                                                      gcgtcatcccgccagcccgaacggggacatggaacggtgcaagccaggagacctgcctcgtcacgttttcgactttcaaCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAg                                                                                   >  2:368003‑M1/80‑144 (MQ=255)                                                      gcgtcatcccgccagcccgaacggggacatggaacggtgcaagccaggagacctgcctcgtcacgttttcgactttcaaCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAg                                                                                   <  1:368003‑M1/65‑1 (MQ=255)                                                             ggCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTg                                                                      <  2:359346/150‑1 (MQ=255)                                                                         cATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATg                                                          >  2:29024/1‑150 (MQ=255)                                                                           tGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTc                                                        <  2:368000/150‑1 (MQ=255)                                                                                        ttGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGc                                                                                     >  2:359060/1‑108 (MQ=255)                                                                                        ttGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGc                                                                                     <  1:359060/108‑1 (MQ=255)                                                                                                        ggCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTc                           <  2:311428/150‑1 (MQ=255)                                                                                                           aTGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTc                                                              <  2:21432/112‑1 (MQ=255)                                                                                                           aTGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTc                                                              >  1:21432/1‑112 (MQ=255)                                                                                                                                 cAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCCCGTCGGCGCAGTGGCTgtgt  >  1:384037/1‑150 (MQ=255)                                                                                                                                     |                                                                                                                                                  GCGCAGCAAGGTGGTCTTGCCCGAGCCGAGGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGT  >  NC_002947/3981671‑3981948

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