New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1164698 =32 (1.040)4 (0.150) 4/114 2.2 11.5% coding (604/1026 nt) nagZ beta N‑acetyl‑glucosaminidase
?NC_000913 1165027 = 34 (1.290)coding (933/1026 nt) nagZ beta N‑acetyl‑glucosaminidase
Rejected: Frequency below/above cutoff threshold.

AGATAATCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTACATG  >  NC_000913/1164961‑1165088
                                                                  |                                                             
aGATAATCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGAc                                                              >  2:518955/1‑68 (MQ=255)
  aTAATCTGTCGCCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAg                                                            >  2:232272/1‑68 (MQ=255)
  aTAATCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAg                                                            >  1:264723/1‑68 (MQ=255)
  aTAATCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAg                                                            >  2:746134/1‑68 (MQ=255)
    aaTCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGa                                                          >  2:63360/1‑68 (MQ=255)
    aaTCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGa                                                          >  2:517349/1‑68 (MQ=255)
     aTCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGaa                                                         >  1:278656/1‑68 (MQ=255)
            gccGGTCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAGCTGATgg                                                  >  2:269806/2‑68 (MQ=255)
             ccGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGa                                                 <  1:1111454/68‑1 (MQ=255)
             acatcactgtagatcacatgcgcaggcataatggcgtcgagttcattttcgcgACAGGAACTGATGGa                                                 <  1:257813‑M2/15‑1 (MQ=255)
                aTCAAGGCAGAACGTGTTTCACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTc                                              >  1:968235/1‑68 (MQ=255)
                       agatcacatgcgcaggcataatggcgtcgagtttattttcgcgACAGGAACTGATGGACTCGGCTCGc                                       <  2:159794‑M2/25‑1 (MQ=255)
                           aCGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGa                                   >  1:299598/1‑68 (MQ=255)
                           aCGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGa                                   >  1:31031/1‑68 (MQ=255)
                           aCGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGa                                   >  1:895487/1‑68 (MQ=255)
                           aCGTGTTACACGATAGTATCATAAAGGTTCATTTTCGCGACAGGAACAGATGGACTCGGCTCGCTGGa                                   >  2:389526/1‑68 (MQ=255)
                            cGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGaa                                  <  2:113983/68‑1 (MQ=255)
                            cGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCCGGag                                  >  2:562138/1‑67 (MQ=255)
                                  gcaggcataatggcgtcgagtttattttcgcgACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGAt                            >  1:844284‑M2/33‑68 (MQ=255)
                                  acacGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGAt                            <  1:291988/68‑1 (MQ=255)
                                    acGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCa                          >  1:629334/1‑68 (MQ=255)
                                     cGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAg                         >  2:896230/1‑68 (MQ=255)
                                            tggcgtcgagttcattttcgcgACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGt                  <  2:618067‑M2/46‑1 (MQ=255)
                                            aTCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGt                  <  2:805261/68‑1 (MQ=255)
                                            aTCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGt                  <  2:1010977/68‑1 (MQ=255)
                                            aTCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGt                  <  1:969376/68‑1 (MQ=255)
                                            aTCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGt                  >  1:45383/1‑68 (MQ=255)
                                               aTAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTg               >  2:536780/1‑68 (MQ=255)
                                                 aaaGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGaa             >  2:374293/1‑68 (MQ=255)
                                                 aaaGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGaa             >  2:558589/1‑68 (MQ=255)
                                                 aaaGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGaa             >  2:83384/1‑68 (MQ=255)
                                                   aGGTTCATTTTCGCGACAGGAACTGAAGGCATCGGCTCGCTGGAAAGCGAACAGCACCCGTCTGAATc           >  2:309613/1‑68 (MQ=255)
                                                    ggTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCTATCAGCACCCGTCTGAATCa          >  1:579453/1‑68 (MQ=255)
                                                    ggTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCa          >  2:305377/1‑68 (MQ=255)
                                                    ggTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCa          <  1:411881/68‑1 (MQ=255)
                                                     gTTCATTTTCGCGACAGTAACTGATGGACTCGTCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAg         <  1:1142898/68‑1 (MQ=255)
                                                     gTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAg         <  1:877748/68‑1 (MQ=255)
                                                        cATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTa      <  2:991152/68‑1 (MQ=255)
                                                           tttcgcgACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTACAt   >  1:1137299‑M2/8‑68 (MQ=255)
                                                            ttcgcgACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCGCCCGTCTGAATCAGTTACATg  >  1:219331‑M2/7‑68 (MQ=255)
                                                                  |                                                             
AGATAATCTGTCACCGATCAAGGCAGAACGTGTTACACGTTTGTATCATAAAGGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTACATG  >  NC_000913/1164961‑1165088

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