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.

CGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAA  >  NC_000913/3642796‑3643080                                                                                                                                               |                                                                                                                                               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                                                          <  1:24051/127‑1 (MQ=255)                                                                                                       ccACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTaa                                                          >  2:24051/1‑127 (MQ=255)                                                                                                                  ccGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTc                         >  4:63232/1‑149 (MQ=255)                                                                                                                          accGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGcccc                  >  4:48941/1‑148 (MQ=255)                                                                                                                          accGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATAGTGCTCCAGccccc                 >  5:2065/1‑149 (MQ=255)                                                                                                                            cgccgAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAg               <  2:63425/149‑1 (MQ=255)                                                                                                                            cgccgAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAg               >  1:63425/1‑149 (MQ=255)                                                                                                                               cgAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCgg            >  5:68454/1‑149 (MQ=255)                                                                                                                                            gTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGaa  >  7:50689/1‑146 (MQ=255)                                                                                                                                               |                                                                                                                                               CGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGTCCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAA  >  NC_000913/3642796‑3643080

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.