New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	3437896 =	18 (0.600)	4 (0.150)	4/262	3.6	20.8%	intergenic (+2/‑128)	trkA/mscL	NAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
	?	NC_000913_3_bme_pgi	3437924 =	14 (0.510)					intergenic (+30/‑100)	trkA/mscL	NAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Rejected: Coverage evenness skew score above cutoff.

GTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGA  >  NC_000913_3_bme_pgi/3437761‑3438040                                                                                                                                            |                                                                                                                                                 gTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTCATTCACGTT‑G‑TCA‑‑CAGATAGTATGGGTATTACCGACGTCGAAAGAATCTTCCAGCCAAGCCCTTGCTTCTTGTAATTAATAAGGCGTCTa                                                                                                                                     >  1:518757/1‑149 (MQ=255)       acgacgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACg                                                                                                                               <  2:231418/149‑1 (MQ=255)         gaagTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGcc                                                                                                                             >  2:12591/4‑149 (MQ=255)            gggatgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtta                                                                                                                          <  2:183216/147‑1 (MQ=255)                        aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTaa                                                                                                                                    <  1:294788/127‑1 (MQ=255)                        aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTaa                                                                                                                                    >  2:294788/1‑127 (MQ=255)                        aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACAGACGTCGAAAGACTCTTCAAGCAAAGCCCTTTCTTCTTGTAATTAATATGGAGTCGAATGACGCCTTATTGTTTACCttt                                                                                                              >  2:415214/1‑149 (MQ=255)                                  tGCGCATTGAGCAAGGCGATCCAGTAATG‑AGGTTCCTTCATCCAGCTATACAG‑‑‑AGTATCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                    <  1:514554/149‑1 (MQ=255)                                                   gATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattcatt                                                                                   <  1:162336/149‑1 (MQ=255)                                                           aaTT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTATTTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAg                                                                           >  2:480969/1‑149 (MQ=255)                                                                  gTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAAt                                                                     <  1:451762/149‑1 (MQ=255)                                                                              aGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACa                                                            >  1:359127/1‑149 (MQ=255)                                                                                   aaaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGt                                                       >  2:241661/1‑149 (MQ=255)                                                                                     aaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                     <  2:359127/149‑1 (MQ=255)                                                                                        tttATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGAc                                                  <  2:127417/149‑1 (MQ=255)                                                                                                       gAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTGGTTAGACGTATGGTTGTAGGCt                                   >  1:397855/1‑149 (MQ=255)                                                                                                                ttgcAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGga                          <  2:100931/146‑1 (MQ=255)                                                                                                                                        taattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGa  >  2:384190/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                 GTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATT‑ATGTTCC‑TCA‑‑CAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGA  >  NC_000913_3_bme_pgi/3437761‑3438040

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