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.

AAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGT  >  NC_000913_3_bme_pgi/2729295‑2729587                                                                                                                                                     |                                                                                                                                                 aaaTATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGa                                                                                                                                                  >  1:286365/1‑149 (MQ=255)   atatatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGaca                                                                                                                                                <  2:620233/149‑1 (MQ=255)    tttttCAAACAAGACTCAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGacac                                                                                                                                               <  1:646103/145‑1 (MQ=255)        tCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGaaa                                                                                                                                           >  1:395521/1‑149 (MQ=255)           aaCAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAAtt                                                                                                                                        <  2:190260/149‑1 (MQ=255)                               tttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGCCACGAAAATTGCTTCTCACGCGTTGCGTGa                                                                                                                    >  1:320904/1‑149 (MQ=17)                                      tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                             <  2:399876/149‑1 (MQ=32)                                        ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgt                                                                                                           <  2:545043/149‑1 (MQ=32)                                                                                                                              tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                     <  1:431778/149‑1 (MQ=32)                                                                                                                                   ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                >  1:639297/1‑149 (MQ=35)                                                                                                                                    gCCTGTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGTTGAGCTTTCATTAATAGAAAGCGAACGGGCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc               <  2:378147/149‑1 (MQ=255)                                                                                                                                     ccTTTAGACGAAGGGGACACGAAAATTGCTGATCACGCGTTGCGCGATATTTTCGTATAGTGTGAGCTTTCATTAATAGAACGCGAACGGCCTTATTCTCTTCAGCCTCTCACCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca              <  2:587652/149‑1 (MQ=11)                                                                                                                                     ccTCTAGACGAAGGGGACACGAAAATTGCTTATCCCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca              <  1:250449/149‑1 (MQ=38)                                                                                                                                       tctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa            >  2:224902/1‑149 (MQ=37)                                                                                                                                           gACGAAGGGGACACGCAAATTGCTTATCACGCGTTGAGTGATATTTTCGTGTAGGGTCAGCTTGCCGTCATAGATAGCGAACGCGCTTGTTCTCTTAGGCCTCACTACCAACGCGTAAACGCCTTGCTATTCACTTTTCATACGACaac        >  2:288221/1‑148 (MQ=11)                                                                                                                                           gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt        >  2:672691/1‑149 (MQ=38)                                                                                                                                                 ggggACACGATAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt  <  1:342734/149‑1 (MQ=38)                                                                                                                                                     |                                                                                                                                                 AAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGT  >  NC_000913_3_bme_pgi/2729295‑2729587

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.