New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4578928	128 (1.620)	6 (0.080)	6/284	6.3	4.5%	coding (410/1380 nt)	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
	?	NC_000913	= 4579012	133 (1.750)					coding (326/1380 nt)	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4578775‑4578928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aggcaaCTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTCTATGCTGGACAACATAATCAACGACTATAAATTAATATTTAATTCTGGCAAGAG  <  NC_000913/4579012‑4578868                                                                                                                                                                                                                                                                                                             tGTATCTTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGA                                                                                                                                                        <  1:269288/148‑1  tTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATATAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATTACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAG                                                                                                                                                       <  1:1258881/148‑1  GTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAG                                                                                                                                                       >  2:1105890/1‑149    ATCGTCTCTTTTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGC                                                                                                                                                     <  1:580034/149‑1    ATCGTCTCTATTGTGGGTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTTTGATTCAGTTTTTGACATAGGTGGAATAACACTCATGCAAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGC                                                                                                                                                     <  2:772647/149‑1     TCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATATTCGTTGATTATGTTGTCCAGCATAGAGGCA                                                                                                                                                    >  2:527334/1‑149       GTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTTTCACGCATAGAGtctaa                                                                                                                                                  >  2:936883/1‑144           CTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAACTTCG                                                                                                                                              >  2:74503/1‑149            TATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAACTTCGG                                                                                                                                             >  2:748846/1‑149                                                                                          AACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGTATAGAGGCAA‑ATCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATT                                                              <  1:514778/149‑1                                                                                                CTTGCCAGAATTAAATATAAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCT                                                         <  1:146026/149‑1                                                                                                                   AATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTCTATGCTGGACAACATAAT                                      <  1:748846/149‑1                                                                                                                    ATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAACTTCGGGTGAATATCCTAAAAAATACCGACAGTCCTATTACGCCTGTTCCCAAAAAGACTCACAGGGTATTGATTACACCCGATTTGCCTCTATGCTGGACAACATAATC                                     <  1:74503/149‑1                                                                                                                                                     AGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTCTATGCTGGACAACATAATCAACGACTATAAATTAATATTTAATTCTGGCAAG    >  1:471920/1‑149                                                                                                                                                       GCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTCTATGCTGGACAACATAATCAACGACTATAAATTAATATTTAATTCTGGCAAGAG  >  2:192578/1‑149                                                                                                                                                                                                                                                                                                             AGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4578775‑4578928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aggcaaCTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTCTATGCTGGACAACATAATCAACGACTATAAATTAATATTTAATTCTGGCAAGAG  <  NC_000913/4579012‑4578868

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