New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 781688	0 (0.000)	16 (0.370)	16/380	1.0	100%	intergenic (+36/‑397)	lysQ/nadA	tRNA‑Lys/quinolinate synthase
	?	NC_000913	781845 =	0 (0.000)					intergenic (+193/‑240)	lysQ/nadA	tRNA‑Lys/quinolinate synthase

GTGGTTACTGGTAGAGAACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATT  >  NC_000913/781453‑781860                                                                                                                                                                                                                                            |                                                                                                                                                                             gtggtTACTGGTAGAGAACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaat                                                                                                                                                               >  1:129031‑M1/1‑236 (MQ=255)          ggTAGAGAACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttg                                                                                                                                                      >  1:219949‑M1/1‑227 (MQ=255)              gagaACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaa                                                                                                                                                  <  2:92291‑M1/251‑29 (MQ=255)                                        aGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactg                                                                                                                        <  2:129031‑M1/251‑55 (MQ=255)                                                    aaGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttccctttgtaaacctgatgacatcgtaagagcttactgtgaaagcaactc                                                                                                            >  1:550908‑M1/1‑185 (MQ=255)                                                    aaGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactc                                                                                                            <  1:248021‑M1/251‑67 (MQ=255)                                                       ggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctat                                                                                                         >  1:65145‑M1/1‑182 (MQ=255)                                                                                      ccaccaTCCTGAATGATTAAGGCAGCAGAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatc                                                                          <  1:541117‑M1/251‑101 (MQ=255)                                                                                            tCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcac                                                                    >  2:147642‑M1/1‑145 (MQ=255)                                                                                                        aaGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaat                                                        >  2:355297‑M1/1‑133 (MQ=255)                                                                                                          ggCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattc                                                      <  2:269726‑M1/251‑121 (MQ=255)                                                                                                               caTAATCCCGCAAGGAGTCGATAGCTCAGCTGGTAGAGCAGTTGACTTTTAATCAACTGCTCGAAGGTTCGAATCCTGCACGACCCACCCATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtcaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtg                                                 >  2:273211‑M1/1‑126 (MQ=255)                                                                                                                                 cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttccccaccattcg                               <  1:291984‑M1/251‑144 (MQ=255)                                                                                                                                      gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttccccaccattctactat                          <  2:288213‑M1/251‑149 (MQ=255)                                                                                                                                                  gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcacttagtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttccccaccattcgactatcttgtttagcat              <  2:143714‑M1/251‑161 (MQ=255)                                                                                                                                                               ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttccccaccattcgactatcttgtttagcatataaaacaaatt  <  1:395211‑M1/250‑173 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                             GTGGTTACTGGTAGAGAACGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATT  >  NC_000913/781453‑781860

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