New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5750393 =15 (0.870)4 (0.240) 4/268 1.6 20.8% coding (1226/1455 nt) thiI tRNA sulfurtransferase
?NC_002947 5750415 = 16 (0.970)coding (1204/1455 nt) thiI tRNA sulfurtransferase

GCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGA  >  NC_002947/5750256‑5750513
                                                                                                                                         |                                                                                                                         
gCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAgcg                                                                                                              >  1:178303/1‑150 (MQ=255)
    gTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTc                                                                                                          <  2:228255/150‑1 (MQ=255)
                atggagcacgaagaaaagcagttcgacatggccgtgctggagcgcgctcttgagcgcgccaagttcatttccatcgatcatgtgatcgatgagctgggcaaggacatcgaaatcgaggaagTGGCCTGGCAGCGCCTCGGCCACTTACTc                                                                                              >  1:231190‑M1/122‑150 (MQ=255)
                   gTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGAt                                                                                           <  2:191714/150‑1 (MQ=255)
                         gTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGAt                                                                                     <  1:277262/150‑1 (MQ=255)
                                   ttGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGccc                                                                           >  2:63332/1‑150 (MQ=255)
                                        gTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTc                                                                      >  1:259423/1‑150 (MQ=255)
                                                  tAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATcaca                                                            >  2:165183/1‑150 (MQ=255)
                                                                cTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATgagctg                                              <  1:63332/150‑6 (MQ=255)
                                                                                   ttcACCAGAGGTTCGTCTTCCTGTGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                           <  2:135276/148‑1 (MQ=255)
                                                                                      ccatcgatcatgtgatcgatgagctgggcaaggacatcgaaatcgaggaagTGGCCTGGCAGCGCCTCGGCCACTTCCt                                                                                               <  1:398107‑M1/28‑1 (MQ=255)
                                                                                      ccatcgatcatgtgatcgatgagctgggcaaggacatcgaaatcgaggaagTGGCCTGGCAGCGCCTCGGCCACTTCCt                                                                                               >  2:398107‑M1/52‑79 (MQ=255)
                                                                                                  tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATga                                                  <  1:199522/112‑2 (MQ=255)
                                                                                                  tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATga                                                  >  2:199522/1‑111 (MQ=255)
                                                                                                     tCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGcc         <  1:187452/150‑1 (MQ=255)
                                                                                                          ggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTc    <  1:165183/150‑1 (MQ=255)
                                                                                                            gctgggcaaggacatcgaaatcgaggaagTGGCCGAGGCGCTGCCAGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGa  <  2:130283‑M1/121‑1 (MQ=255)
                                                                                                                 gggTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGAt                                                    >  2:75107/1‑95 (MQ=255)
                                                                                                                 gggTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGAt                                                    <  1:75107/95‑1 (MQ=255)
                                                                                                                                         |                                                                                                                         
GCAAGTACTGGCGCGTGGGGTCCAGGTGCTTGAACTTGCTGTTGATGGCGTAGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGA  >  NC_002947/5750256‑5750513

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