New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 4168341	33 (0.740)	6 (0.140)	6/266	5.0	17.7%	intergenic (+141/‑31)	rrsB/gltT	16S ribosomal RNA of rrnB operon/tRNA‑Glu
	?	NC_000913_3_bme_pgi	= 4168363	25 (0.600)					intergenic (+163/‑9)	rrsB/gltT	16S ribosomal RNA of rrnB operon/tRNA‑Glu
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGAT  >  NC_000913_3_bme_pgi/4168216‑4168508                                                                                                                                                    |                                                                                                                                                  cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCaa                                                                                                                                                  <  2:728452/149‑1 (MQ=255)       aGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCg                                                                                                                                            >  2:648778/1‑149 (MQ=38)        gTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                                                                                                           <  1:159251/149‑1 (MQ=37)        gTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                                                                                                           <  2:461547/149‑1 (MQ=38)        gTGCTCACACAGATCGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                                                                                                           <  1:554852/149‑1 (MQ=38)          gCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgtgt                                                                                                                                         <  2:251041/149‑1 (MQ=38)            tCACACAGATTGGCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                                       <  2:225026/149‑1 (MQ=38)                  aGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCg                                                                                                                                 >  1:112342/1‑149 (MQ=39)                      tGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTa                                                                                                                             <  1:774195/149‑1 (MQ=38)                              aGAAAGTGAAAAGCAAGGCATGGACTCTGTTTGAGTGAGGGTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTGACCCTACACGAAAATATCACGGAACGCGTGATAAGCGATTTTCGTGTCCCCTGCGTCTAGAGGCCCa                                                                                                                     <  2:748629/149‑1 (MQ=11)                               gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAg                                                                                                                    <  2:653028/149‑1 (MQ=33)                                 aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTagag                                                                                                                          >  2:783443/1‑141 (MQ=32)                                 aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTagag                                                                                                                          <  1:783443/141‑1 (MQ=32)                                 aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGa                                                                                                                  <  1:692883/149‑1 (MQ=32)                                    tGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGacac                                                                                                               >  1:508710/1‑149 (MQ=32)                                    tGAAAAGCAAGGCGGTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCGATTAATGAAAGCTCACCCTACACGAAAAGATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAAAGGCCCAGGacac                                                                                                               >  2:683779/1‑149 (MQ=255)                                      aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                                                                                                           >  1:154186/1‑119 (MQ=32)                                      aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                                                                                                           <  2:154186/119‑1 (MQ=32)                                       aaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGc                                                                                                            >  1:640595/1‑149 (MQ=31)                                                          tttggAGTGAGGCTGAAGAGAATAAGGCCGGTCGCTTTCTATTAATGAAAGCTCACCCTACCCGAACATATCACGCAACGCGTGATAAGCCATTTTCGTTGCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAg                                                                                         <  1:274603/146‑1 (MQ=1)                                                                                                                      cacGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCGCAGGAGACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCCACCTTa                             >  2:329715/1‑149 (MQ=9)                                                                                                                         gAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata                          >  1:447167/1‑149 (MQ=32)                                                                                                                         gAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCATTAGTCAAAAGGCCAAGGACTCTGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAATCGCCGCCCTTAata                          >  1:385970/1‑149 (MQ=11)                                                                                                                              tatCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGACCTTTTACGGCGGTAACAGGGGTTCGAATCCCCTAGGGAACGCCACTTGCTGGTTTGTGAGGGAAAGTCGCCGACCTTAATGTCTCa                     >  1:222637/1‑149 (MQ=11)                                                                                                                                  aCGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAAc                 >  2:89359/1‑149 (MQ=32)                                                                                                                                  aCGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAAc                 >  2:89440/1‑149 (MQ=32)                                                                                                                                  aCGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTACGGGACGCCACTTGCTGGTTTGTGAATGAAAGTAGCCGACCTTAATATCTaaaacc                 >  2:761408/1‑147 (MQ=255)                                                                                                                                     cAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCTCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCa              >  1:519371/1‑149 (MQ=25)                                                                                                                                         gcgTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCtt          >  2:443692/1‑149 (MQ=38)                                                                                                                                            tGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCgg       <  2:503097/149‑1 (MQ=40)                                                                                                                                              aTAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGGGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGt     >  1:198431/1‑149 (MQ=255)                                                                                                                                                 aGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATCTCTCAACACTCCTCTTCGCGTgct  >  2:473929/1‑147 (MQ=255)                                                                                                                                                 aGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGAt  >  1:431954/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGAT  >  NC_000913_3_bme_pgi/4168216‑4168508

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