New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1978503 =	NA (NA)	18 (1.810)	17/320	0.0	90.0%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 2406787	2 (0.200)					intergenic (‑146/‑774)	lrhA/alaA	DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA

GTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCC  >  NC_000913/2406589‑2407027                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatcgtgttttatgttcagatcatgcccgatg                                                                                                                                                                                              <  1:7400‑M2/251‑53 (MQ=255)                                            cACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATTAGTTGCACTAAAAAAGCgttaatgactccaacttattgatagtgttttatgtttagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgact                                                                                                                                                   >  1:101647‑M2/1‑156 (MQ=255)                                                                 tACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACGCACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGGAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccgccgattttgagaacgacagcgacttccgtcccagccgtgccag                                                                                                                                >  1:143552‑M2/1‑135 (MQ=255)                                                                     ggcggcTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGGAATTCCATTTGTAATGGGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacgtattgatagtgtgtgatgtttagataattcccgatgactgtgtcatgcagctccaccgattttaagaacgacagcgacttccgtcccagccgtgccaggtgc                                                                                                                            >  1:194798‑M2/1‑131 (MQ=255)                                                                     ggcggcTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagctacttccgtcccagn                                                                                                                                        >  1:174328‑M2/1‑131 (MQ=255)                                                                              tACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCCATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGAAGAACGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattg                                                                                                                                                                                                                             <  1:198699‑M2/143‑22 (MQ=255)                                                                                              gCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCCGTTTCATCCGGCACCCGGACGAACGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctc                                                                                                  <  1:93457‑M2/250‑145 (MQ=255)                                                                                                                         aTTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCgttaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtc                                                                                                                                                                                      >  1:109885‑M2/1‑79 (MQ=255)                                                                                                                                    gCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgcttaattcgctgcgtatatcgcttgctgattacgtgcagct                                                            >  1:86736‑M2/1‑68 (MQ=255)                                                                                                                                     ccTCCCTGGCAGTTTCATCCGGAATCCGGGCGAAAGTACAAATGCATATGAGATGCACTAAAAAAGCggtaatgactacaacttattgatagtgttttatgttcagataatgcccgaagactttgtcatgaagctccaccaattttgagaacgacagcgacttcagtcccagccgtgccaggtgctgcctcagattcaggttatgcc                                                                                                      <  1:209707‑M2/207‑141 (MQ=255)                                                                                                                                             gCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAACAGAggtagtgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttc                                                    <  1:36320‑M2/249‑191 (MQ=255)                                                                                                                                              cAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTCAACAAGCggtattgactcccacttattgatagtgttttttgatcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccacccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgat                                                                       <  1:101681‑M2/229‑172 (MQ=255)                                                                                                                                                                ggACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgncttnctcangcngctcccccgattttgngaacgacagcgacttccgtcccagccgtgccaggtgctgcctcanntnnngntgatgccg                                                                                                     >  2:154888‑M2/1‑40 (MQ=255)                                                                                                                                                                ggACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccg                                                                                                     <  1:154888‑M2/181‑142 (MQ=255)                                                                                                                                                                    gAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttc                                                    <  1:175691‑M2/226‑191 (MQ=255)                                                                                                                                                                                           tcaccaaaaaaGCGCCTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTccc                                                                                                      >  1:10810/6‑152 (MQ=255)                                                                                                                                                                                           gCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccg                                                                                                                                      <  1:210643‑M2/121‑109 (MQ=255)                                                                                                                                                                                             accaaaaaaGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTGCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTAccn  <  1:131054/248‑2 (MQ=255)                                                                                                                                                                                              cTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagcgccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgctt                                                                             >  1:85900‑M2/1‑10 (MQ=255)                                                                                                                                                                                              cTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataaggcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcggatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatc                       >  1:158073‑M2/1‑10 (MQ=255)                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 GTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCC  >  NC_000913/2406589‑2407027

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