New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4498181	1 (0.030)	9 (0.320)	8/476	1.5	72.1%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	4499518 =	6 (0.210)					intergenic (+18/‑1742)	insD6/yjgZ	IS2 insertion element protein InsB/uncharacterized protein YjgZ

AAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTG  >  NC_000913/4497932‑4498392                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                    aaaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGAAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTACGGAGTGCTGCAGGGATTTGATGGTCCCTTGTg                                                                                                                                                                                                                    >  1:352241‑M1/1‑250 (MQ=255)                                                              aaaGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagt                                                                                                                                                       <  1:183711‑M1/251‑63 (MQ=255)                                                                       cGGTATTAACTGATGTATTTACTATTTACCGCTCACTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGAGTGATGCGGGATAACTGACGTAAGGAGTGCTGCTGGGATTTGATGGTCCCTTGTgctgttatgaatacgaatattttatacagcatacttttcagtacaagcgactgtattacagtctgccatat                                                                                                                                              >  1:148‑M1/1‑180 (MQ=255)                                                                                                    cGCTCCCTGCAGGATGGATTATCAGAATCAGCGGTTATTCATTATTTATATCCATTCAATAAAAAAGTAACCTCATAAAAAGTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCATTGTgctgatatgaatacggatattttatacagcatacttttcagaacaagcgactgtattacagtctgccatcaatgattaccggtgatcaggaagtataagg                                                                                                                 >  2:281019‑M1/1‑151 (MQ=255)                                                                                                        cccttcagtggtgattaccaatatccgcggttattcattatttatACACATTCAATAAAAAAGTAACCTCAAAAAAAATCCTGTACTCCCAGGGTTATGCGTGATAACTGACTTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgatttcctgtgatcaggaagaattagtgg                                                                                                               <  2:269781‑M1/227‑103 (MQ=255)                                                                                                             cAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggt                                                                                                        <  1:347078‑M1/251‑110 (MQ=255)                                                                                                                                                    tataCATTCAATAAAGAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGG‑GATAAGATGATTGATGTCTTAGGGCCGGTGAAACGCAGACGGCGTACCACACAGGAAAAGATCGc                                                                >  2:195296/1‑251 (MQ=255)                                                                                                                                                     ataCAGTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtcttccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaa                                                                <  1:183119‑M1/251‑150 (MQ=255)                                                                                                                                                                                             cTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataaccggatagttatcttcgtctgtccgcggattga                         >  2:183872‑M1/1‑62 (MQ=255)                                                                                                                                                                                                                     gACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaa  <  1:353643‑M1/249‑212 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                    AAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTG  >  NC_000913/4497932‑4498392

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