New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4648496	28 (0.570)	6 (0.130)	4/242	NT	16.8%	coding (643/1461 nt)	PP_4112	sulfate transporter
	?	NC_002947	= 4648526	33 (0.710)					coding (613/1461 nt)	PP_4112	sulfate transporter

CACACACTCGCGGTCCGGGAAGCCACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAG  >  NC_002947/4648350‑4648640                                                                                                                                                    |                                                                                                                                                   cacacaCTCTCGGTCCGGGAACCCACGGAAATCGGTGATTTCGAACGTCACGATGAGGGTGAGCAGGGTGTCCAGCAGGCCGACCATGGCCAAAAGCACCGCGTAGGGCGCGATGATCTTCAGGGTTTCCAGATTCCAAGACACAT‑CCg                                                                                                                                                  >  4:276725/1‑149 (MQ=255)  acacacTCGCGGTCCGGGAAGCCACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCgg                                                                                                                                                 >  2:244999/1‑149 (MQ=255)  acacacTCGCGGTCCGGGAAGCCACGGCTCTCGGTGATTTAGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCgg                                                                                                                                                 >  7:51973/1‑149 (MQ=255)                                      aTTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTgcgc                                                                                                             <  1:244999/149‑1 (MQ=255)                                       tttCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑Cgccggtggcctgccggggctggcgctgccggatgtgcctt                                                                                                            <  5:280648‑M1/149‑41 (MQ=255)                                       tttCGTCGGGCAGGTTGAGGGTGATCAGTGTTTCCATCAGGCCGACCATGTCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑Cgccggtggcctgccggggctggcgctgccggatgtgcctt                                                                                                            <  6:105525‑M1/149‑41 (MQ=255)                                                         gggTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑Cgccggtggcctgccggggctggcgctgccggatgtgccttggaatctg                                                                                                    <  6:235668‑M1/139‑49 (MQ=255)                                                         gggTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑Cgccggtggcctgccggggctggcgctgccggatgtgccttggaatctg                                                                                                    >  5:235668‑M1/1‑91 (MQ=255)                                                          ggTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑ggg                                                                                       <  2:181426/149‑1 (MQ=255)                                                          ggTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAA‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑ggg                                                                                       >  4:17715/1‑149 (MQ=255)                                                          ggTGAGCAGGGTGTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑ggg                                                                                       <  4:122724/149‑1 (MQ=255)                                                             gAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAg                                                                                    <  7:210021/149‑1 (MQ=255)                                                                cAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAgg                                                                                   >  4:171373/1‑147 (MQ=255)                                                                  gggTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCTGCCGGGGCTGGCG‑CTGCCggatgtgccttggaatctg                                                                                                   >  8:15121‑M2/1‑111 (MQ=255)                                                                  gggTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCTGCCGGGGCTGGCG‑CTGCCggatgtgccttggaatctg                                                                                                   <  7:15121‑M2/130‑20 (MQ=255)                                                                         cagcagGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCa                                                                        >  7:160895/1‑149 (MQ=255)                                                                          agcagGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCAc                                                                       <  8:160895/149‑1 (MQ=255)                                                                               gCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGAGCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTaca                                                                  >  6:255783/1‑149 (MQ=255)                                                                                 cGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAgg                                                                                   <  3:106377/130‑1 (MQ=255)                                                                                 cGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAgg                                                                                   >  4:106377/1‑130 (MQ=255)                                                                                      aTGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGt                                                                                      >  6:316399/1‑122 (MQ=255)                                                                                      aTGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGt                                                                                      <  5:316399/122‑1 (MQ=255)                                                                                                 cACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAacac                                                   <  4:238414/146‑1 (MQ=255)                                                                                                 cACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAacac                                                   >  3:238414/1‑146 (MQ=255)                                                                                                      cgTAGGGGGGGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGGCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATc                                           <  8:51973/149‑1 (MQ=255)                                                                                                               cGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGAACGTACACCACCAAGCCAACACCGAGGATCGccacacag                                  >  5:86398/1‑145 (MQ=255)                                                                                                                      cTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCcggcg                           >  7:186687/1‑149 (MQ=255)                                                                                                                         cAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCcggc                            <  3:17715/145‑1 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACATCCggcagcgccagccccggcaggcgctgccggatgtgccttggaatctggaaaccctgaagatcatcgccccctacgcggtgctgatggccatggtcggcctgctggaaaccc                                     <  7:15121‑M1/130‑112 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACATCCggcagcgccagccccggcaggcgctgccggatgtgccttggaatctggaaaccctgaagatcatcgccccctacgcggtgctgatggccatggtcggcctgctggaaaccc                                     >  8:15121‑M1/1‑19 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCgatgtgccttggaatctggaaaccctgaagatcatcgccccctacgcggtgctgatggccatggtcggcctgctggaaaccctgctcaccc                           <  6:235668‑M2/139‑92 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCgatgtgccttggaatctggaaaccctgaagatcatcgccccctacgcggtgctgatggccatggtcggcctgctggaaaccctgctcaccc                           >  5:235668‑M2/1‑48 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCTTGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACcgcg                   <  6:9338/145‑1 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACcgcg                   >  5:9338/1‑145 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCgg                >  1:289775/1‑148 (MQ=255)                                                                                                                                  cAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCgg                <  2:289775/148‑1 (MQ=255)                                                                                                                                   aGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGt              >  7:48601/1‑149 (MQ=255)                                                                                                                                    gATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGt              >  7:162996/1‑148 (MQ=255)                                                                                                                                     aTTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCa            <  5:255783/149‑1 (MQ=255)                                                                                                                                          aaGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCgatgtgccttggaatctggaaaccctgaagatcatcgccccctacgcggtgctgatggccatggtcggcctgctggaaaccctgctcaccctcaacctgaccgacgaaa         <  5:280648‑M2/149‑110 (MQ=255)                                                                                                                                          aaGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCgatgtgccttggaatctggaaaccctgaagatcatcgccccctacgcggtgctgatggacatggtcggcctgatggaaacactgatcaccctcaacctgcccgacgaaa         <  6:105525‑M2/149‑110 (MQ=255)                                                                                                                                          aaGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGTCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTc       <  3:276725/149‑1 (MQ=255)                                                                                                                                          aaGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTc       <  8:162996/149‑1 (MQ=255)                                                                                                                                               acaT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGgcag  >  3:252028/1‑149 (MQ=255)                                                                                                                                                 aT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGAGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGgcag  >  1:249149/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CACACACTCGCGGTCCGGGAAGCCACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACAT‑CCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGC‑‑GGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAG  >  NC_002947/4648350‑4648640

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