New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4580793 =32 (1.370)3 (0.130) 3/244 NT 8.5% coding (2679/3321 nt) treSB bifunctional trehalose synthase B/maltokinase
?NC_002947 4580833 = 34 (1.520)coding (2639/3321 nt) treSB bifunctional trehalose synthase B/maltokinase

TGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTGCGCTTCGGCGTCGTTGCGGGTGGGC  >  NC_002947/4580686‑4580961
                                                                                                                                                   |                                                                                                                                
tGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAg                                                                                                                                <  2:261470/150‑1 (MQ=255)
    gctgcAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCg                                                                                                                            <  2:125655/150‑1 (MQ=255)
     cactgacgcattggccgaactcagcggctttgccgccttgcttggccagcgcctgggcgagatgcacctgctgctggccgcacccaccgaagacccggccttCCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCgg                                                                                                                           >  1:247230‑M1/103‑150 (MQ=255)
      tgcAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCggg                                                                                                                          >  1:208384/1‑150 (MQ=255)
         aGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCt                                                                                                                       >  1:569624/1‑150 (MQ=255)
              gtcgACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCggt                                                                                                                  <  2:523676/150‑1 (MQ=255)
                    cAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCg                                                                                                                    <  1:176773/142‑1 (MQ=255)
                    cAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCg                                                                                                                    >  2:176773/1‑142 (MQ=255)
                           tGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCcagc                                                                                                     >  2:486362/1‑150 (MQ=255)
                             gCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCg                                                                                                                            >  2:248959/1‑125 (MQ=255)
                             gCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCg                                                                                                                            <  1:248959/125‑1 (MQ=255)
                             gCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCcagcag                                                                                                   >  1:347811/1‑150 (MQ=255)
                                       tCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCAtc                                                                                         >  1:317221/1‑150 (MQ=255)
                                          aGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCg                                                                                      >  2:513180/1‑150 (MQ=255)
                                               gtCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAg                                                                                 <  1:503717/150‑1 (MQ=255)
                                               gtCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAg                                                                                 <  2:347811/150‑1 (MQ=255)
                                                  gcgATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACGAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAACGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCCGATGCATCTCGCACAGgca                                                                              >  2:14181/1‑149 (MQ=255)
                                                      tGCTGGGCCAGCAGGTCGAGCGCGTGCGTCCGTTCAGTGCTGATCTGTGCCCACCAGCGCTCGCTCTCCTGCGCGTCGCGGGGGCGCCGCGCGAAGCCCCGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTCCAGCTCGCCCAGCCGCTgg                                                                          <  2:506933/150‑1 (MQ=255)
                                                          gggTCAGCAGGGCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCcaa                                                                      <  2:246346/150‑1 (MQ=255)
                                                              cagcagGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCcaagcaa                                                                  >  2:253737/1‑150 (MQ=255)
                                                                  agGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATTGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCGGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAggcg                                                              <  1:205013/150‑1 (MQ=255)
                                                                         gggCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGc                                                       <  1:513180/150‑1 (MQ=255)
                                                                                 gTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGt                                               <  2:605563/150‑1 (MQ=255)
                                                                                      ttcagtGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGc                                          <  1:486362/150‑1 (MQ=255)
                                                                                           tGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAggc                                                               <  1:331384/124‑1 (MQ=255)
                                                                                           tGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAggc                                                               >  2:331384/1‑124 (MQ=255)
                                                                                              tGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGt                                  <  1:508004/150‑1 (MQ=255)
                                                                                                aTCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCa                                <  1:253737/150‑1 (MQ=255)
                                                                                                   ccggccttccagccgcgccccagcgatgcgcaggacagcgagcgctggAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAgtg                             >  1:247230‑M2/49‑150 (MQ=255)
                                                                                                    gtgCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAgtgt                            <  2:569624/150‑1 (MQ=255)
                                                                                                        acacCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTgcgc                        >  1:516377/1‑150 (MQ=255)
                                                                                                           ccagccgcgccccagcgatgcgcaggacagcgagcgctggAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGa                                                 >  1:164061‑M2/41‑122 (MQ=255)
                                                                                                           ccagccgcgccccagcgatgcgcaggacagcgagcgctggAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGa                                                 <  2:164061‑M2/82‑1 (MQ=255)
                                                                                                                    gctGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTGCGCTTCGGCGTCGtt            >  1:100370/1‑150 (MQ=255)
                                                                                                                     ctGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTGCGCTTCGGCGTCGTTg           >  1:290049/1‑150 (MQ=255)
                                                                                                                          cTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTGCGCTTCGGCGTCGTTGCGGGt      >  1:500070/1‑150 (MQ=255)
                                                                                                                              gcATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAgt                              <  1:538098/122‑1 (MQ=255)
                                                                                                                              gcATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAgt                              >  2:538098/1‑122 (MQ=255)
                                                                                                                              gcATCGCTGGGGCGAGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTGCGCTTCGGCGTCGTTGCGGGTGGGc  >  2:578927/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                
TGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGCGCTCGCTGTCCTGCGCATCGCTGGGGCGCGGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTGTGCGCTTCGGCGTCGTTGCGGGTGGGC  >  NC_002947/4580686‑4580961

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