New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 153686640 (1.340)4 (0.140) 4/266 3.6 9.9% coding (444/696 nt) narW nitrate reductase 2 (NRZ), delta subunit (assembly subunit)
?NC_000913_3_bme_pgi = 1536879 36 (1.290)coding (431/696 nt) narW nitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCT  >  NC_000913_3_bme_pgi/1536732‑1537020
                                                                                                                                                   |                                                                                                                                             
ttCCCCTGTTCCTCTTCCCACACTGCATCAATCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTc                                                                                                                                              <  2:1880/149‑1 (MQ=255)
  cACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgc                                                                                                                                            <  1:194431/149‑1 (MQ=255)
    ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGAGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTCGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                                                                          <  1:487770/149‑1 (MQ=255)
                cccACACGGCATCAAGCCGCTGGCGGGAGTGCGCGCGCTCTTCGGTGTTTCCTTGTTTTGTGATACTGTCTCTTGAAGGGTAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCgc                                                                                                                              <  2:347674/149‑2 (MQ=255)
                   acacTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCg                                                                                                                           >  1:88697/1‑149 (MQ=255)
                      cTGCATCAAGCGCCTGGCGGGTGTCGTCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCGCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGGGCATCAAACAACGCGTACCACGGTCAATCGCGTTGTTTTAAGCGACCGCCa                                                                                                                        >  1:10409/1‑149 (MQ=255)
                         cATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc                                                                                                                     <  2:522330/149‑1 (MQ=255)
                         cATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc                                                                                                                     <  1:523757/149‑1 (MQ=255)
                                cgcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGa                                                                                                              <  2:201231/149‑1 (MQ=255)
                                 gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCCGGGCGAg                                                                                                             >  2:452102/1‑149 (MQ=255)
                                   cTGGCGGGTGTCATCGCGCTCTTCGCAGTTTACTTGTTTTATGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACCACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCCACCGCCAAGCAGGGCGAggg                                                                                                           >  2:401986/1‑148 (MQ=255)
                                    tGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAt                                                                                                          <  2:229155/149‑1 (MQ=255)
                                      gCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCg                                                                                                        <  1:67382/149‑1 (MQ=255)
                                        gggTGTCACCGCGCCTTTCGCTGTTTACTTGTTTTGGGAAAATGTCACTTGAAAGGGAGCTTCCCGCCCGTTGCAGCAGAGCATCAAACCACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgc                                                                                                      <  2:24743/149‑1 (MQ=255)
                                        ctgtgtCATCGCGGCCTTCGCAGAATACTTGTTATGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCCCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgc                                                                                                      <  1:14107/147‑1 (MQ=255)
                                               aTCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCCACCGCCAAGCCGGGCGAGGATCCGCGCCACAt                                                                                               >  2:47775/1‑149 (MQ=255)
                                                  gcgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGa                                                                                            <  2:205278/149‑1 (MQ=255)
                                                             tGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttct                                                                                 <  1:223106/149‑1 (MQ=255)
                                                                           gtgCCAGTGGCACGTGAAAGGGAGCTTCCCGCCATTTGCAGCAGAGCATCAACCAACGCGTACCACGGTGCCTCGCTGTGTTTCACGCGACCGCCAAGCAGGGCGAGGATCAGCGCTACATGGATCAATCCTTTTTTCGCCTGATCGTc                                                                   <  2:43482/149‑1 (MQ=255)
                                                                            tgACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCa                                                                  >  1:6496/1‑149 (MQ=255)
                                                                                     cACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAATCGCCCGCCAAGCAGGGCGAGGTTCGGCGCAACATTGAGCAATCCTTCTTTCATCTGATCGGCAGGCAGtcag                                                         >  1:55427/1‑146 (MQ=255)
                                                                                           aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                   <  1:355421/149‑1 (MQ=255)
                                                                                                gAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTc                                              <  2:359922/149‑1 (MQ=255)
                                                                                                    ttCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCCGTACGCCTAAATACCCCAGa                                          >  2:407165/1‑149 (MQ=255)
                                                                                                    ttCCCGCCAGGTGCAGCAGAGCATCAAACAACGCGTACCCCGGTCCCTCGGGGTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGa                                          <  2:505383/149‑1 (MQ=255)
                                                                                                           cAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCg                                   <  1:268426/149‑1 (MQ=255)
                                                                                                              ttGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTCAATACTCCAGATACAGCGGTa                                <  1:473237/149‑1 (MQ=255)
                                                                                                              ttGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGGCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTa                                <  2:247680/149‑1 (MQ=255)
                                                                                                                 cagcagAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTCCGCTTAAATACTCCAGATACCGCGGTACAt                             >  2:243730/1‑149 (MQ=255)
                                                                                                                     agagCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGACCGTCAGGCAGTACGCTTAAATACCCCAGATACAGCGGTAGATAATc                         >  2:492513/1‑149 (MQ=255)
                                                                                                                          ttaaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCGCCATTGAGCAATCCTTCTTTCCCCCGATCGTCGGGCAGTACGCTTAAATCCTCCAGCGACAGAGGTAGATACACGGGCa                    >  2:295198/4‑149 (MQ=21)
                                                                                                                          aTCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                    <  1:402663/149‑1 (MQ=255)
                                                                                                                                caacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc              >  2:439612/1‑149 (MQ=255)
                                                                                                                                 aacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCg             >  2:145948/1‑149 (MQ=255)
                                                                                                                                      gTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGAGCAGTTCCCGACAAt        >  1:210358/1‑149 (MQ=255)
                                                                                                                                       tACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAg                               >  2:435371/1‑125 (MQ=255)
                                                                                                                                       tACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAg                               <  1:435371/125‑1 (MQ=255)
                                                                                                                                          cACGGTGCCTCGCGTTGTTTGAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAg    <  1:47775/149‑1 (MQ=255)
                                                                                                                                            cGGTGCCTCGCGTTGTTTTAAGCGGCCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTTAGGCAGTACGCTTAAATTCTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCt  <  2:296502/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCT  >  NC_000913_3_bme_pgi/1536732‑1537020

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