New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 589531624 (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

GAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAACACGCCGATGGCGGTCGCGCCGG  >  NC_002947/5895216‑5895507
                                                                                                                                                   |                                                                                                                                                
gAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGtt                                                                                                                                                <  2:265913‑M2/150‑3 (MQ=255)
gAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGtt                                                                                                                                                <  2:309684‑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)
                   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)
                             cGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGAt                                                                                                                       <  2:415904/146‑1 (MQ=255)
                             cGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGAt                                                                                                                       >  1:415904/1‑146 (MQ=255)
                             cGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAAc                                                                                                                   >  1:342543/1‑150 (MQ=255)
                                 cttcAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGa                                                                                                               <  2:268616/150‑1 (MQ=255)
                                 cttcAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGa                                                                                                               <  2:342543/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:266549/150‑1 (MQ=255)
                                                   cACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCaa                                                                                             <  1:265913/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)
                                                                                                  caAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGAt                                              >  1:154213/1‑150 (MQ=255)
                                                                                                   aaGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttct                                                           <  1:381245‑M2/136‑88 (MQ=255)
                                                                                                   aaGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGttggctggctggtgttcgccctgttcctgctgctgccgctggtgatcgtggtatcgcaggcgctgaagaacggttttggcaccttct                                                           >  2:381245‑M2/1‑49 (MQ=255)
                                                                                                        cGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATcagca                                                        <  1:331459/134‑1 (MQ=255)
                                                                                                        cGATCAGGATACGCAGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATcagca                                                        >  2:331459/1‑134 (MQ=255)
                                                                                                                       gTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCaa                         <  1:309684/150‑1 (MQ=255)
                                                                                                                         gACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAaca                       <  2:68368/150‑1 (MQ=255)
                                                                                                                            gTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAACACGc                    >  1:190439/1‑150 (MQ=255)
                                                                                                                                  cTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAACACGCCGATgg              >  1:158054/1‑150 (MQ=255)
                                                                                                                                              gcggcgTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAACACGCCGATGGCGGTCGCGCCgg  >  2:64323/1‑150 (MQ=255)
                                                                                                                                              gcggcgTTGGCTGCGGCACTTGCAAGCATCGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAACACGCCGATGGCGGTCGCGCCgg  <  1:112910/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GAATATCGCCTCGAAGAAGGTGCCAAAACCGTTCTTCAGCGCCTGCGATACCACGATCACCAGCGGCAGCAGCAGGAACAGGGCGAACACCAGCCAGCCAAGGCCGATCAGGATACGCCGTGACGTGGCGCTGCCACGGCGGGCGGCGTTGGCTGCGGCACTTGCACTCAGGGATGAACTGGACATGCCGGCCTCCTTCAAGGGGTTTCGATGCGGCGCTGCAGCAGGTTGATCAGCAGCAGCAGGATGAAGGAAACCACCAGCATCAACACGCCGATGGCGGTCGCGCCGG  >  NC_002947/5895216‑5895507

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