New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3897448 =39 (1.320)4 (0.140) 4/244 2.8 8.2% intergenic (+9/‑58) yieH/cbrB phosphoenolpyruvate and 6‑phosphogluconate phosphatase/PRK09823 family inner membrane protein, creBC regulon
?NC_000913 3897501 = 53 (1.920)intergenic (+62/‑5) yieH/cbrB phosphoenolpyruvate and 6‑phosphogluconate phosphatase/PRK09823 family inner membrane protein, creBC regulon
Rejected: Frequency below/above cutoff threshold.

TCGACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAT  >  NC_000913/3897300‑3897596
                                                                                                                                                    |                                                                                                                                                    
agaaaaattggatgtattatcacccaaagaaccagcggttcttttgcgaagaatatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaa                                                                                                                                                      >  4:118935‑M1/149‑149 (MQ=255)
  gACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACt                                                                                                                                                    >  7:41449/1‑149 (MQ=255)
    cGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTcc                                                                                                                                                  <  5:78297/149‑1 (MQ=255)
      cAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCt                                                                                                                                                <  2:80007/149‑1 (MQ=255)
      attggatgtattatcacccaaagaaccagcggttcttttgcgaagaatatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaACTCCCt                                                                                                                                                <  4:87286‑M1/7‑1 (MQ=255)
       aGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCtt                                                                                                                                               <  2:87101/149‑1 (MQ=255)
            tGGAAGTGTTCTACTTCTGCGCCGCCCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACt                                                                                                                                           >  6:49907/1‑148 (MQ=255)
                            gaaccagcggttcttttgcgaagaatatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaACTCCCTTCACTTACCCCGCTTAAATTgg                                                                                                                          <  4:64093‑M1/29‑1 (MQ=255)
                                               gaagaatatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaACTcc                                                                                                                                                  <  6:75131‑M1/5‑1 (MQ=255)
                                               gaagaatatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaACTcc                                                                                                                                                  >  5:75131‑M1/102‑106 (MQ=255)
                                               gaagaatatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                                                       <  5:31860‑M1/48‑1 (MQ=255)
                                                   cGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAg                                                                                                   <  8:41449/149‑1 (MQ=255)
                                                    atatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgtcgtgaaacagacatatcaaACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAg                                                                                                   >  6:31860‑M1/97‑148 (MQ=255)
                                                     tatgtagaggacaagaaacagaacaccaattaacggtcctaaaactgcaaaataaagtccgtgatgtattacccgacgtgaaacagacatatcaaACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGtt                                                                                                 >  3:64093‑M1/96‑149 (MQ=255)
                                                             cccGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                                                         >  2:120016/1‑149 (MQ=255)
                                                              ccGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTc                                                                                        >  3:47806/1‑149 (MQ=255)
                                                                    gtacccaccTTTCACCATATTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTCCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAgg                                                                                                      <  8:61415/125‑1 (MQ=255)
                                                                    gTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTc                                                                                  <  1:120016/149‑1 (MQ=255)
                                                                    gTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTCCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAgg                                                                                                      >  7:61415/1‑129 (MQ=255)
                                                                     tCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCa                                                                                                                                             >  7:91475/1‑89 (MQ=255)
                                                                     tCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCa                                                                                                                                             <  8:91475/89‑1 (MQ=255)
                                                                                   cATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACa                                                                   <  1:11976/149‑1 (MQ=255)
                                                                                   cATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACa                                                                   >  7:64534/1‑149 (MQ=255)
                                                                                    aTCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                                        <  4:7688/111‑1 (MQ=255)
                                                                                    aTCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                                        >  3:7688/1‑111 (MQ=255)
                                                                                     tCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACAt                                                                  >  8:116977/1‑148 (MQ=255)
                                                                                         tCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATc                                                                 <  5:49907/145‑1 (MQ=255)
                                                                                             aGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACtt                                                         >  1:31131/1‑149 (MQ=255)
                                                                                              gTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACtt                                                         >  3:79651/1‑148 (MQ=255)
                                                                                                   cTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCAc                                                                                >  4:71636/1‑120 (MQ=255)
                                                                                                   cTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCAc                                                                                <  3:71636/120‑1 (MQ=255)
                                                                                                    tGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTg                                                  >  5:56927/1‑149 (MQ=255)
                                                                                                     gAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTg                                                  >  5:70015/1‑148 (MQ=255)
                                                                                                          gtgGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGttt                                            >  2:27668/1‑149 (MQ=255)
                                                                                                             gAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                          <  2:31131/148‑1 (MQ=255)
                                                                                                                                 ttACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCt                      <  4:47806/148‑1 (MQ=255)
                                                                                                                                     agCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTc                 >  8:97590/2‑149 (MQ=255)
                                                                                                                                       cATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCgg                                                                           >  2:12813/1‑89 (MQ=255)
                                                                                                                                       cATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCtt               <  8:64534/149‑1 (MQ=255)
                                                                                                                                            ttcttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                                          >  2:150995/1‑85 (MQ=255)
                                                                                                                                            ttcttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                                          <  1:150995/85‑1 (MQ=255)
                                                                                                                                               ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTc            <  4:5383/144‑1 (MQ=255)
                                                                                                                                               ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTc            >  3:5383/1‑144 (MQ=255)
                                                                                                                                                 cacaCTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTAc      <  7:97590/148‑1 (MQ=255)
                                                                                                                                                  acacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTAc      >  6:31718/1‑147 (MQ=255)
                                                                                                                                                  acacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTAc      <  5:31718/147‑1 (MQ=255)
                                                                                                                                                   cacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACata   >  8:39728/1‑149 (MQ=255)
                                                                                                                                                    ACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTc            >  6:60572‑M1/1‑139 (MQ=255)
                                                                                                                                                    ACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTc            <  5:60572‑M1/139‑1 (MQ=255)
                                                                                                                                                    ACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACatat  <  8:111465‑M1/149‑1 (MQ=255)
                                                                                                                                                    ACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACatat  <  4:106403‑M1/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAT  >  NC_000913/3897300‑3897596

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.