New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3941549 =	19 (0.220)	9 (0.110)	3/168	NT	33.4%	intergenic (‑222/‑259)	yieP/rrsC	putative transcriptional regulator/16S ribosomal RNA of rrnC operon
	?	NC_000913	4035245 =	NA (NA)					intergenic (+92/‑286)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

CGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAA  >  NC_000913/3941450‑3941648                                                                                                    |                                                                                                    cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     >  2:928256‑M1/100‑100 (MQ=255) cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     >  2:63200‑M1/100‑100 (MQ=255) cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     >  1:1826019‑M1/100‑100 (MQ=255) cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     <  1:2384189‑M1/1‑1 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   >  2:657825‑M1/98‑100 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   <  2:1727637‑M1/3‑1 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   >  1:1561257‑M1/98‑100 (MQ=255)    cttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCa                                                                                                  <  1:1178503‑M1/4‑1 (MQ=255)    cttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCa                                                                                                  <  2:135634‑M1/4‑1 (MQ=255)    cttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCa                                                                                                  <  1:858457‑M1/4‑1 (MQ=255)    cttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCa                                                                                                  <  1:649775‑M1/4‑1 (MQ=255)     ttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCAg                                                                                                 >  1:1318885‑M1/96‑100 (MQ=255)     ttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCAg                                                                                                 >  1:913117‑M1/96‑100 (MQ=255)                                                            gaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGt                                          <  1:357815‑M1/60‑1 (MQ=255)                                                              accccgctgacccggcggcgtgtttgccgttgttccgtGTCAGTGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAg                                        <  1:354602‑M1/62‑1 (MQ=255)                                                               ccccgctgacccggcggcgtgtttgccgttgttccgtGTCAGTGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGc                                       <  1:2141064‑M1/63‑1 (MQ=255)                                                               ccccgctgacccggcggcgtgtttgccgttgttccgtGTCAGTGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGc                                       <  1:2274170‑M1/63‑1 (MQ=255)                                                               ccccgctgacccggcggcgtgtttgccgttgttccgtGTCAGTGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGc                                       <  1:956739‑M1/63‑1 (MQ=255)                                                               ccccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGc                                       <  1:1797007‑M1/63‑1 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGTGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCg                                      <  1:1956259‑M1/64‑1 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCg                                      <  1:963145‑M1/64‑1 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCg                                      <  2:2058030‑M1/64‑1 (MQ=255)                                                                        cacGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCg                              >  2:2049181/1‑100 (MQ=33)                                                                          cGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGta                            >  1:1190025/1‑100 (MQ=33)                                                                          cGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTAACTCTGTAGCGGGAAAGCg                              >  2:2003535/1‑98 (MQ=17)                                                                           ggAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtat                           >  1:38329/1‑100 (MQ=33)                                                                           ggAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGgat                           >  1:1035791/1‑100 (MQ=18)                                                                            gAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatt                          <  2:663553/100‑1 (MQ=33)                                                                               caacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATg                       >  1:955776/1‑100 (MQ=32)                                                                               caacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATg                       >  2:691179/1‑100 (MQ=33)                                                                                aacGGCAAACACGCCGCCGGGTCATCGGAGTTCTCCTGAGAACTCCGTCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGc                      <  1:2380062/100‑1 (MQ=255)                                                                                aacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGc                      <  1:1492686/100‑1 (MQ=32)                                                                                aacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGc                      <  1:2290101/100‑1 (MQ=33)                                                                                aacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGc                      <  1:1759004/100‑1 (MQ=33)                                                                                 acGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGca                     >  2:2298155/1‑100 (MQ=33)                                                                                  cGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcac                    >  2:463117/1‑100 (MQ=33)                                                                                   ggCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcaca                   >  1:2313435/1‑100 (MQ=33)                                                                                    gCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCacac                  >  1:1439648/1‑100 (MQ=33)                                                                                      aaaCACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAccc                <  1:810145/100‑1 (MQ=33)                                                                                       aaCACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAcccc               <  1:2049199/100‑1 (MQ=33)                                                                                         cacGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCcgc             >  2:2043667/1‑100 (MQ=32)                                                                                            gccgccGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCG            >  2:2144025‑M1/1‑98 (MQ=255)                                                                                                    gTCGGCGGTGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCTTATTATGCACCCCCCGCGCCGtt       <  2:290445‑M1/95‑3 (MQ=255)                                                                                                    GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGaa  <  2:1138860‑M1/100‑1 (MQ=255)                                                                                                    |                                                                                                    CGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAA  >  NC_000913/3941450‑3941648

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