New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5895316	24 (1.480)	5 (0.330)	5/242	NT	17.3%	coding (86/873 nt)	cysW	sulfate ABC transporter permease CysW
	?	NC_002947	= 5895363	25 (1.630)					coding (39/873 nt)	cysW	sulfate ABC transporter permease CysW

ACCACGGCGAGCAGGGTCAGCTTCAGCGCCGACAACGCATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATG  >  NC_002947/5895169‑5895464                                                                                                                                                    |                                                                                                                                                     aCCACGGCGAGCAGGGTCAGCTTCAGCGCCGACAACGCATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAgg                                                                                                                                                    >  1:80708/1‑150 (MQ=255)          aGCAGGGTCAGCTTCAGCGCCGACAACGCATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAgg                                                                                                                                           >  2:44669/1‑150 (MQ=255)                           cgcCGACAACGCATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGgc                                                                                                                          <  2:166092/150‑1 (MQ=255)                                  aaCGCATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCa                                                                                                                   >  2:30673/1‑150 (MQ=255)                                      cATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGTCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCAcggc                                                                                                               <  1:26498/150‑1 (MQ=255)                                             ttCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGggcggc                                                                                                        <  1:44669/150‑1 (MQ=255)                                                gAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGtt                                                                                                     <  2:309684‑M2/150‑3 (MQ=255)                                                gAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGtt                                                                                                     <  2:265913‑M2/150‑3 (MQ=255)                                                     tCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctg                                                                                                >  2:80708‑M2/1‑143 (MQ=255)                                                     tCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGGGGCGCTGCCACGGCGGGCGGCGttggctg                                                                                                <  2:49678‑M2/150‑8 (MQ=255)                                                        ccTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCgg                                                                                             >  1:168898/1‑150 (MQ=255)                                                         cTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAcgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggccttggctggctg                                                                                            <  1:412572‑M1/150‑59 (MQ=255)                                                              agaagGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAcgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggcctt                                                                                                     >  2:381245‑M1/1‑87 (MQ=255)                                                              agaagGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAcgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggcctt                                                                                                     <  1:381245‑M1/136‑50 (MQ=255)                                                                   gTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGc                                                                                               >  2:345752/1‑137 (MQ=255)                                                                   gTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGc                                                                                               <  1:345752/137‑1 (MQ=255)                                                                    tGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCa                                                                                 >  1:268616/1‑150 (MQ=255)                                                                     gCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAg                                                                                <  1:170575/150‑1 (MQ=255)                                                                        aaaaCCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAcgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggccttggctggctgg                                                                                           >  1:210964‑M1/1‑77 (MQ=255)                                                                        aaaaCCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAcgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggccttggctggctgg                                                                                           <  2:210964‑M1/136‑60 (MQ=255)                                                                             cGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGAt                                                                            <  2:415904/146‑1 (MQ=255)                                                                             cGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGAt                                                                            >  1:415904/1‑146 (MQ=255)                                                                             cGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAAc                                                                        >  1:342543/1‑150 (MQ=255)                                                                                 cttcAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGa                                                                    <  2:342543/150‑1 (MQ=255)                                                                                 cttcAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGa                                                                    <  2:268616/150‑1 (MQ=255)                                                                                           tGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGcc                                                          >  1:68368/1‑150 (MQ=255)                                                                                               aTACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGcctcct                                                      <  1:30673/150‑1 (MQ=255)                                                                                                  ccACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGTGATGAACTGGACATGCCGGCCTCCTTCa                                                   >  2:112910/1‑150 (MQ=255)                                                                                                   cACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCaa                                                  <  1:265913/150‑1 (MQ=255)                                                                                                   cACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCaa                                                  <  1:266549/150‑1 (MQ=255)                                                                                                    aCGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAg                                                 >  1:306902/1‑150 (MQ=255)                                                                                                                               agcgCGAACACCAGTCAGCCAAGGCCGATCAGGATACGCCGTTACGTGTCGCTGCCACGGCGGGAGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCGCCGGGTTTCGATGCGGCGCTGCAGCAGGt                      <  1:216617/147‑1 (MQ=255)                                                                                                                                  gCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGa                   >  1:40317/1‑150 (MQ=255)                                                                                                                                         ccagccagCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggtttt                          <  2:210964‑M2/136‑78 (MQ=255)                                                                                                                                         ccagccagCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggtttt                          >  1:210964‑M2/1‑59 (MQ=255)                                                                                                                                          cagccagcCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttcttcgag           <  1:412572‑M2/150‑93 (MQ=255)                                                                                                                                          cagccagcCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATcagcagc           >  2:144391/1‑150 (MQ=255)                                                                                                                                              cagccaAcgccgcccgccgtggcagcgccccgtcacggcgtatcctgatcggccttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttcttcgaggcga       <  2:49678‑M1/150‑144 (MQ=255)                                                                                                                                              cagccaAcgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggccttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttcttcgaggcga       >  2:80708‑M1/1‑7 (MQ=255)                                                                                                                                                  caAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGAt   >  1:154213/1‑150 (MQ=255)                                                                                                                                                   aacgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggccttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttcttcgaggcgatattc  <  2:309684‑M1/150‑149 (MQ=255)                                                                                                                                                   aacgccgcccgccgtggcagcgccacgtcacggcgtatcctgatcggccttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttcttcgaggcgatattc  <  2:265913‑M1/150‑149 (MQ=255)                                                                                                                                                   aaGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttct                <  1:381245‑M2/136‑88 (MQ=255)                                                                                                                                                   aaGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttct                >  2:381245‑M2/1‑49 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACCACGGCGAGCAGGGTCAGCTTCAGCGCCGACAACGCATCGGGTTCGAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATG  >  NC_002947/5895169‑5895464

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