New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3839257 =	27 (0.620)	2 (0.050)	4/240	NT	8.1%	coding (137/1158 nt)	PP_3389	membrane protein
	?	NC_002947	3839280 =	20 (0.490)					coding (160/1158 nt)	PP_3389	membrane protein

GAATGGGGATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGC  >  NC_002947/3839133‑3839425                                                                                                                                                     |                                                                                                                                                  gAATGGGGATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAaca                                                                                                                                                 >  2:221208/1‑150 (MQ=255) gAATGGGGATTCTGGCAGATGGCAACATATTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAaca                                                                                                                                                 >  1:870629/1‑150 (MQ=255)  aaTGGGGATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAacac                                                                                                                                                <  1:668472/150‑1 (MQ=255)  aaTGGGGATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAacac                                                                                                                                                <  2:440304/150‑1 (MQ=255)   aTGGGGATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGAACGCAAACAGAGATGTTTATGAGGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACAcc                                                                                                                                               >  1:569242/1‑150 (MQ=255)     ggggATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACAcccc                                                                                                                                             <  2:1103039/150‑1 (MQ=255)     ggggATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACAcccc                                                                                                                                             <  1:221208/150‑1 (MQ=255)                     ggCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAACCACCCCTCCTACGCGGCCCTGt                                                                                                                             <  2:1093204/150‑1 (MQ=255)                         aCTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTAt                                                                                                                         >  1:489566/1‑150 (MQ=255)                             agtcagAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCAt                                                                                                                     <  1:447974/150‑1 (MQ=255)                                          aCGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCca                                                                                                        >  2:319156/1‑150 (MQ=255)                                            gTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCcaca                                                                                                      >  2:207793/1‑150 (MQ=255)                                             tAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACAt                                                                                                     <  1:450931/150‑1 (MQ=255)                                                 ttctcttctgagccgccacagcgcatttccatgtgggtttcattgaatgcataaacagggccgcgtatgaggggtgtttcacttttcccgggtaagtgaaACAcc                                                                                                                                               >  1:121884‑M2/101‑105 (MQ=255)                                                 ttctcttctgagccgccacagcgcatttccatgtgggtttcattgaatgcataaacagggccgcgtatgaggggtgtttcacttttcccgggtaagtgaaACAcc                                                                                                                                               <  2:121884‑M2/5‑1 (MQ=255)                                                   gTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTc                                                                                                                                           <  1:238792/106‑1 (MQ=255)                                                   gTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTc                                                                                                                                           >  2:238792/1‑106 (MQ=255)                                                    tACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAAt                                                                                               >  2:479567/1‑149 (MQ=255)                                                    tACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAAt                                                                                               <  1:479567/149‑1 (MQ=255)                                                        aaGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCt                                                                                          >  1:1039132/1‑150 (MQ=255)                                                             gggATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACACCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGACATGCGCTGTggc                                                                                     >  2:204492/1‑150 (MQ=255)                                                                         gcagAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCagaagag                                                                         <  1:834369/150‑1 (MQ=255)                                                                                            aagaTTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGc                                                      <  2:1039132/149‑1 (MQ=255)                                                                                              aatgcataaacagggccgcgtatgaggggtgtttcacttacccgggaaaagtgaaACACCCCTCATAcg                                                                                                                                      <  1:678632‑M2/14‑1 (MQ=255)                                                                                              aatgcataaacagggccgcgtatgaggggtgtttcacttacccgggaaaagtgaaACACCCCTCATAcg                                                                                                                                      >  2:678632‑M2/56‑69 (MQ=255)                                                                                               aTTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAg                                                   >  2:1120283/1‑150 (MQ=255)                                                                                                               cgtatgaggggtgtTTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCAtt                                                                                                                    <  1:678632‑M1/55‑1 (MQ=255)                                                                                                               cgtatgaggggtgtTTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCAtt                                                                                                                    >  2:678632‑M1/15‑69 (MQ=255)                                                                                                                    atgaaGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGtgctgc                              >  2:939902/1‑150 (MQ=255)                                                                                                                        ggtgtTTCACTTA‑CCCGGGAAAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCagaagagaa                                                                       <  2:121884‑M1/100‑1 (MQ=255)                                                                                                                        ggtgtTTCACTTA‑CCCGGGAAAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCagaagagaa                                                                       >  1:121884‑M1/6‑105 (MQ=255)                                                                                                                            aTTCACTTTTCCCGGGTAA‑TTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTc                                                             >  2:1029767/1‑111 (MQ=255)                                                                                                                            aTTCACTTTTCCCGGGTAA‑TTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTc                                                             <  1:1029767/111‑1 (MQ=255)                                                                                                                             TTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATaa                     >  1:939902‑M1/1‑150 (MQ=255)                                                                                                                             TTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATaa                     >  1:207793‑M1/1‑150 (MQ=255)                                                                                                                                  TTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTc                <  2:710519‑M1/150‑1 (MQ=255)                                                                                                                                        cGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGc                                                                >  2:134502/1‑96 (MQ=255)                                                                                                                                        cGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAAcc          >  2:311974/1‑150 (MQ=255)                                                                                                                                        cGGGTAA‑GTGAAACACACCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGc                                                                <  1:134502/96‑1 (MQ=255)                                                                                                                                               agtcgaACACCCCTCCTACGCGGCCCTGTTTATGCATTCAATGAATCCCGCGTGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCTGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATa    <  2:401534‑M2/144‑1 (MQ=255)                                                                                                                                                gtgaaACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCa                                                                               >  2:104851‑M2/6‑74 (MQ=255)                                                                                                                                                gtgaaACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCa                                                                               <  1:104851‑M2/69‑1 (MQ=255)                                                                                                                                                 tgaaACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGc  >  2:1070055‑M2/5‑150 (MQ=255)                                                                                                                                                 tgaaACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGc  >  1:880878‑M2/5‑150 (MQ=255)                                                                                                                                                 tgaaACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAAGCGGCTGTGGCGGCTTAGAAGACAAGCGCGGGCCCATTACGATAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGc  >  1:522981‑M2/5‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GAATGGGGATTCTGGCAGATGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAA‑GTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGC  >  NC_002947/3839133‑3839425

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.