New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1743497 =	9 (0.280)	9 (0.280) +G	8/486	1.4	50.1%	coding (41/1374 nt)	mdtK	multidrug efflux pump MdtK
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

CTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGC  >  NC_000913/1743251‑1743721                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 gggctgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcccTAGc                                                                                                                                                                                                                               <  2:152911‑M1/4‑1 (MQ=255)                gataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcatcttcccTAGCAATCCCGGTGATtctc                                                                                                                                                                                                               >  1:159055‑M1/232‑251 (MQ=255)                        cTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTCCCAGTGCAGGGGAAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAACTAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTattagcattagggaaggtgcgaacaagtccctgat                                                                                                                                                                                                        >  2:215617/1‑226 (MQ=255)                                       gCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGAt                                                                                                                                                                                        <  1:6451/251‑1 (MQ=255)                                             ggccaacctgtttcgggcggaccaaatgatacatcagtgggagagctctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcacggcatttacgggagaaaaaatcggctcaaaaatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcccTAGCAATCCCGGTGATTCTCGCGCAACTCGCCCAAACTGCGATGGGttt                                                                                                                                                                                  <  2:238898‑M1/49‑1 (MQ=255)                                                                           aaTGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCtat                                                                                                                                                      >  1:116605/1‑249 (MQ=255)                                                                              ggCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTg                                                                                                                                                  <  1:3199/250‑1 (MQ=255)                                                                                                     aataaATCCTCGCTACAATAGACTGACTTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGcc                                                                                                                           >  2:148588/1‑249 (MQ=255)                                                                                                              tCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACtt                                                                                                                 >  1:131119/1‑251 (MQ=255)                                                                                                                                    acggtctaaataggccgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcacctgcccTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCCCCGACATGGCGGCGTTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCct                                                                                           >  1:3623‑M1/116‑251 (MQ=255)                                                                                                                                     ccTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCctc                                                                                          <  2:188820/251‑1 (MQ=255)                                                                                                                                             ataggctgattaaaggcatttacgggagaaaaaatcggctcaaacttgaagaaatgaaatgactgagtaagccgagaagaatttccccgattcttcgcaccttcctTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGCTACCGTGATGGCGGGCGGCTATAGTGCCACCGCCATTGCGGCGGACGACATCGGTACTTATATTTGGCTTCCGGCGATCCTCtttgaac                                                                                   <  1:160152‑M1/144‑4 (MQ=255)                                                                                                                                                                             aatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcccTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATcg                                                  >  1:152911‑M1/75‑251 (MQ=255)                                                                                                                                                                                cggatcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcccTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACCGCTGCTGGCATTAACCCCGGTTATCgcgc                                               <  2:93003‑M1/180‑1 (MQ=255)                                                                                                                                                                                           tACAAATAAAAGGTGTTCACGTGCAGAATTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCCTTAACGCCGGTTATCGCGACATTAAATGGt                                    >  1:10399/1‑251 (MQ=255)                                                                                                                                                                                                aaatgaaaggactgagtctgccgagaagaatttccccgcttattcgcaccctcccTAGCAATCCCGGTGGTTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGGGATGGCGGGCGGCTATAGGTCCACCGACATGGCGGCGGTCGCTATCGCTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTAGCGCGCAATGAAATGGttcctg                               >  2:158652‑M1/56‑249 (MQ=255)                                                                                                                                                                                                aTAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGATATCGCGCCATTAAATGGTTCCgg                               >  1:28969/1‑251 (MQ=255)                                                                                                                                                                                                                             atttccccgcttattcgcaccttcccTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCTCAAGTAAATGGGTTCGGTCGACGTGAGCGCATTGCGCATCAGGTGc  >  1:287432‑M1/27‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 CTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCTGCCCATTCAGGCGGCTTTTTAGTCTCTCATATAACTACAAATAAAAGGTGTTCACGTGCAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGC  >  NC_000913/1743251‑1743721

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