New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4200145 =28 (0.910)4 (0.160) 4/110 2.1 18.9% intergenic (+51/‑136) yjaG/hupA DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
?NC_000913 4200182 = 11 (0.430)intergenic (+88/‑99) yjaG/hupA DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Rejected: Frequency below/above cutoff threshold.

TAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAG  >  NC_000913/4200078‑4200212
                                                                   |                                                                   
taATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCt                                                                     <  2:300683/68‑1 (MQ=255)
  aTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCttt                                                                   >  2:839742/1‑68 (MQ=255)
  aTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCttt                                                                   >  1:651633/1‑68 (MQ=255)
    ttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCg                                                                 >  2:114922/1‑68 (MQ=255)
    atttCAGCAATAAGACCAGAAAACGTGATATAACGCATGATTTGTCGTACCTGGAGTCTTCCCTTTCg                                                                 >  1:186284/2‑68 (MQ=255)
    atcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCg                                                                 <  2:315605‑M1/5‑1 (MQ=255)
    atcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCg                                                                 <  1:15055‑M1/5‑1 (MQ=255)
    atcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCg                                                                 <  2:499941‑M1/5‑1 (MQ=255)
         gtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGccccc                                                            >  2:900710‑M1/59‑68 (MQ=255)
         gtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGccccc                                                            >  1:315605‑M1/59‑68 (MQ=255)
          gCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccccc                                                           >  1:914694/1‑68 (MQ=255)
                aGACCAGTATACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCTCCGTCTgg                                                     <  2:606928/68‑1 (MQ=255)
                aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGATTCTTCCCTTTCGCCCCCCGTCTgg                                                     <  2:160551/68‑1 (MQ=255)
                aGACCAGAAAACGCGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGCCTgg                                                     <  2:1053574/68‑1 (MQ=255)
                 ggccAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGt                                                    >  2:1101507/3‑68 (MQ=255)
                 gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGt                                                    >  1:444294/1‑68 (MQ=255)
                 gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGt                                                    <  2:739371/68‑1 (MQ=255)
                   aatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGCCCCCCGTCTGGTCt                                                  >  1:1128879‑M1/49‑68 (MQ=255)
                        aaaCGTGATTTAACGCCTGATTTGTCGTACCTGGCGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAtt                                             >  2:355484/1‑68 (MQ=255)
                               aTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                     <  2:946776/53‑1 (MQ=255)
                               aTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                     >  1:946779/1‑53 (MQ=255)
                                   aaCGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                  <  1:1101512/68‑1 (MQ=255)
                                                cgcaccgatagccacttttTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt                     <  2:913129‑M1/49‑1 (MQ=255)
                                                 cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACa                                               >  2:475540/1‑41 (MQ=255)
                                                 cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACa                                               <  1:475540/41‑1 (MQ=255)
                                                 cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATc                    >  1:1094645/1‑68 (MQ=255)
                                                   accgatagccacttttTTTCGCAACCCGTCAGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCgg                  <  2:49498‑M1/52‑1 (MQ=255)
                                                      cTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGc               >  2:15055/1‑68 (MQ=255)
                                                       cggAGTCTTCCCTTTCGCCCCCCGTCTTGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGcgcg              <  2:670808/67‑1 (MQ=255)
                                                          aGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt           >  2:721509/1‑68 (MQ=255)
                                                          aGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGCGCgcgt           <  2:30093/68‑1 (MQ=255)
                                                            tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt         <  2:717562/68‑1 (MQ=255)
                                                            tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt         >  1:1090300/1‑68 (MQ=255)
                                                            tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt         <  1:1062266/68‑1 (MQ=255)
                                                            tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGCGCGTGTAt         <  2:1099555/68‑1 (MQ=255)
                                                            tCTTCCCTTTCGCCCCCCGGCTGTTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt         <  2:294464/68‑1 (MQ=255)
                                                             acttttTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATg        <  2:792517‑M1/62‑1 (MQ=255)
                                                                   TTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag  >  1:1059639‑M1/1‑68 (MQ=255)
                                                                   |                                                                   
TAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAG  >  NC_000913/4200078‑4200212

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