New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 2137517	32 (1.080)	4 (0.140)	4/234	NT	11.0%	coding (334/624 nt)	gstB	glutathione S‑transferase
	?	NC_002947	= 2137523	34 (1.210)					coding (340/624 nt)	gstB	glutathione S‑transferase

GGAATCCAACACCATCGTCCGCTACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGATCGGCATGGGCGACA  >  NC_002947/2137375‑2137670                                                                                                                                                     |                                                                                                                                                    ggAATCCAACACCATCGTCCGCTACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTg                                                                                                                                                    >  2:332968‑M2/1‑149 (MQ=255) ggAATCCAACACCATCGTCCGCTACCTGTGCGCGGAGTACGGCACGGAGCAGGGCCGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTg                                                                                                                                                    <  2:381199‑M2/150‑2 (MQ=255)  gAATCCAACACCATCGTCCGCTACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTgg                                                                                                                                                   >  1:50159‑M2/1‑148 (MQ=255)  gAATCCAACACCATCGTCCGCTACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTgg                                                                                                                                                   <  2:667236‑M2/150‑3 (MQ=255)                       tACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGaa                                                                                                                              <  1:679012/150‑1 (MQ=255)                       tACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGaa                                                                                                                              <  2:775039/150‑1 (MQ=255)                         ccTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCgg                                                                                                                                >  2:500400/1‑146 (MQ=255)                         ccTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCgg                                                                                                                                <  1:500400/146‑1 (MQ=255)                          cTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGAc                                                                                                                           <  1:577412/150‑1 (MQ=255)                                      tACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTgggg                                                                                                                                                 >  2:281484‑M2/1‑112 (MQ=255)                                      tACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTgggg                                                                                                                                                 <  1:281484‑M2/116‑5 (MQ=255)                                          gCACGGAGCAGCGCTGGTATCTGGCAGCTGAACGGCAGCGCTCCCTCGGTGACAAATGCATTGGCGTGAGTACATCTTCCGGCGCTACTCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTgg                                                                                                           >  2:209869/1‑150 (MQ=255)                                                     ggCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAAcgcc                                                                                                >  2:238677/1‑150 (MQ=255)                                                          gTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGca                                                                                           <  2:723090/150‑1 (MQ=255)                                                                     aCCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCAc                                                                                                                                     >  2:165278/1‑97 (MQ=255)                                                                     aCCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCAc                                                                                                                                     <  1:165278/97‑1 (MQ=255)                                                                       ccACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTgcgcg                                                                              >  2:80630/1‑150 (MQ=255)                                                                       ccACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCAGAACAggggccggaatggg                                                                                                                                       >  2:453544‑M2/1‑79 (MQ=255)                                                                       ccACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCAGAACAggggccggaatggg                                                                                                                                       <  1:453544‑M2/93‑15 (MQ=255)                                                                       ccACGCCAGCGCGCCCTGGCCGACAAATGGAAGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGa                                                                                                        <  2:288679/124‑1 (MQ=255)                                                                       ccACGCCAGCGCGCCCTGGCCGACAAATGGAAGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGa                                                                                                        >  1:288679/1‑124 (MQ=255)                                                                          cGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCaa                                                                           <  2:22827/150‑1 (MQ=255)                                                                                       tGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGTTAGGGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTg                                                              <  1:209869/150‑1 (MQ=255)                                                                                       tGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACt                                                                      >  1:674967/1‑142 (MQ=255)                                                                                       tGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACt                                                                      <  2:674967/142‑1 (MQ=255)                                                                                                  tggatggaCTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGTGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAAcc                                                   <  1:28825/150‑1 (MQ=255)                                                                                                  tggatggaCTCGACCACATAATCCTTCGCTGCACCATTACGGTAACAGTTGGTGGGTCTTCTGCCCACGGCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAAcc                                                   <  2:372803/150‑1 (MQ=255)                                                                                                   ggatggaCTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAAccc                                                  >  2:22867/1‑150 (MQ=255)                                                                                                           tggacCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTgg                                                                                                               >  2:205113/1‑81 (MQ=255)                                                                                                           tggacCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTgg                                                                                                               <  1:205113/81‑1 (MQ=255)                                                                                                                    tCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGTCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTTGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTAc                                 <  2:266781/150‑1 (MQ=255)                                                                                                                      gTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCt                               <  2:235365/150‑1 (MQ=255)                                                                                                                        ccTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTa                                  <  1:579835/145‑1 (MQ=255)                                                                                                                        ccTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTa                                  >  2:579835/1‑145 (MQ=255)                                                                                                                          ttCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTAc                                 >  1:688821/1‑144 (MQ=255)                                                                                                                          ttCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTAc                                 <  2:688821/144‑1 (MQ=255)                                                                                                                                aCCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCGCGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTgg                                                                    >  2:351036/1‑103 (MQ=255)                                                                                                                                aCCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCGCGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTgg                                                                    <  1:351036/103‑1 (MQ=255)                                                                                                                                     aTTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGAt                <  1:667236/150‑1 (MQ=255)                                                                                                                                       tCCGCCCCCTGTTCTGGGGCCTTCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACCAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGATCg              <  1:381199/150‑1 (MQ=255)                                                                                                                                         cGGCCCCTGTTCTggggccggaatggggtagcgaaggacgatgtggtccagtcc                                                                                                            >  2:302190‑M2/1‑13 (MQ=255)                                                                                                                                         cGGCCCCTGTTCTggggccggaatggggtagcgaaggacgatgtggtccagtcc                                                                                                            <  1:302190‑M2/54‑42 (MQ=255)                                                                                                                                           gCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGATCGGCAt          >  1:70540/1‑150 (MQ=255)                                                                                                                                                tGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGATCGGCATGGGCg     >  1:572943/1‑150 (MQ=255)                                                                                                                                                   tCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGATCGGCATGGGCGACa  <  1:22867/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGAATCCAACACCATCGTCCGCTACCTGTGCGCGGAGTACGGCACGGAGCAGGGCTGGTATCTGGCAGACCCACGCCAGCGCGCCCTGGCCGACAAATGGATGGACTGGACCACATCGTCCTTCGCTACCCCATTCCGGCCCCTGTTCTGGGGCCTGCTGCGCACGCCGGAAGACCAGCGCGACTGGGTGGCGATCAACGCCGCGCACAAGCAGTGCGCGCAATTACTGGCCATTGCCGACGAAACCCTGGCCAGTCAGCCGTACCTGTCCGGTGACCAGATCGGCATGGGCGACA  >  NC_002947/2137375‑2137670

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