New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4164667 =	34 (1.110)	7 (0.250)	3/122	2.7	21.7%	coding (1029/1845 nt)	btuB	vitamin B12/cobalamin outer membrane transporter
	?	NC_000913	4164716 =	19 (0.670)					coding (1078/1845 nt)	btuB	vitamin B12/cobalamin outer membrane transporter

CTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACG  >  NC_000913/4164660‑4164779                                                         |                                                                cTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGAcc                                                      >  2:952209/1‑68 (MQ=255)    cgacgCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCggg                                                   <  1:78870/68‑1 (MQ=255)      acgCCGGGTACAGGTTATGTTGAGGATGGATATGATCAGCGTAATACCGGCATCTATCTGACCGGGCt                                                 >  2:1127192/1‑68 (MQ=255)                  cgttgatcatatccatcctcaacataccctgtacccggcATCTATCTGACCGGGCTGCAACAAGTCgg                                     >  2:93369‑M2/40‑68 (MQ=255)                  cgttgatcatatccatcctcaacataacctgtacccggcATCTATCTGACCGGGCTGCAACAAGTCgg                                     >  2:94822‑M2/40‑68 (MQ=255)                   gttgatcatatccatcctcaacataacctgtacccggcATCTATCTGACCGGGCTGCAACAAGTCGGc                                    >  2:850426‑M2/39‑68 (MQ=255)                   gttgatcatatccatcctcaacataacctgtacccggcATCTATCTGACCGGGCTGCAACAAGTCGGc                                    >  2:200678‑M2/39‑68 (MQ=255)                   gTTATGTAGAGGATGGATAAGATCAACGTATTAAAGGCATCAATCAGACCGGGCTGCAACTAGTCgga                                    >  1:1044790/1‑67 (MQ=255)                    ttATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCg                                   >  2:569110/1‑68 (MQ=255)                       tGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGAtt                                >  2:73638/1‑68 (MQ=255)                         ttGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGCGCGGGCTGCAACAAGTCGGCGAttgt                              >  2:124322/1‑66 (MQ=255)                         ttGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGAtttt                              >  1:1058122/1‑68 (MQ=255)                         ttGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGAtttt                              >  1:177321/1‑68 (MQ=255)                         ttGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGAtttt                              >  2:1109156/1‑68 (MQ=255)                               atggatATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCttt                        >  1:422295/1‑68 (MQ=255)                                  gatATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGaa                     <  2:691590/68‑1 (MQ=255)                                    tatGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAgg                   >  1:863290/1‑68 (MQ=255)                                    tatGATCAACGTAATACCGGAATCTATCTGAACGGGCTGCAACAAGTCGGCGATTTTACCTTTGACgg                   >  1:1022584/1‑68 (MQ=255)                                     atGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGgc                  <  2:630196/68‑1 (MQ=255)                                         tCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCgcag              >  2:670901/1‑68 (MQ=255)                                         ataacctgtacccggcATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCgcag              <  1:200678‑M2/52‑1 (MQ=255)                                         ataacctgtacccggcATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCgcag              <  1:94822‑M2/52‑1 (MQ=255)                                         ataacctgtacccggcAGCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCgcag              <  1:93369‑M2/52‑1 (MQ=255)                                              gTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGc         >  1:503082/1‑68 (MQ=255)                                              gTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGc         >  1:605109/1‑68 (MQ=255)                                                 aTACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGt      <  2:863287/68‑1 (MQ=255)                                                    ccggcATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGAc   >  1:575338‑M2/6‑68 (MQ=255)                                                     cggcATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACg  <  2:575338‑M2/64‑1 (MQ=255)                                                     cggcATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACg  <  1:1127197‑M2/64‑1 (MQ=255)                                                         |                                                                CTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACG  >  NC_000913/4164660‑4164779

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