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

CGAGCCGAGGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGT  >  NC_002947/3981692‑3981963                                                                                                                                                      |                                                                                                                           tcgacacagccactgcgccaacgcgcgggaaggcgtcatcccgccagcccgaacggggacatggaacggtgcaagccaggagacctgcctcgtcacgttttcgactttcaaCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGt                                                                                                                            >  1:232838‑M1/112‑150 (MQ=255)                                 gcgtcatcccgccagcccgaacggggacatggaacggtgcaagccaggagacctgcctcgtcacgttttcgactttcaaCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAg                                                                                                  <  1:368003‑M1/65‑1 (MQ=255)                                 gcgtcatcccgccagcccgaacggggacatggaacggtgcaagccaggagacctgcctcgtcacgttttcgactttcaaCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAg                                                                                                  >  2:368003‑M1/80‑144 (MQ=255)                                        ggCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTg                                                                                     <  2:359346/150‑1 (MQ=255)                                                    cATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATg                                                                         >  2:29024/1‑150 (MQ=255)                                                      tGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTc                                                                       <  2:368000/150‑1 (MQ=255)                                                                   ttGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGc                                                                                                    <  1:359060/108‑1 (MQ=255)                                                                   ttGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGc                                                                                                    >  2:359060/1‑108 (MQ=255)                                                                                   ggCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTc                                          <  2:311428/150‑1 (MQ=255)                                                                                     ctgcctcgtcacgttttcgactttcaaccgggcggggtgatccggtggcgaacaagcccggccggTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGc                                              >  2:368003‑M2/66‑144 (MQ=255)                                                                                     ctgcctcgtcacgttttcgactttcaaccgggcggggtgatccggtggcgaacaagcccggccggTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGc                                              <  1:368003‑M2/79‑1 (MQ=255)                                                                                      aTGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTc                                                                             >  1:21432/1‑112 (MQ=255)                                                                                      aTGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTc                                                                             <  2:21432/112‑1 (MQ=255)                                                                                                            cAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCCCGTCGGCGCAGTGGCTgtgt                 >  1:384037/1‑150 (MQ=255)                                                                                                               accgggcggggtgatccggtggcgaacaagcccggccggTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGa              >  1:232838‑M2/40‑150 (MQ=255)                                                                                                                     cggGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGtt        <  1:434194/150‑1 (MQ=255)                                                                                                                        gCTTGTTCGCCGCCGGAGCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTc     <  2:400541/150‑1 (MQ=255)                                                                                                                           tGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGt  <  1:29024/150‑1 (MQ=255)                                                                                                                                     cGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCtgt                   >  2:114625/1‑123 (MQ=255)                                                                                                                                     cGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCtgt                   <  1:114625/123‑1 (MQ=255)                                                                                                                                                      |                                                                                                                           CGAGCCGAGGAAGCCGGTGACGATGGTGACGGGGAGCTTGGCCAGTGTTTTCATGTCGCATTGCCCTTGGCAGTGTGGCGCGGGCATGCAGGACGACAGCCGCGGGCCAGGCCGGCCGGGCTTGTTCGCCACCGGATCACCCCGCCCGGTTGAAAGTCGAAAACGTGACGAGGCAGGTCTCCTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGT  >  NC_002947/3981692‑3981963

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