New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3870040	184 (0.970)	3 (0.020)	3/286	NT	1.6%	coding (664/1065 nt)	cbrA	colicin M resistance protein; FAD‑binding protein, putative oxidoreductase
	?	NC_000913	= 3870083	203 (1.110)					coding (707/1065 nt)	cbrA	colicin M resistance protein; FAD‑binding protein, putative oxidoreductase
Rejected: Frequency below/above cutoff threshold.

GAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3869761‑3870040 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaaagGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTCTTTCATTGGATAGGCACCGCCAAAGATAAAATAACCGTCTTTGCTGATACTCCATGAATAACAGTTAGTTATCGAATTATCAAAGATGCAGGAGTAGAAC  <  NC_000913/3870083‑3869940                                                                                                                                                                                                                                                                                                                                                                                                                                          GAGGATGATAAATGGCATGTCATTTTTCGTTCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGCCATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAAT                                                                                                                                                       <  2:462474/275‑1                                                                                                                                TCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATG                                                                                                                                                      <  5:58668/149‑1                                                                                                                                  CGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAA                                                                                                                                                    <  6:9013/149‑1                                                                                                                                   GTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAA                                                                                                                                                   <  7:731195/149‑1                                                                                                                                   GTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAA                                                                                                                                                   >  7:826681/1‑149                                                                                                                                   GTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAA                                                                                                                                                   >  8:758889/1‑149                                                                                                                                                                         CATCCGGTGCCGTTCTACTTCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAA                                                                                                                                                    >  5:67883/1‑110                                                                                                                                                                         CATCCGGTGCCGTTCTACTTCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAA                                                                                                                                                    <  6:67883/110‑1                                                                                                                                                                                TGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAGGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTC                                                                                                      <  6:694164/149‑1                                                                                                                                                                                                  TCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAGGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTCTTTCATTGGATA                                                                                          <  7:984406/143‑1                                                                                                                                                                                                  TCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAGGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTCTTTCATTGGATA                                                                                          >  8:984406/1‑143                                                                                                                                                                                                                             ATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATG                                                                                                                                                      >  4:80068/1‑56                                                                                                                                                                                                                                 ATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAG                                                                                                                                                  <  2:80068/56‑1                                                                                                                                                                                                                                                                                    GAAAGGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTCTTTCATTGGATAGGCACCGCCAAAGATAAAATAACCGTCTTTGCTGATACTCCATGAATAACAGTTAGTTATCGAATTATCAAAGATGCAGGAGTAGAAC  <  6:561919/149‑1                                                                                                                                                                                                                                                                                    GAAAGGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTCTTTCATTGGATAGGCACCGCCAAAGATAAAATAACCGTCTTTGCTGATACTCCATGAATAACAGTTAGTTATCGAATTATCAAAGATGCAGGAGTAGAAC  <  7:942485/149‑1                                                                                                                                                                                                                                                                                                                                                                                                                                          GAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3869761‑3870040 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaaagGCGCTCATTTTCTCTTTCAGCGTCGTGAAACGCGTCTGACCGTCTTTCATTGGATAGGCACCGCCAAAGATAAAATAACCGTCTTTGCTGATACTCCATGAATAACAGTTAGTTATCGAATTATCAAAGATGCAGGAGTAGAAC  <  NC_000913/3870083‑3869940

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