New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3091309 =	11 (0.550)	4 (0.220)	3/118	1.8	25.1%	coding (176/732 nt)	rsmE	16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
	?	NC_000913	3091468 =	14 (0.790)					coding (335/732 nt)	rsmE	16S rRNA m(3)U1498 methyltransferase, SAM‑dependent

AATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATT  >  NC_000913/3091406‑3091535                                                               |                                                                    aaTGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgc                                                                <  2:513480/68‑1 (MQ=255)           gattcgcgatcgtcgatctggccttccagcaccttcacttccacgcttttttCTGAGCGCTGCGGCGt                                                      <  1:359761‑M2/16‑1 (MQ=255)           gattcgcgatcgtcgatctggccttccagcaccttcacttccacgcttttttCTGAGCGCTGCGGCGt                                                      <  1:573692‑M2/16‑1 (MQ=255)           aCTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                      >  2:650853/1‑68 (MQ=255)              aTCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTCTCTGAGCGCTgcggc                                                        <  2:443552/63‑1 (MQ=255)              aTCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTATCTGAGCGCTgcggc                                                        >  1:443551/1‑63 (MQ=255)                ccAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAAc                                                 >  2:643436/1‑68 (MQ=255)                ccAGAAATCGATCGAACTCGGTGTAAGCCGCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAAc                                                 <  2:3429/68‑1 (MQ=255)                  gatcgtcgatctggccttccagcaccttcacttccacgcttttttCTGAGCGCTgcggc                                                        >  2:197714‑M2/46‑59 (MQ=255)                  gatcgtcgatctggccttccagcaccttcacttccacgcttttttCTGAGCGCTgcggc                                                        <  1:197714‑M2/14‑1 (MQ=255)                       tcgatcgaACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAg                                          >  2:384701/1‑68 (MQ=255)                             gaACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACg                                    <  2:194485/68‑1 (MQ=255)                                          gCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaagaag                       <  1:643435/68‑1 (MQ=255)                                           ccTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGc                      >  1:494095/1‑68 (MQ=255)                                           ccTCATTACGCCACTTTTTTCTGAGAGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGc                      >  2:694673/1‑68 (MQ=255)                                                ttACGCCACTTTTTTCTGAGCGCTGCGGCGTTATACTGGATAGTGAACGTCTGAACAAGAAGCTTcag                 <  2:618701/68‑1 (MQ=255)                                                ttACGCCACTTTATTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcag                 <  2:618702/68‑1 (MQ=255)                                                     ccACTTTTTTCTGAGCGCTGCGGCTTTAATCTGTATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTg            <  1:679671/68‑1 (MQ=255)                                                              tCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAt   >  2:560557‑M2/2‑68 (MQ=255)                                                               cTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAtt  <  2:49757/68‑1 (MQ=255)                                                               |                                                                    AATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATT  >  NC_000913/3091406‑3091535

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