New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4637531	76 (0.780)	3 (0.030)	3/162	NT	3.9%	coding (836/1425 nt)	creC	sensory histidine kinase in two‑component regulatory system with CreB or PhoB
	?	NC_000913	= 4637546	80 (0.900)					coding (851/1425 nt)	creC	sensory histidine kinase in two‑component regulatory system with CreB or PhoB
Rejected: Frequency below cutoff threshold.

GAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4637425‑4637531 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcggaCCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTCCCT  <  NC_000913/4637546‑4637448                                                                                                                                                                                                                GAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAG                                                                                                            >  1:615390/1‑100  AAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGC                                                                                                           <  1:211327/100‑1      TATGCGCGTGAAGCTGGAAGGGAAAATCTATATTGAGCAGTATGTTTACGTATTAACTCATGAGCTAACAAGCCCACTGGCGGCGATTCGTGGAGCGGCG                                                                                                       <  2:1661441/100‑1      TATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCG                                                                                                       <  1:1041109/100‑1                             AAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGACCTTCGCGTAAAATTTCCGCC                                                                                <  1:2063733/100‑1                               AACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGACCTTCGCGTAAAATTTCCGCCGC                                                                              <  1:1897689/100‑1                                    TATTGAGCAGGATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGACCTTCGCGTAAAATTTCCGCCGCTCCAC                                                                         <  1:2608750/100‑1                                                                                                       GCGGACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTT      <  1:2565946/100‑1                                                                                                        CGGACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTC     <  2:2278782/100‑1                                                                                                        CGGACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTC     <  1:2094160/100‑1                                                                                                        CGGACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTC     <  1:1458935/100‑1                                                                                                        CGGACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTC     <  1:1458592/100‑1                                                                                                         GGACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTCC    <  1:1071295/100‑1                                                                                                          GACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTCCC   <  1:1757007/100‑1                                                                                                          GACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTCCC   <  2:2143697/100‑1                                                                                                           ACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTCCCT  <  2:2160154/100‑1                                                                                                           ACCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTCAATGCGTAAACATACTGCTCAATATAGTTTTTC     <  2:1754647/97‑1                                                                                                                                                                                                                GAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4637425‑4637531 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcggaCCTTCGCGTAAAATTTCCGCCGCTCCACGAATCGCCGCCAGTGGGCTTTTTAGCTCATGAGTTAATGCGTAAACATACTGCTCAATATAGTTTTTCCCT  <  NC_000913/4637546‑4637448

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