New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1423890	12 (0.210)	4 (0.070)	4/478	2.9	8.5%	coding (109/1458 nt)	trkG	K(+) transporter TrkG
	?	NC_000913	= 4356677	74 (1.300)					coding (1941/2148 nt)	cadA	lysine decarboxylase 1
Rejected: Frequency below/above cutoff threshold.

AAATCTTAACTACAATATGATTCATTTTGATGATTGTTTCATAAGGAACAGTGAAGTAAGATCTAAGAGGAGTTAAATTTTATACAGTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCTGGAATAACAACAACAGGAGCAACTGTA  >  NC_000913/1423644‑1424138                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         aaaTCTTAACTACAATATGATTCATTTTGATGATTGTTTCATAAGGAACAGTGAAGTAAGATCTAAGAGGAGTTAAATTTTATACAGTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAg                                                                                                                                                                                                                                                      <  2:407366/251‑1 (MQ=255)                                                aCAGTGAAGTAAGAGCTAAGATGAGTTAAATTTTATACAGTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTTGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATtggtg                                                                                                                                                                                                       <  2:258926/251‑1 (MQ=255)                                                          aaGATCTAATAGGAGTTAAATTTTATACAGGATAATCATAATATTGCAGCAAGGTAGGTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCTTCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTAGAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTg                                                                                                                                                                                             <  2:492483/251‑1 (MQ=255)                                                                                                                tggtTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAtatgatccttccgtacccgccgggagttcctctggtaatgccgggtgaaatgatcaccgaagaaagccgtccggttctggagttcctgcagatgctgtgtgaactc                                                                                                                                                >  2:208712‑M1/1‑136 (MQ=255)                                                                                                                tggtTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAtatgatccttccgtacccgccgggagttcctctggtaatgccgggtgaaatgatcaccgaagaaagccgtccggttctggagttcctgcagatgctgtgtgaaatc                                                                                                                                                <  1:208712‑M1/242‑107 (MQ=255)                                                                                                                        aTTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCAAATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAACAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTa                                                                                                                                >  2:367980/1‑250 (MQ=255)                                                                                                                                   atTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGAtt                                                                                                                    >  2:250305/1‑251 (MQ=255)                                                                                                                                             gAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTg                                                                                                          >  1:51814/1‑251 (MQ=255)                                                                                                                                               aTACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGtt                                                                                                        <  1:564137/251‑1 (MQ=255)                                                                                                                                                          gTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCAt                                                                                              <  1:195550/250‑1 (MQ=255)                                                                                                                                                                                                   ctATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTg                                                    >  2:141674/1‑251 (MQ=255)                                                                                                                                                                                                            ttCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCCTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTg                                            <  2:51814/250‑1 (MQ=255)                                                                                                                                                                                                                      ccaccaATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACCACTAAGAAATCTGGCATTCAATTACGTACCCGTGCTGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGt                                 >  2:510984/1‑251 (MQ=255)                                                                                                                                                                                                                                         tGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTTCCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCTGGAATaacaacaac               <  1:510984/250‑1 (MQ=255)                                                                                                                                                                                                                                            tttccgtgcaaatatgatctctccgtacccgccgggagttcctctggtaatgccgggtgaaatgatcaccgaagaaagccggccggttctggagttccttcagatgctgtgtgaaatcggcgctcactatccgggctttgaaaccgatattcacggtgcataccgtcaggctgatggccgctataccgttaaggtattgaaagaagaaagcaaaaaataattagctcgtacaagggaagtggcttgcca             <  2:440964‑M1/240‑238 (MQ=255)                                                                                                                                                                                                                                                     aaatatgatccttccgtacccgccgggagttcctctggtaatgccgggtgaaatgatcaccgaagaaagccgtccggttctggagttcctgcagatgctgtgtgaaatcggcgctcactatccgggctttgaaaccgatattcacggtgcataccgtcaggctgatggccgctataccgttaaggtattgaaagaagaaagcaaaaaataattagctcgtacaagggaagtggcttgccacttcccttttt  <  2:374426‑M1/251‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         AAATCTTAACTACAATATGATTCATTTTGATGATTGTTTCATAAGGAACAGTGAAGTAAGATCTAAGAGGAGTTAAATTTTATACAGTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGGATTTTGTTTTCTGTTATTAGTGCATTCCCTTTATGGATTGACTCAGAACTTAATTTAACGTTCATTGATGCTCTGTTTGAAGGGGTTTCTGGAATAACAACAACAGGAGCAACTGTA  >  NC_000913/1423644‑1424138

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