New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 302638	11 (0.940)	3 (0.260)	3/304	1.0	19.4%	intergenic (‑65/‑353)	paoA/yagU	aldehyde dehydrogenase, Fe‑S subunit/DUF1440 domain‑containing inner membrane protein YagU
	?	NC_000913	= 1924823	14 (1.200)					coding (147/375 nt)	yobA	copper‑binding protein
Rejected: Frequency below/above cutoff threshold.

CAGGTAAACGCCAGCATGATTTATTTCACGCTAAAGGTGTAGTGTCCTTTCGTTTTATGCCCGTCCACCGAAACAACATGCCAGTCAACGGTGTAGGTTCCGGGTTTCAGCGAATCGGCCAGTGGAACAATCAGTTGCTTTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTGAGTCGAGGATCATCA  >  NC_000913/1924575‑1925042                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            cAGGTAAACGCCAGCATGATTTATTTCACGCTAAAGGTGTAGTGTCCTTTCGTTTTATGCCCGTCCACCGAAACAACATGACAGTCAACGGTGTAGGTTCCGGGTTTCAGCGAATCGGCCAGTGGAACAATCAGTTGCTTTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGt                                                                                                                                                                                                                           >  1:43333/1‑251 (MQ=255)             aGCATGATTTATTTCACGCTAAAGGTGTAGTGTCCTTTCGTTTTATGCCCGTCCACCGAAACAACATGCCAGTCAACGGTGTAGGTTCCGGGTTTCAGCGAATCGGCCAGTGGAACAATCAGTTGCTTTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTGTGCGCCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCgn                                                                                                                                                                                                               <  1:54586/251‑2 (MQ=255)                                                 ttCGTTTTATGCCCGTCCACCGAAACAACATGCCAGTCAACGGTGTAGGTTCCGGGTTTCAGCGAATCGGCCAGTGGAACAATCAGTTGCTGTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTTATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATAc                                                                                                                                                                             <  1:133163/249‑1 (MQ=255)                                                             ccGTCCACCGAAACAAAATGCCAGTCAACGGTGTAGGTTCCGGGTTTCAGCGACTCGGCCAGTGGAACAATCAGTTGCTTTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATg                                                                                                                                                               >  1:82574/1‑251 (MQ=255)                                                               gTCCACCGAAACAACATTCCAGTAAACGGTGTAGGCTCCGGGTTTCAGCGCATCGGCCAGTGGAACAATCAGCTGCTTATGGTCCTGCTCATTTCGCTTCGCCGGTACTGTTTTAATATTTTCGTTTTTTGGACCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTAGCCGACGGATACTGATGCGTTAAATgcg                                                                                                                                                             >  1:100101/1‑251 (MQ=255)                                                                                                                        cAGTGGAACAATCAGTTGCTTTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGGTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCGGTTTCAATACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTTATGCGTTAAAGGCGCATGTGCCCAGACGGAATGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGc                                                                                                    <  1:41483/251‑1 (MQ=255)                                                                                                                                            nntGNNNNNGCTNATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGttt                                                                                                                                                                                         <  2:167217/144‑1 (MQ=255)                                                                                                                                                                      ngtANTNNTNNAATNTTTTCGTTTTTTGGCCCCGTGNNTTTTGCACCNCTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAtaaaaaagcatcagcacnctcgcaacgnnctnttatgntnnncnnnagnataccggaggcctggtatgngcnnccaaggngaatcccccgaagataatcgggttgg                                                                                                                   >  2:64182‑M2/2‑84 (MQ=255)                                                                                                                                                                      ggTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAtaaaaaagcatcagcactctcgcaacgcactcttattttcccctttagaataccggaggcctggtatgagcaaccaaggcgaataccccgaagataatcgggttgg                                                                                                                   <  1:64182‑M2/190‑107 (MQ=255)                                                                                                                                                                       gTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATCTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAg                                                                                                                                                                                                          <  1:20722/102‑1 (MQ=255)                                                                                                                                                                                          gTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGTTGTCACCAACGAAGTGGTCAGGATTTCCATCGCGt                                                                                                       <  1:227572/182‑1 (MQ=255)                                                                                                                                                                                                        tGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAATGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGGGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAAACTTGTCTAATGt                    >  1:139629/1‑251 (MQ=255)                                                                                                                                                                                                                   ncacTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTAGTCTAATGTCTGAGTCGAgg         >  1:223152/2‑251 (MQ=255)                                                                                                                                                                                                                        gAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCTCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTg                 >  1:133498/1‑238 (MQ=255)                                                                                                                                                                                                                          aTCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAATCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTGAGTCGAGGAtcatca  <  1:270099/251‑1 (MQ=255)                                                                                                                                                                                                                              tgtgtCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATAc                                                                                                                                                                             <  1:130872/72‑1 (MQ=255)                                                                                                                                                                                                                                                       gtataaaaaagcatcagcactctcgcaacgcactcttattttccccttgagaataccggaggcctggtatgag                                                                                                                                                       >  1:133046‑M2/3‑3 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            CAGGTAAACGCCAGCATGATTTATTTCACGCTAAAGGTGTAGTGTCCTTTCGTTTTATGCCCGTCCACCGAAACAACATGCCAGTCAACGGTGTAGGTTCCGGGTTTCAGCGAATCGGCCAGTGGAACAATCAGTTGCTTTTGGTCCTGCTCATTTCGCTTCGCCGGTAATGTTTTAATATTTTCGTTTTTTGGCCCCGTGATTTTTGCACCACTGAATCCTGTTTCAACACCTTCCGAGAAGTTTAAAGTGATTGCCTGCGGTGCAGCTGTCACTTGCGCGTTTGCCGCAGGATACTGATGCGTTAAATGCGCATGTGCCCAGACGGAAGGTGTCACCAACGAAGTGGTCAGGATTGCCAGCGCGTAGCGAAGGGAGCGTGCAGTTGAAGCCATATTATCTATTCCTTTTTGTAATAACTTTTTTACAGAGCATAACCTTGTCTAATGTCTGAGTCGAGGATCATCA  >  NC_000913/1924575‑1925042

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