New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 2003635	19 (0.440)	3 (0.070)	3/242	NT	13.9%	coding (93/1011 nt)	PP_1787	hypothetical protein
	?	NC_002947	= 2003657	19 (0.460)					coding (115/1011 nt)	PP_1787	hypothetical protein

TCTACATGACCTTGAAAGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGAGCTAGCCGGAGCCTTCGAAAAATTTAGGC  >  NC_002947/2003494‑2003783                                                                                                                                              |                                                                                                                                                       tctACATGACCTTGAAAGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGcc                                                                                                                                                >  2:66747/1‑150 (MQ=255)  ctACATGACCTTGAAAGATAGAACTCGAGTATGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCa                                                                                                                                               <  2:985538/150‑1 (MQ=255)    aCATGACCTTGAAAGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCAtt                                                                                                                                             >  1:572242/1‑150 (MQ=255)           cTTGAAAGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATgg                                                                                                                                      >  2:1072594/1‑150 (MQ=255)           cTTGAAAGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGAGAAGTCCTAACTCCTGATTGAATCTCTGTACGACGTTGTTTATGCATTCATGGAGGTTAGCCAGTCGCAGCTTGAGGTGCCATTCAAATgg                                                                                                                                      <  2:680237/150‑1 (MQ=255)                 aGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAac                                                                                                                                >  2:292114‑M2/1‑148 (MQ=255)                 aGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAac                                                                                                                                <  1:654723‑M2/150‑3 (MQ=255)                          tGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTtcacccatttgaatggcacctcacgctgcgcct                                                                                                                       <  2:898996‑M1/150‑34 (MQ=255)                                  gaaggaagCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTcagca                                                                                                             >  2:894071/1‑150 (MQ=255)                                              taAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGa                                                                                                 <  1:1072594/150‑1 (MQ=255)                                               aaaaTGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAg                                                                                                <  2:582278/150‑1 (MQ=255)                                               aaaaTGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTACAAATGTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAg                                                                                                <  1:912587/150‑1 (MQ=255)                                                    gAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGCTC‑‑CGCAGCAGGTAGCCATTGAGGTCaa                                                                                           <  1:426928/150‑1 (MQ=255)                                                           aCAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCag                                                                                     <  1:221426/149‑1 (MQ=255)                                                           aCAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCag                                                                                     >  2:221426/1‑149 (MQ=255)                                                                      gCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGt                                                                         >  1:459707/1‑150 (MQ=255)                                                                             tCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTtcacccatttgaatggcacctcacgctgcgcctg                                                                                                                      <  1:694967‑M1/100‑35 (MQ=255)                                                                             tCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTtcacccatttgaatggcacctcacgctgcgcctg                                                                                                                      >  2:694967‑M1/1‑66 (MQ=255)                                                                                       ttGAATTTGCGAGCGTGAATGCTTCCAAGTGTATATGTTGTAGCCGGGGGGAACGTGAGGTCGCAATCACATGGAGGATGTCTTGTGTTC‑‑CGCGTCAGGGAGCCATTGTGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCtgc                                                        <  2:459707/150‑3 (MQ=255)                                                                                                       cAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCa                                        <  1:755603/150‑1 (MQ=255)                                                                                                          tGCTTCCAAATTTATAAATCTTAGCCAGTCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGa                                     <  1:19977/150‑1 (MQ=255)                                                                                                                cAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCAATTCAAATGGGTGAACGCTGGAGCTC‑‑CTCAGCAGGTGGCCATTGTGGTTGGGATCAGAGAGTTGTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGa                               >  2:971224/1‑150 (MQ=255)                                                                                                                 aaaTTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGAg                              >  1:1036078/1‑150 (MQ=255)                                                                                                                                  ccAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAGCG‑TTGCAGTCAGGACGGCA‑GCAACCAGTGAGGGCAAGAGCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGAGCTAGCCGGAGCCTTCGa             >  1:570670/1‑150 (MQ=255)                                                                                                                                   cAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAacgctgcgcctggctaagatttataaatttggaagcattgacgctcgcaaattcaatcacgactga                                                                >  2:694967‑M2/1‑34 (MQ=255)                                                                                                                                   cAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAacgctgcgcctggctaagatttataaatttggaagcattgacgctcgcaaattcaatcacgactga                                                                <  1:694967‑M2/100‑67 (MQ=255)                                                                                                                                    aGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAacgctgcgcctggctaagatttataaatttggaagcattgacgctcgcaaattcaatcacgactgaggacagccctatatttgtaactttcattttataagcttccttcctactgca             <  2:898996‑M2/150‑118 (MQ=255)                                                                                                                                             gTtcacccatttgaatggcacctcacgctgcgcctggctaagatttataaatttggaagcattgacgctcgcaaattcaatcacgactgaggacagccctatatttgtaactttcattttataagcttccttcctactgcatcgccatct    >  2:292114‑M1/1‑2 (MQ=255)                                                                                                                                             gTtcacccatttgaatggcacctcacgctgcgcctggctaagatttataaatttggaagcattgacgctcgcaaattcaatcacgactgaggacagccctatatttgtaactttcattttataagcttccttcctactgcatcgccatct    <  1:654723‑M1/150‑149 (MQ=255)                                                                                                                                             gTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGAGCTAGCCGGAGCCTTCGAAAAATTTAGGc  <  1:66747/150‑1 (MQ=255)                                                                                                                                             gTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGAGCTAGCCGGAGCCTTCGAAAAATTTAGGc  <  2:1036078/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                       TCTACATGACCTTGAAAGATGGCGATGCAGTAGGAAGGAAGCTTATAAAATGAAAGTTACAAATATAGGGCTGTCCTCAGTCGTGATTGAATTTGCGAGCGTCAATGCTTCCAAATTTATAAATCTTAGCCAGGCGCAGCGTGAGGTGCCATTCAAATGGGTGAACGCTGGCGATC‑‑CTCAGCAGGTAGCCATTGAGGTCAATATCAGAGATTTTTCGGTCTATGAAAGTTTGCTACTGACTTCTGACGAGCATGAGTTGGAGCTAGCCGGAGCCTTCGAAAAATTTAGGC  >  NC_002947/2003494‑2003783

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