New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4467063 =37 (1.010)8 (0.230) 7/234 2.7 18.0% coding (184/948 nt) treR trehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
?NC_000913 4467094 = 37 (1.040)coding (153/948 nt) treR trehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
Rejected: Frequency below/above cutoff threshold.

TTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAA  >  NC_000913/4466920‑4467210
                                                                                                                                               |                                                                                                                                                   
ttCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCttt                                                                                                                                                <  5:92482/149‑1 (MQ=255)
      aaCTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCAc                                                                                                                                          <  8:77950/149‑1 (MQ=255)
         tcaaagatatcgcgcgcttaagcggcgtggggaaatctacagtttcccgggtgctgaataacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCa                                                                                                                                       <  1:82435‑M1/15‑1 (MQ=255)
                   gAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAgcgcgcg                                                                                                                              >  5:54059/1‑148 (MQ=255)
                         cttaagcggcgtggggaaatctacagtttcccgggtgctgaataacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCg                                                                                                                       >  2:82435‑M1/119‑149 (MQ=255)
                           aCTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCgg                                                                                                                      >  1:42942/1‑148 (MQ=255)
                              ttCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATa                                                                                                                                     >  1:413/1‑130 (MQ=255)
                              ttCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATa                                                                                                                                     <  2:413/130‑1 (MQ=255)
                                ccATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgg                                                                                                                 <  6:54059/148‑1 (MQ=255)
                                 catcatGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAggg                                                                                                                >  2:59967/1‑148 (MQ=255)
                                   tcatGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGa                                                                                                              >  5:30577/1‑148 (MQ=255)
                                           gggTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCAt                                                                                                      <  2:42942/148‑1 (MQ=255)
                                              tacagtttccctggtgctgaataacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCt                                                                                                                                                  >  8:30656‑M1/98‑101 (MQ=255)
                                              tacagtttcccgggtgctgaataacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCt                                                                                                                                                  <  7:30656‑M1/4‑1 (MQ=255)
                                                 agtttcccgggtgctgaataacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgg                                                                                                                 <  4:93832‑M1/37‑1 (MQ=255)
                                                 agtttcccgggtgctgaataacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgg                                                                                                                 >  3:93832‑M1/95‑131 (MQ=255)
                                                      ttGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATg                                                                                                     <  2:108999/138‑1 (MQ=255)
                                                      ttGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATg                                                                                                     >  1:108999/1‑138 (MQ=255)
                                                            aTAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGATGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCtt                                                                                    >  5:3516/1‑149 (MQ=255)
                                                            aTAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCtt                                                                                    <  1:59967/149‑1 (MQ=255)
                                                             tAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTc                                                                                   <  6:30577/149‑1 (MQ=255)
                                                              aGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTc                                                                                   >  1:204778/1‑148 (MQ=255)
                                                                   taacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgg                                                                                                                 >  8:130117‑M1/77‑113 (MQ=255)
                                                                   taacgaaagcggcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgg                                                                                                                 <  7:130117‑M1/37‑1 (MQ=255)
                                                                       gcagcaTGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTcg                                                                         >  6:41754/1‑149 (MQ=255)
                                                                         agcaTGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTcgcg                                                                       <  2:204778/149‑1 (MQ=255)
                                                                            aTGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCgg                                                                      <  3:134253/147‑1 (MQ=255)
                                                                            aTGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGt                                                                    >  4:150216/1‑149 (MQ=255)
                                                                            aTGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCAATGCTGATTCATCAAAGCTTCAACACGCTCGCgg                                                                      >  4:134253/1‑147 (MQ=255)
                                                                              gcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCgg                                                                                                                      <  6:14902‑M1/32‑1 (MQ=255)
                                                                              gcgtgagccagctcacccgcgagcgtgttgaagcagtgatgaatcagcatggattttccccttccCGCTTTGCCCACGCATAGCGCGCGCAGAGCgg                                                                                                                      >  5:14902‑M1/66‑97 (MQ=255)
                                                                                         cGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCtt                                                                                    <  2:3828/120‑1 (MQ=255)
                                                                                         cGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCtt                                                                                    >  1:3828/1‑120 (MQ=255)
                                                                                          gagaTTTTCTGACAACGAATCCAGATGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCCAGCTCAccc                                                      <  6:3516/149‑3 (MQ=255)
                                                                                                      cAACGAATCCAGCCGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGATCTGGCTCACGCCGCTTTAGTTAt                                          >  3:45905/1‑149 (MQ=255)
                                                                                                               cAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAccc                                 >  1:38328/1‑149 (MQ=255)
                                                                                                                aGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCg                                >  5:81536/1‑149 (MQ=255)
                                                                                                                      gTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgggg                                                                                                               <  1:128823/64‑1 (MQ=255)
                                                                                                                      gTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAgggg                                                                                                               >  2:128823/1‑64 (MQ=255)
                                                                                                                      gTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCt                                                               >  2:77806/1‑112 (MQ=255)
                                                                                                                      gTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCt                                                               <  1:77806/112‑1 (MQ=255)
                                                                                                                          cAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                      <  6:81536/149‑1 (MQ=255)
                                                                                                                                          tttATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgccgc      >  8:82534/1‑149 (MQ=255)
                                                                                                                                            tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                 >  2:176206/1‑104 (MQ=255)
                                                                                                                                            tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                 <  1:176206/104‑1 (MQ=255)
                                                                                                                                              tCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTaa  <  3:150216/149‑1 (MQ=255)
                                                                                                                                              tCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTaa  <  2:38328/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTCGGCAACTAATTGCGGGGAAAACTGACTTTCCATCATGATTGGGTCGTAACCTTGTTCATAGAACGCTGGCAGCATGGTTTGAACGGCGAGATTTTCTGACAACGAATCCAGACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAA  >  NC_000913/4466920‑4467210

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