New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1061357	30 (0.860)	4 (0.120)	4/252	3.9	12.3%	coding (2102/2547 nt)	torA	trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
	?	NC_000913	= 1061374	29 (0.890)					coding (2119/2547 nt)	torA	trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1061205‑1061357 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tttccgactCACATAACTGCGAGTGAAGTCGGAAATCCGGAT  <  NC_000913/1061374‑1061342                                                                                                                                                                                            GATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCC                                       >  6:98358/1‑149                                  GAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGACTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTCACATAACTGCGAGTGAAGTCGGAAATCC      >  4:109504/1‑149                                        CGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGACTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTCACATAACTGCGAGTGAAGTCGGAAATCCGGAT  <  3:109504/147‑1                                            GATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTCACATAACTGCGAGTGAAGTCGGAAATC       >  1:84896/1‑138                                            GATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTCACATAACTGCGAGTGAAGTCGGAAATC       <  2:84896/138‑1                                                     GGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGAC                                    <  5:157243/100‑1                                                     GGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGAC                                    >  6:157243/1‑100                                                                                                                                                                                            GATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1061205‑1061357 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tttccgactCACATAACTGCGAGTGAAGTCGGAAATCCGGAT  <  NC_000913/1061374‑1061342

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