New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1243939 =	17 (0.610)	4 (0.150)	4/224	3.1	17.1%	coding (589/735 nt)	ycgR	flagellar velocity braking protein, c‑di‑GMP‑regulated
	?	NC_000913	1243983 =	23 (0.890)					coding (545/735 nt)	ycgR	flagellar velocity braking protein, c‑di‑GMP‑regulated
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGCACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCT  >  NC_000913/1243799‑1244084                                                                                                                                             |                                                                                                                                                  gcgcACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTg                                                                                                                                            <  4:126851/148‑1 (MQ=255)                                           gaAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGACTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAacac                                                                                                 >  6:16856/1‑149 (MQ=255)                                            aaaaTAATCCGCTGTAATTGCCGCTCCACCGTCGGACTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACAcc                                                                                                <  5:16856/149‑1 (MQ=255)                                                            aTTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGt                                                                                                          <  1:83354/123‑1 (MQ=255)                                                                               gCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGccc                                                                                       <  7:60331/123‑1 (MQ=255)                                                                               gCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGGTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGccc                                                                                       >  8:60331/1‑123 (MQ=255)                                                                                        aaGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAgc                                                             <  4:28979/140‑1 (MQ=255)                                                                                        aaGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAgc                                                             >  3:28979/1‑140 (MQ=255)                                                                                                  aaGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCAtt                                                                                           >  7:40272/1‑100 (MQ=255)                                                                                                  aaGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCAtt                                                                                           <  8:40272/100‑1 (MQ=255)                                                                                                     cTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTaa                                             >  7:61554/1‑143 (MQ=255)                                                                                                     cTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTaa                                             <  8:61554/143‑1 (MQ=255)                                                                                                                           tttCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGtt                                                                                 <  7:4937/85‑1 (MQ=255)                                                                                                                           tttCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGtt                                                                                 >  8:4937/1‑85 (MQ=255)                                                                                                                            ttCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGgag                                                                                                                             <  5:117864/40‑1 (MQ=255)                                                                                                                            ttCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGgag                                                                                                                             >  6:117864/1‑40 (MQ=255)                                                                                                                                           aatcaCTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCt  >  1:154866/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                  GCGCACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCT  >  NC_000913/1243799‑1244084

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