New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3642916 =	31 (0.940)	3 (0.090)	3/268	4.4	8.6%	coding (217/753 nt)	rsmJ	16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
	?	NC_000913	3642938 =	34 (1.070)					coding (195/753 nt)	rsmJ	16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTC  >  NC_000913/3642779‑3643057                                                                                                                                          |                                                                                                                                              aCGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAAt                                                                                                                                     <  3:38004/148‑1 (MQ=255)  cGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAAtt                                                                                                                                    <  2:74581/148‑1 (MQ=255)                  cgcgGCAGCCGACTGAAACCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATc                                                                                                                   <  8:53303/149‑1 (MQ=255)                             aCTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCcaa                                                                                                         >  2:4787/1‑148 (MQ=255)                                 aaGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACaaaa                                                                                                    <  4:74196/149‑1 (MQ=255)                                    ccAGTACAAAGGCATCTCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGTCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATc                                                                                                  >  8:84533/1‑148 (MQ=255)                                                        ccAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGttt                                                                             >  5:3566/1‑149 (MQ=255)                                                              ccTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtc                                                                       <  6:3566/149‑1 (MQ=255)                                                               cTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtca                                                                      >  5:123755/1‑149 (MQ=255)                                                                  cAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGAGACCACCGCCGAATTTGCGTCGGTACGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcatca                                                                   >  2:84085/1‑149 (MQ=255)                                                                     tGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgcc                                                                                      <  6:86645/127‑1 (MQ=255)                                                                     tGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgcc                                                                                      >  5:86645/1‑127 (MQ=255)                                                                       gCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCtt                                                              <  1:4787/149‑1 (MQ=255)                                                                            caccacATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTgcgc                                                          >  8:79356/1‑148 (MQ=255)                                                                                  aTCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCa                                                   <  7:79356/149‑1 (MQ=255)                                                                                   tCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTc                                                     >  3:65662/1‑146 (MQ=255)                                                                                     cGGCAAATAATCGCCTTTACTGCCCCCCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGAt                                                <  8:142187/149‑1 (MQ=255)                                                                                     cGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCCCCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCACCAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGAt                                                >  5:27700/1‑149 (MQ=255)                                                                                     cGGCAAATAATCGCCTTTAATGCCCACCGCGTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGAt                                                <  5:2774/149‑1 (MQ=255)                                                                                              atCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCg                                       >  8:14345/1‑149 (MQ=255)                                                                                               tCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCTCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGt                                      <  7:84533/149‑1 (MQ=255)                                                                                                  ccTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTaa                                   >  7:36336/1‑149 (MQ=255)                                                                                                    tttAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAaca                                 >  3:108209/1‑149 (MQ=255)                                                                                                             cACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCa                        >  7:28250/1‑149 (MQ=255)                                                                                                             cACCGCTTTCGCCACCGCCTCACCGCGACCACCGCAGAATTTGCGTCGGTGAGCCATCGCTCCGCCAACAAAATAAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCATATGTTCCGGCGTTAACACCAGCGCCa                        >  8:119523/1‑149 (MQ=255)                                                                                                              aCCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCGCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCAt                       >  3:19957/1‑149 (MQ=255)                                                                                                                     tCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGtt                 <  4:19957/148‑1 (MQ=255)                                                                                                                     tCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGtt                 <  8:28250/148‑1 (MQ=255)                                                                                                                        ccACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTaa                                   >  2:24051/1‑127 (MQ=255)                                                                                                                        ccACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTaa                                   <  1:24051/127‑1 (MQ=255)                                                                                                                                   ccGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTc  >  4:63232/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                              ACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTC  >  NC_000913/3642779‑3643057

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