New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	2729443 =	17 (0.380)	5 (0.120)	5/266	5.4	18.2%	noncoding (2/76 nt)	gltW	tRNA‑Glu
	?	NC_000913_3_bme_pgi	= 4168331	29 (0.700)					intergenic (+131/‑41)	rrsB/gltT	16S ribosomal RNA of rrnB operon/tRNA‑Glu
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAG  >  NC_000913_3_bme_pgi/4168188‑4168421                                                                                                                                                |                                                                                           gATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGCTTGTCTGATAGCAAGTGAAAAGCAAGTCGTTTACGCGTTGGACGCGAGGCTGAATAAAAAAAGCCCGTTCTCTTTATATTAATGAAAACACACGCTgcaag                                                                                       >  2:34371/1‑147 (MQ=14)  aTCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGCTCGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCATTGGGAGTGAGGCTGAAGAGAACAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGa                                                                                      >  2:30486/1‑149 (MQ=255)   tCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGaa                                                                                     >  2:557158/1‑149 (MQ=255)     aCCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCGGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAGGGGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTGACCCGTCACGacaa                                                                                   >  2:4701/1‑149 (MQ=14)      cctcctTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGACGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCAGTTCGCTTTCTATTACTGAAAACTCACCCTACACGaaact                                                                                  >  2:306788/1‑147 (MQ=255)          cttaccttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                              >  1:728452/1‑149 (MQ=255)               cttaAAGAAGGTGCCTTTGTAGGGCTCACACAGATTGTCTGCTAGAAAGGTAAAAGCAAGGCGTTTACGCGTTGTGATTTAGTCCGAAGAGAACAACCCAGATCTCTTCCCATCACTGCACGCTCTCCCTACACGAAAATATCACGCaa                                                                         <  2:151808/149‑1 (MQ=9)                ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                        >  1:714305/1‑149 (MQ=255)                 taAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCCCGCAAcg                                                                       >  2:74802/1‑149 (MQ=255)                   aagaagCGGACTTTGTAGTGCTCACACAGATTGTCTGATAGAAACTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGACAATAAGGCCGTTCGCTTTCCATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAcgcg                                                                     <  2:624323/149‑1 (MQ=255)                    agaagCGTACTGTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGTTAGTGAGTCTGACGAGAATAAGCCAGATCGCCTTCCATTAATGAAAGCTTACCCTACACGAAAATATCACGCAACGCGt                                                                    <  2:773155/149‑1 (MQ=14)                     gaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTg                                                                   >  2:503019/1‑149 (MQ=255)                           tACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGc                                                             >  1:415303/1‑149 (MQ=255)                             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCaa                                                           <  2:728452/149‑1 (MQ=255)                              ttttTAGTGCTCACACCGATTGTCTGATACCAAGTGAAACGCAAGACGTTTACGCGATAGGAGTGATAATGGAGAGAATAAGGCCGTTCGCTTTCTATTAGTGTAAGCTCACCCTACACGTACATATCCCGCACCGCGTGCTactacat                                                          >  2:3673/1‑143 (MQ=2)                                   aGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCg                                                     >  2:648778/1‑149 (MQ=38)                                    gTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                    <  2:461547/149‑1 (MQ=38)                                    gTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                    <  1:159251/149‑1 (MQ=37)                                    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                                                    <  2:154186/119‑1 (MQ=32)                                                                  aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCgt                                                    >  1:154186/1‑119 (MQ=32)                                                                   aaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGc                     >  1:640595/1‑149 (MQ=31)                                                                                      tttggAGTGAGGCTGAAGAGAATAAGGCCGGTCGCTTTCTATTAATGAAAGCTCACCCTACCCGAACATATCACGCAACGCGTGATAAGCCATTTTCGTTGCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAg  <  1:274603/146‑1 (MQ=1)                                                                                                                                                |                                                                                           GATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAG  >  NC_000913_3_bme_pgi/4168188‑4168421

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