New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2731302 =	44 (0.360)	8 (0.070)	6/196	14.6	15.5%	intergenic (‑145/+298)	rrsG/clpB	16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
	?	NC_000913	= 4035358	NA (NA)					intergenic (+205/‑173)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

GAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGG  >  NC_000913/2731212‑2731401                                                                                           |                                                                                                    gacaCTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCcgcttc                                                                                            <  2:67679/97‑1 (MQ=1)                                                               aGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGCGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTc                              <  1:50857/100‑1 (MQ=2)                                                                gATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                             <  1:921242/100‑1 (MQ=14)                                                                  ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                           >  1:2666104/1‑100 (MQ=18)                                                                  ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                           <  1:636069/100‑1 (MQ=17)                                                                   tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                          <  2:388180/100‑1 (MQ=18)                                                                      cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                       <  2:2783021/100‑1 (MQ=255)                                                                      cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                       >  2:136081/1‑100 (MQ=255)                                                                      cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                       >  1:846291/1‑100 (MQ=255)                                                                       tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                      >  1:1126106/1‑100 (MQ=255)                                                                       tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                      >  1:712010/1‑100 (MQ=255)                                                                       tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                      >  2:1202134/1‑100 (MQ=255)                                                                        gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTTTTTGAATGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                     >  1:533713/1‑100 (MQ=9)                                                                        gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGCGTCAAGCATTTCTTTGCtttt                     >  2:1786980/1‑100 (MQ=255)                                                                        gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                     >  2:1823534/1‑100 (MQ=255)                                                                        gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                     >  1:508455/1‑100 (MQ=255)                                                                        gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                     >  1:361477/1‑100 (MQ=255)                                                                          tAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct                   >  1:1141438/1‑100 (MQ=255)                                                                           aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                  >  1:1492645/1‑100 (MQ=255)                                                                           aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                  <  2:622320/100‑1 (MQ=255)                                                                            aaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                 >  2:2014876/1‑100 (MQ=255)                                                                             aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                <  2:802293/100‑1 (MQ=255)                                                                               tggtAAGGAGCAGTGCCGCTTTACTTTTTCTCTGCGGCGCGGGGTGTTAAGAATTCGCCTTTCCGCTACAGAGGCAAGCATTTCTTTGCTTTTCTCTGt               <  1:1890497/96‑1 (MQ=17)                                                                               tGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTCTCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtt              >  1:2184729/1‑100 (MQ=17)                                                                                tctAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATATGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg             <  1:906505/98‑1 (MQ=255)                                                                                gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCTTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg             <  1:2258765/100‑1 (MQ=255)                                                                                gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg             <  1:1168765/100‑1 (MQ=255)                                                                                gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg             >  1:269697/1‑100 (MQ=255)                                                                                gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg             >  2:781944/1‑100 (MQ=255)                                                                                 ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga            <  2:2717481/100‑1 (MQ=255)                                                                                  tAAAGAGCAGTGCCGCTTCACTTTTTCTCGGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACGGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag           <  1:1845954/100‑1 (MQ=17)                                                                                  tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag           <  2:285319/100‑1 (MQ=255)                                                                                   aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATATGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGGTgag           >  2:906508/1‑99 (MQ=255)                                                                                   aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga          >  1:2918030/1‑100 (MQ=255)                                                                                   aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga          >  1:2961591/1‑100 (MQ=255)                                                                                      gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc       <  1:503179/100‑1 (MQ=255)                                                                                       agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct      >  1:154605/1‑100 (MQ=255)                                                                                        gCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc     <  1:1823527/100‑1 (MQ=255)                                                                                        gCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc     <  1:399583/100‑1 (MQ=255)                                                                                         cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa    <  2:1399833/100‑1 (MQ=255)                                                                                         cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa    <  1:11662/100‑1 (MQ=255)                                                                                         cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa    <  1:1951995/100‑1 (MQ=255)                                                                                         cAGTGCCGCTTCACCTTTTCTCAGCGGCGCGGGGTGTGCATACTACGACTTCCCGCTACAGAGTCAAGCATTCCTTTCCTTGTCTCTGTTGAGATtctta    <  1:1786973/100‑3 (MQ=17)                                                                                          aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg   <  2:399525/100‑1 (MQ=255)                                                                                           gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg  >  2:2796528/1‑100 (MQ=255)                                                                                           gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg  <  2:834044/100‑1 (MQ=255)                                                                                           gTGCCGCTCCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg  <  1:712249/100‑1 (MQ=255)                                                                                           |                                                                                                    GAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGG  >  NC_000913/2731212‑2731401

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.