New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4035244	66 (0.590)	4 (0.040)	3/268	NT	7.8%	intergenic (+91/‑287)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
	?	NC_000913	= 4166399	31 (0.290)					intergenic (+114/‑260)	murI/rrsB	glutamate racemase/16S ribosomal RNA of rrnB operon

TCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCAC  >  NC_000913/4166251‑4166514                                                                                                                                                     |                                                                                                                    tCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGg                                                                                                                    <  2:1054693‑M2/150‑2 (MQ=255) tCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGg                                                                                                                    >  1:1809615‑M2/1‑149 (MQ=255) tCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGg                                                                                                                    >  1:1255889‑M2/1‑149 (MQ=255)  cGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGCCACGGAACAACGGCAAACACGCCGCCGGgt                                                                                                                   <  2:655525‑M2/150‑3 (MQ=255)  cGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgt                                                                                                                   <  2:603294‑M2/150‑3 (MQ=255)  cGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgt                                                                                                                   <  1:1194003‑M2/150‑3 (MQ=255)   gAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgtc                                                                                                                  >  1:27188‑M2/1‑147 (MQ=255)   gAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgtc                                                                                                                  <  2:1734163‑M2/150‑4 (MQ=255)   gAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgtc                                                                                                                  >  1:1482297‑M2/1‑147 (MQ=255)       cGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCg                                                                                                              >  1:583399/1‑150 (MQ=255)        gCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCgg                                                                                                             <  1:977756/150‑1 (MQ=255)        gCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCAATATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCTGGTCAGCgg                                                                                                             >  2:175430/1‑150 (MQ=255)         cTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCggg                                                                                                            >  1:1252768/1‑150 (MQ=255)         cTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCggg                                                                                                            >  1:504657/1‑150 (MQ=255)           cGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGt                                                                                                          >  2:459716/1‑150 (MQ=255)           cGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGt                                                                                                          <  1:1393816/150‑1 (MQ=255)           cGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGt                                                                                                          >  2:19383/1‑150 (MQ=255)           cGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGt                                                                                                          <  1:724407/150‑1 (MQ=255)            gAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGtt                                                                                                         >  1:1394402/1‑150 (MQ=255)             aaaaaCTGGCAGTTTTAGGCTGATTTGGTTGACTGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGGCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTtc                                                                                                        <  1:1032248/150‑1 (MQ=255)                aaCTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTcc                                                                                                     <  1:1175157/150‑1 (MQ=255)                  cTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTg                                                                                                   <  1:821933/150‑1 (MQ=255)                    ggCAGTTTTAGGCTGCTTTGGTGGAATGTTGCGCGGTCAGAAAATAATTTTAAATTTCATCTTGTCAGTCCGGAATAACTACCTATAATGCGCCACCACTGACACGGCACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTgag                                                                                                 <  2:1077258/150‑1 (MQ=17)                     gCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTgaga                                                                                                >  2:1848465/1‑150 (MQ=255)                     gCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATAATTTTAAATTTCCTCTTCTCAGCCCGCAATAACTCCCTATAATGCGCCACCACTGACACCGAACAACGGCAACCACGCCGCCGGGTCAGCGGGGATCTCCTgaga                                                                                                >  2:1899610/1‑150 (MQ=17)                          tttAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCTGGGTCAGCGGGGTTCTCCTGAGAACTcc                                                                                           >  1:710029/1‑150 (MQ=255)                           ttAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCg                                                                                          >  2:645189/1‑150 (MQ=255)                           ttAGGCTGATTTGGTTGAATGTGGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCCGGCCGGAATAACTCCCTATAATGCGCCACCCCTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGAGTTCTCAGGAGAACCCCg                                                                                          <  1:1401681/150‑1 (MQ=14)                              ggCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGACTAACTCCCTATAATGCGCCACAACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCa                                                                                       <  2:1751633/150‑1 (MQ=255)                                cTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCaga                                                                                     >  1:1690195/1‑150 (MQ=255)                                 tGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCagag                                                                                    <  1:881616/150‑1 (MQ=255)                                 tGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCAACCACTGACACGGAACAACTGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCagag                                                                                    >  1:273680/1‑150 (MQ=255)                                  gATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAgaga                                                                                   <  2:212318/150‑1 (MQ=255)                                  gATTTGGTTGAATGTTGCGCGGGCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACTCGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAgaga                                                                                   <  1:1255304/150‑1 (MQ=255)                                    tttGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGaaa                                                                                 >  1:757856/1‑150 (MQ=40)                                       ggTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCa                                                                              >  2:1160401/1‑150 (MQ=38)                                       ggTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATATCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACCCGCCGCCGGGTCAGCCGGGTTCTCCTGAGAACGCCGGCAGAGAAAGCa                                                                              >  2:267114/1‑150 (MQ=255)                                        gTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACTCGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGACCTCCGGCAGAGAAAGCaa                                                                             >  2:376785/1‑150 (MQ=255)                                        gTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCaa                                                                             >  1:1622200/1‑150 (MQ=37)                                        gTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCaa                                                                             >  1:1560954/1‑150 (MQ=37)                                        gTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCATGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCaa                                                                             <  1:1787591/150‑1 (MQ=37)                                            aaTGTTTCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGTAAGCAAaaat                                                                         <  1:948626/150‑1 (MQ=17)                                            aaTGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAaaat                                                                         <  2:1290296/150‑1 (MQ=33)                                            aaTGGTGCGCGGGCAGAAAATTATTTTAAATTTCCTCTTGTCACGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACACCGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAaaat                                                                         <  2:77398/150‑1 (MQ=9)                                             aTGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAaaata                                                                        <  1:1074134/150‑1 (MQ=32)                                                                                                               gcCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcattttttattctttatgcttatttcagcgtcggtttgt      <  1:837123‑M2/150‑112 (MQ=255)                                                                                                                  accacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcattttttattctttatgcttatttcagcgtcggtttgtcgg   <  2:1742294‑M2/150‑115 (MQ=255)                                                                                                                  accacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcattttttattctttatgcttatttcagcgtcggtttgtcgg   <  1:528175‑M2/150‑115 (MQ=255)                                                                                                                   ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcattttttattctttatgcttatttcagcgtcggtttgtcggt  <  1:1195262‑M2/150‑116 (MQ=255)                                                                                                                                                     |                                                                                                                    TCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCAC  >  NC_000913/4166251‑4166514

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