New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4389180	124 (0.960)	3 (0.020)	3/264	NT	2.6%	coding (191/3324 nt)	mscM	mechanosensitive channel protein, miniconductance
	?	NC_000913	= 4389203	111 (0.900)					coding (168/3324 nt)	mscM	mechanosensitive channel protein, miniconductance
Rejected: Frequency below/above cutoff threshold.

GGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4389025‑4389180 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgcgctTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGACGCGC  <  NC_000913/4389203‑4389055                                                                                                                                                                                                                                                                                                                   GGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGA                                                                                                                                                             >  2:1985060/1‑150     GCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCG                                                                                                                                                         >  1:1546262/1‑150      CTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGC                                                                                                                                                        <  2:1550487/150‑1                                       GTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCA                                                                                                                       <  1:707872/150‑1                                          GACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAAC                                                                                                                    <  2:1193441/150‑1                                              TTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGAAAATTATCGATAACTTGCTGATATTGTTTGATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATA                                                                                                                <  1:1755853/150‑1                                                                                                                                                      atgagaTTGAGGAACGAAAAGTTTCACTTTAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGATCCGCGCAGCGTGTCGCCGGGAATGTCTACCG        <  1:1285603/144‑1                                                                                                                                                        GCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGcc      >  2:1396130/1‑148                                                                                                                                                        GCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGAC      <  2:1550789/150‑1                                                                                                                                                        GCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGAC      >  2:2000142/1‑150                                                                                                                                                        GCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACAGAC      >  2:189685/1‑150                                                                                                                                                         CGCTTGAGTAACGAAAAGGTGCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGACACTCACCGCTACTCTCCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGACG     <  2:1372252/150‑1                                                                                                                                                         CGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGACG     >  1:216852/1‑150                                                                                                                                                            TTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGACGCGC  <  1:2088142/150‑1                                                                                                                                                                                                                                                                                                                   GGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4389025‑4389180 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgcgctTGAGGAACGAAAAGGTTCCCTTGAGCGCATCAAACAATATCAGCAAGTTATCGATAATTATCCGAAACTCTCCGCTACTCTGCGCGCACAATTAAACAACATGCGTGACGAGCCGCGCAGCGTGTCGCCGGGAATGTCTACCGACGCGC  <  NC_000913/4389203‑4389055

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.