New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4583475 =36 (0.600)4 (0.080) 3/226 NT 15.8% intergenic (‑48/+123) treSB/glgE bifunctional trehalose synthase B/maltokinase/alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase
?NC_002947 4583576 = 12 (0.240)intergenic (‑149/+22) treSB/glgE bifunctional trehalose synthase B/maltokinase/alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase

GTGCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCT  >  NC_002947/4583428‑4583713
                                                                                                                                                    |                                                                                                                                         
gtgCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCTGGTGTTCGCGGATAAACCCGCTCCTACgggggg                                                                                                                                           >  3:267753/1‑149 (MQ=255)
   cGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACggggggggg                                                                                                                                        >  3:25712/1‑149 (MQ=255)
          tttCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAg                                                                                                                                 >  1:348662/1‑149 (MQ=255)
          tttCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTAGCGGATAAACCAGCTCCTACGGGGGGGGGGGcgac                                                                                                                                  >  3:250781/1‑145 (MQ=255)
                                 cGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGGCTGAACCTTGGTAACGGGCttt                                                                                                          >  6:81619/1‑149 (MQ=255)
                                  gAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCCCACGGTGAACGGTTGTAAGGGGCtttt                                                                                                         >  4:102557/1‑149 (MQ=37)
                                          tCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGGTTTTTGGATgcg                                                                                                 >  1:52832/1‑149 (MQ=255)
                                                  CCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAAACTTTTTAACCGGCGTTTTAGATgcgccaaaaga                                                                                         <  6:284183‑M1/149‑6 (MQ=255)
                                                     aggagcgggtttatccgcgaacaccggcgaagccggtgccatgcacctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacggggggGGGGGGCACAGACTGAACCTTTTTAACGGGCTTTTTCGATGCGCCACATGCCg                                                                                       >  8:15825‑M2/96‑148 (MQ=255)
                                                      aggagcgggtttatccgcgaacaccggcgaagccggtgccatgcacctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacgggggGGGGGGCCCCGAGGGACACTTGTTACCGGGTTTTTTCGATGCCCCAGATCCGGaa                                                                                     >  6:326039‑M2/95‑149 (MQ=255)
                                                        GAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACCTTCTAAACATTGTTAACGGGCTTTTACAATGCGCAAGATGCCGAACg                                                                                   <  4:261169‑M1/149‑1 (MQ=255)
                                                        GAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTACATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGGACAGGCGGCAACATGTTAACGGGCGTTTTGGAGGCCCCAAGAGCCGAACg                                                                                   <  8:14446‑M1/149‑1 (MQ=255)
                                                          gcgggtttatccgcgaacaccggcgaagccggtgccatgcacctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacgggggGGGGGGCACCGTGTGAACCTTGTTAACGGGGTTTTTCGATGCGCCAGATGCCGAACgg                                                                                  >  6:199151‑M2/91‑148 (MQ=255)
                                                           cgggtttatccgcgaacaccggcgaagccggtgccatgcacctcggcaccgcggtgtcatcttcgcgggtaaacccgctcctacgggggGGGGGGCAAAGTATGAACCTTGTTAACGGGGTTTTTAGATGCGCCAGATGCCGCACgg                                                                                  >  6:326870‑M2/90‑147 (MQ=255)
                                                                    ccGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAggg                                                                      >  7:16927/1‑149 (MQ=255)
                                                                    ccGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGGATGAACCTTGTTAACGGGCTTTTTCGATGCGCCAAATGCCGAACGCAAGGTGCCAggg                                                                      >  8:202958/1‑149 (MQ=38)
                                                                                     cgGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCAtt                                                     <  8:16927/149‑1 (MQ=255)
                                                                                       gTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGCACAGTCTGAACCTTGTTAACTGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTgg                                                   <  7:202958/149‑1 (MQ=255)
                                                                                           cGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCCCCGGCTGAACCGTGTTAACGGGc                                                                                                             >  5:262869/1‑88 (MQ=21)
                                                                                                        aCCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCAc                                    >  2:266234/1‑148 (MQ=255)
                                                                                                        aCCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACg                                   >  3:229169/1‑149 (MQ=255)
                                                                                                                              ggcaaccccgcccccacgggggGGCGGCCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAg             <  5:199151‑M2/127‑1 (MQ=255)
                                                                                                                              gATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAg             <  1:266234/149‑1 (MQ=255)
                                                                                                                                      ccgctcctacggggGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCAc     >  1:216399‑M2/15‑149 (MQ=255)
                                                                                                                                         tcctacgggggGGGGGACACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCt  >  7:283152‑M2/12‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
GTGCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCT  >  NC_002947/4583428‑4583713

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