New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4370205 =	13 (0.440)	4 (0.140)	4/280	2.3	18.1%	coding (208/1257 nt)	yjeH	putative transporter
	?	NC_000913	4370316 =	24 (0.860)					coding (97/1257 nt)	yjeH	putative transporter
Rejected: Frequency below cutoff threshold.

TGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAA  >  NC_000913/4370078‑4370352                                                                                                                                |                                                                                                                                                    tGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCa                                                                                                                              >  2:222869/1‑151 (MQ=255)     cGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAAt                                                                                                                          <  2:368562/151‑1 (MQ=255)                       gTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACAcaatcgcaat                                                                                                        >  1:124655/1‑151 (MQ=255)                               gcaggcaAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAac                                                                                                 <  1:318868/150‑1 (MQ=255)                                           aCGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCa                                                                                    <  2:137695/151‑1 (MQ=255)                                                    aCCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGcc                                                                           <  1:35774/151‑1 (MQ=255)                                                        aTAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACg                                                                        >  2:374056/1‑150 (MQ=255)                                                                                 atcgacgtcattattaggcactggcgtgtttgccggtcctgcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATTGCAATCGGGGACACTAAGATAATCAAAACGGGCCACCCCCACAGGCTGTTATTTCCCGCTAc                                               <  1:3246‑M1/103‑1 (MQ=255)                                                                                           ttattaggcactggcgtgtttgccgttcctgcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTa                                    <  1:145285‑M1/114‑1 (MQ=255)                                                                                            aaCCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGGACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTa                                    <  2:159895/150‑1 (MQ=255)                                                                                                aaaCGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCa                               <  2:385636/151‑1 (MQ=255)                                                                                                   cGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGa                            >  2:402174/1‑151 (MQ=255)                                                                                                       ggcgtgtttgccgttcctgcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACgg                        <  1:116184‑M1/126‑1 (MQ=255)                                                                                                                  gTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCa             >  1:142890/1‑151 (MQ=255)                                                                                                                         gcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCta      <  2:159025‑M1/144‑1 (MQ=255)                                                                                                                            ccgccTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGAGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCtaat    <  1:449748/150‑1 (MQ=255)                                                                                                                             cgccTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCtaataa  <  2:41741/151‑1 (MQ=255)                                                                                                                                |                                                                                                                                                    TGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAA  >  NC_000913/4370078‑4370352

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