New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3091097	30 (0.770)	5 (0.140)	5/238	4.1	15.9%	intergenic (+43/‑37)	endA/rsmE	DNA‑specific endonuclease I/16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
	?	NC_000913	= 3091103	25 (0.690)					intergenic (+49/‑31)	endA/rsmE	DNA‑specific endonuclease I/16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCGGGG  >  NC_000913/3090978‑3091244                                                                                                                        |                                                                                                                                                    gcgATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAtt                                                                                                                        >  7:53574/1‑149 (MQ=255)  cgATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAttt                                                                                                                       <  8:53574/149‑1 (MQ=255)   gATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAtttt                                                                                                                      >  1:195101/1‑149 (MQ=255)   gATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAtttt                                                                                                                      >  4:50332/1‑149 (MQ=255)    aTGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAttttt                                                                                                                     <  5:129924/149‑1 (MQ=255)       aaCGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTgg                                                                                                                                        >  4:132893/1‑127 (MQ=255)       aaCGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTgg                                                                                                                                        <  3:132893/127‑1 (MQ=255)                   aGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGt                                                                                                                                          <  3:22685/113‑1 (MQ=255)                   aGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGt                                                                                                                                          >  4:22685/1‑113 (MQ=255)                                tCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTAc                                                                                                                                                  >  3:103469/1‑92 (MQ=255)                                tCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTAc                                                                                                                                                  <  4:103469/92‑1 (MQ=255)                                                cAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAAACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCt                                                                                      >  8:13525/1‑136 (MQ=255)                                                cAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAAACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCt                                                                                      <  7:13525/136‑1 (MQ=255)                                                             aGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGc                                                                                                             >  1:187847/1‑100 (MQ=255)                                                             aGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGc                                                                                                             <  2:187847/100‑1 (MQ=255)                                                                          cTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTcac                                                              >  2:10818/1‑134 (MQ=255)                                                                          cTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTcac                                                              <  1:10818/134‑1 (MQ=255)                                                                          cTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCg                                                 >  2:41921/1‑147 (MQ=255)                                                                           tAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTAACATCGCGCTTTGCg                                                 >  5:51219/1‑146 (MQ=255)                                                                             aCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGa                                             <  1:41921/148‑1 (MQ=255)                                                                               cTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTa                                                                                             >  2:98891/1‑98 (MQ=255)                                                                                         gggATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCat                                <  3:50332/149‑1 (MQ=255)                                                                                           gATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGCCCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCatat                              >  2:124350/1‑149 (MQ=255)                                                                                                 ttttGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGgcg                        >  2:61749/1‑149 (MQ=255)                                                                                                      ttAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACACCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTAc                   >  1:37712/1‑149 (MQ=255)                                                                                                         aCCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTgc                <  2:195101/149‑1 (MQ=255)                                                                                                                    aCGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACACCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCg     <  2:37712/149‑1 (MQ=255)                                                                                                                        gTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCAc                                                                               >  5:103299/1‑71 (MQ=255)                                                                                                                        gTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCAc                                                                               <  6:103299/71‑1 (MQ=255)                                                                                                                        gTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCgggg  >  7:126259/1‑148 (MQ=255)                                                                                                                        |                                                                                                                                                    GCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCGGGG  >  NC_000913/3090978‑3091244

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