New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	2322337 =	5 (0.100)	36 (0.740)	29/156	0.5	88.2%	coding (139/963 nt)	menC	o‑succinylbenzoate synthase
	?	NC_012967	2772857 =	NA (NA)					intergenic (+378/‑54)	iap/ECB_RS13760	alkaline phosphatase isozyme conversion aminopeptidase/IS4‑like element IS421 family transposase

ATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAT  >  NC_012967/2772784‑2772928                                                                          |                                                                        cttttGCTTCATCGTGCTGCCCGGTTTATCCCCGCTGGCGCGGGGCACACCCTGAACGGTGGGGAAAAACCCCCTCAGCGCt                                                                 >  1:2325244/2‑82 (MQ=25)  ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                                <  1:1903327/82‑1 (MQ=255)          cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                        >  1:878019/1‑82 (MQ=255)          cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                        >  1:3832954/1‑82 (MQ=255)                gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                                  >  1:2413080/1‑82 (MQ=255)                 atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                                 <  1:1511114/82‑1 (MQ=255)                  tgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                >  1:4388564/1‑82 (MQ=255)                    aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                              >  1:130110/1‑82 (MQ=255)                     ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                             <  1:3278001/82‑1 (MQ=255)                                ccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                                  <  1:553724/82‑1 (MQ=255)                                       cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                           >  1:1597135/1‑82 (MQ=255)                                       cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGATGTGGTATTACGCCTGATATGATTTaa                           >  1:2522808/1‑82 (MQ=255)                                        gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                          >  1:3011795/1‑82 (MQ=255)                                        gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                          >  1:4134651/1‑82 (MQ=255)                                        gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                          >  1:3011187/1‑82 (MQ=255)                                         cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                         >  1:36444/1‑82 (MQ=255)                                            ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                      <  1:2775683/82‑1 (MQ=255)                                              aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                    <  1:1242347/82‑1 (MQ=255)                                                    cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt              <  1:2570490/82‑1 (MQ=255)                                                        aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct          >  1:2662355/1‑82 (MQ=255)                                                          ggTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct        >  1:3050896/1‑82 (MQ=255)                                                            tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa      >  1:97396/1‑82 (MQ=39)                                                            tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa      >  1:97397/1‑82 (MQ=39)                                                              gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg    <  1:1515172/82‑1 (MQ=37)                                                               ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGa   >  1:4368999/1‑82 (MQ=35)                                                                gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATAAATTACTCTCACGAt  >  1:4152542/1‑82 (MQ=18)                                                                          |                                                                        ATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAT  >  NC_012967/2772784‑2772928

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