New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1759424 =30 (0.940)6 (0.200) 5/242 2.9 17.6% coding (1100/1221 nt) sufS cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
?NC_000913 1759462 = 28 (0.930)coding (1062/1221 nt) sufS cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

CGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGG  >  NC_000913/1759276‑1759563
                                                                                                                                                    |                                                                                                                                           
cGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGCTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATg                                                                                                                                             >  6:81668/1‑149 (MQ=255)
   acagaatctgatgcattatgcgctatcacagctggaatctgtaccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCa                                                                                                                                          <  6:100202‑M1/4‑1 (MQ=255)
   acagaatctgatgcattatgcgctatcacagctggaatctgtaccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCa                                                                                                                                          <  4:78787‑M1/4‑1 (MQ=255)
    cagaatctgatgcattatgcgctatcacagctggaatctgtaccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCAt                                                                                                                                         >  1:22401‑M1/145‑149 (MQ=255)
       aaaGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTgc                                                                                                                                      >  3:116146/1‑149 (MQ=255)
            gatgcattatgcgctatcacagctggaatctgtaccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGt                                                                                                                                 <  3:73841‑M1/13‑1 (MQ=255)
            gatgcattatgcgctatcacagctggaatctgtaccggatatcactctctatgtcccacaaaacaggcttggagttattgattttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGt                                                                                                                                 <  5:106470‑M1/13‑1 (MQ=255)
                  gCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTc                                                                                                                           <  4:116146/149‑1 (MQ=255)
                                 cAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGc                                                                                                              >  4:3/1‑147 (MQ=255)
                                    cAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGc                                                                                                         >  3:54982/1‑149 (MQ=255)
                                      aatctgtaccggatctcactctctatggcccacaaaacaggctttgcgttattgcttttaatctcgggaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGTGATGTCCGGTACTCACAGCAATGCCg                                                                                                       >  6:106470‑M1/111‑149 (MQ=255)
                                        tctgtaccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGcc                                                                                                        >  2:133888‑M1/109‑146 (MQ=255)
                                        tctgtaccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGcc                                                                                                        <  1:133888‑M1/38‑1 (MQ=255)
                                         ggTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCATaa                                                                                                    <  3:3/149‑1 (MQ=255)
                                              ccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATc                                                                                               >  4:73841‑M1/103‑149 (MQ=255)
                                              aTACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGCTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATc                                                                                               <  5:81668/149‑1 (MQ=255)
                                                                   aGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAAcatca                                                                            >  7:16765/1‑147 (MQ=255)
                                                                                      gggAGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGAATTGCGCAGTGATGTCCGGTACGCACAGAAATGCAGTAATTATCGAGAAAAATGCCAACATCATATGCGTGGTGTTTAcc                                                           >  8:24660/1‑145 (MQ=255)
                                                                                           ttATACATGGCCAGCGACGCGCGACACTTCGCAGGGACGTTGTAATAGGATATCAATGGCATTTCGCAGTTATGTCCGGTACGCACAGCAATGGCGTAATTATCGAGAAAACTTCCAACATCATAGGCGTTGTGTTTACCGAGATTaa                                                   <  7:24660/148‑1 (MQ=255)
                                                                                               aCATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGc                                               >  7:6398/1‑148 (MQ=255)
                                                                                                cATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTTCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCaa                                             >  5:54744/1‑149 (MQ=255)
                                                                                                cATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCaa                                             <  8:16765/149‑1 (MQ=255)
                                                                                                   ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATa                                                                          <  6:106467/117‑1 (MQ=255)
                                                                                                     ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCa                                                                                  <  7:95441/107‑1 (MQ=255)
                                                                                                     ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCa                                                                                  >  8:95441/1‑107 (MQ=255)
                                                                                                     ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACg                                        >  3:84256/1‑149 (MQ=255)
                                                                                                      cAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGc                                       >  2:61903/1‑149 (MQ=255)
                                                                                                              cccGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTTCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTg                               <  6:54744/149‑1 (MQ=255)
                                                                                                               ccGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGt                              >  1:74661/1‑149 (MQ=255)
                                                                                                                          cAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCa                   >  5:44321/1‑149 (MQ=255)
                                                                                                                           aGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCa                   >  5:100202/1‑148 (MQ=255)
                                                                                                                                 gTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATagagag            >  7:85182/1‑149 (MQ=255)
                                                                                                                                     tAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGt           >  3:78787/1‑146 (MQ=255)
                                                                                                                                     tAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgag        <  1:61903/149‑1 (MQ=255)
                                                                                                                                      aaTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAg                                                >  2:35472/1‑108 (MQ=255)
                                                                                                                                      aaTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAg                                                <  1:35472/108‑1 (MQ=255)
                                                                                                                                           gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg  <  6:44321/149‑1 (MQ=255)
                                                                                                                                           gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg  <  2:74661/149‑1 (MQ=255)
                                                                                                                                           gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg  <  8:6398/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGG  >  NC_000913/1759276‑1759563

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