New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4607286	107 (0.910)	5 (0.040)	5/268	NT	4.5%	coding (415/1032 nt)	rsmC	16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
	?	NC_000913	= 4607339	107 (0.940)					coding (362/1032 nt)	rsmC	16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below/above cutoff threshold.

AGCAGCAACTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4607135‑4607286 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgCCGAGCAGATGCTGGCAGATTATGCGCCGTTGAATAAAGTCGACAGCGCTCGTCGCTGTGGCCTCTATTTTGGTCGTCTGGAAAAACAGCCGGTATTTGATGCGGAGAAATTCTGGGGCGAATACAGCGTCGATGGCCTGACGGTCA  <  NC_000913/4607339‑4607193                                                                                                                                                                                                                                                                                                             AGCAGCAACTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTGTTTCCAGACGCCCAAAATAGAGGCCACAG                                                                                                                                                       >  1:1788687/1‑150  GCAGCAACTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGC                                                                                                                                                      <  1:1683673/150‑1  GCAGCAACTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGC                                                                                                                                                      <  2:703334/150‑1   CAGCAACTGGCTACCGACATCCGGACCGTAGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGCCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGGTTTTACAGACGACCAAAATAGAGGCCACAGCG                                                                                                                                                     >  2:1372236/1‑150   CAGCAACTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCG                                                                                                                                                     >  1:1326241/1‑150         CTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGCCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGTCCACAGCGCCGAGC                                                                                                                                               >  2:223227/1‑150           GGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTATTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGCCGAGCgc                                                                                                                                             >  2:426657/1‑148                                                           TGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGCCGAGCAGATGCTGGCAGATTATGCGCCGTTGAATAAAGTCGACAGCGCTCGTCGC                                                                                             <  1:96965/150‑1                                                               CGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGCCGAGCAGATGTTGGCAGATTATGCGCCGTTGAATAAAGTCGACAGCGCTCGTCGCTGTG                                                                                         <  1:173828/150‑1                                                                           GACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGGTGTTCCAGACGACCAAAATAGAGGCCGCAGCGCCGAGCAGAAGCTGGCAGATTATGCGCCGTTGAATAAAGGCGACAGCGCTCGTCGCTGTGGCCTCTATTTTG                                                                             <  1:27260/150‑1                                                                                                                                                      GCGCCGAGCAGATGCTGGCAGATTATGCGCCGTTGAATAAAGTCGACAGCGCTCGTCGCTGTGGCCTCTATTTTGGTCGTCTGGAAAAACAGCCGGTATTTGATGCGGAGAAATTCTGGGGCGAATACAGCGTCGATGGCCTGACGGTCA  >  1:787966/1‑150                                                                                                                                                                                                                                                                                                             AGCAGCAACTGGCTACCGACATCCAGACCGTCGCGGCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4607135‑4607286 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgCCGAGCAGATGCTGGCAGATTATGCGCCGTTGAATAAAGTCGACAGCGCTCGTCGCTGTGGCCTCTATTTTGGTCGTCTGGAAAAACAGCCGGTATTTGATGCGGAGAAATTCTGGGGCGAATACAGCGTCGATGGCCTGACGGTCA  <  NC_000913/4607339‑4607193

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.