New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4607578	82 (0.960)	6 (0.070)	4/172	NT	6.8%	coding (123/1032 nt)	rsmC	16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
	?	NC_000913	= 4607601	84 (1.010)					coding (100/1032 nt)	rsmC	16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

CGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4607478‑4607578 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  <  NC_000913/4607601‑4607503                                                                                                                                                                                                          CGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCT                                                                                                      <  2:867659/100‑1   ACGTCATCCGACGTGGCGACCAGACTGAAACGGGCGTTAGCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG                                                                                                     <  2:182101/99‑1                            TGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAG                                                                           <  1:939041/100‑1                             GAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  2:880969/100‑1                              AAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGC                                                                          <  2:1665226/99‑1                               AACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGCCG                                                                        <  1:400677/100‑1                               AACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGCCG                                                                        <  1:406663/100‑1                                                ATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAAT                                                       <  1:857135/100‑1                                                                                                    TGCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAG   >  1:2084063/1‑100                                                                                                     GCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  >  1:1475980/1‑100                                                                                                     GCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  <  2:2069060/100‑1                                                                                                     GCCCGCGCGTTTAGAAACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAG   <  2:2245565/99‑1                                                                                                                                                                                                          CGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4607478‑4607578 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgCCCGCGCGTTTAGATACCGCGGCCAGCCGTGCTCATACCCAGCAATTCCACCACTGGCAGGTATTAAGCCGCCAGATGGGGGATAACGCCCGTTTCAGT  <  NC_000913/4607601‑4607503

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