New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 670713 =18 (0.900)3 (0.180) 3/110 1.6 19.4% coding (893/1032 nt) holA DNA polymerase III, delta subunit
?NC_000913 670903 = 10 (0.600)coding (703/1032 nt) holA DNA polymerase III, delta subunit
Rejected: Frequency below/above cutoff threshold.

CTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAC  >  NC_000913/670648‑670774
                                                                 |                                                             
ctgactgcCCGTAATCTTGTTTGAGGGTGAGTTCCGATCGTGTCAGGAGTAGCACGGCCTGTCGTAAc                                                             >  2:144313/1‑68 (MQ=255)
  gactgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTg                                                           >  2:354878/1‑68 (MQ=255)
   ttttcattgggctgatgctttgttgatgggaaaaagtaagcgcgcattgcatattcttcagcAACTGc                                                          >  1:356687‑M1/63‑68 (MQ=255)
   actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGc                                                          <  2:499097/68‑1 (MQ=255)
    ctgacCGTAATCTTGTTTGAGGGTGAGTTCCGTCCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCg                                                         <  2:160032/68‑1 (MQ=255)
    ctgacCGTAATCTTGTTAGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCg                                                         <  2:718925/68‑1 (MQ=255)
      tcattgggttgatgctttgttgatgggaaaaagtaagagcgcatagcatattcttcagcAAATGCGTc                                                       <  2:181142‑M1/9‑1 (MQ=255)
        attgggttgatgctttgttgatgggaaaaagtaagcgcgcattgcatattcttcagcAACTGCGTCTg                                                     <  1:70245‑M1/11‑1 (MQ=255)
               cTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAg                                              >  1:418392/1‑68 (MQ=255)
                 tGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCg                                            >  2:534981/1‑68 (MQ=255)
                   tttGAGGGCGAGTTCCGTTCGTGTCAGGAGCTGCACGGCCTGACGTAACTGCGGCTGACTTAAGCGAt                                          <  1:100220/68‑1 (MQ=255)
                    ttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAtt                                         >  1:82903/1‑68 (MQ=255)
                    ttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAtt                                         >  2:187603/1‑68 (MQ=255)
                        gttgatgggaaaaagtaagcgcgcattgcatattcttcagcAACTGCGTCTGACTTAAGCGATTTAAc                                     <  2:470999‑M1/27‑1 (MQ=255)
                        gggTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAAc                                     >  2:496910/1‑68 (MQ=255)
                         ttgatgggaaaaagtaagcgcgcattgcatattcttcagcAACTGCGTCTGACTTAAGCGATTTAACg                                    >  2:70245‑M1/41‑68 (MQ=255)
                         ggTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACg                                    <  1:126641/68‑1 (MQ=255)
                         ggTGAGTTCCGGTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACg                                    <  1:219357/68‑1 (MQ=255)
                         ggTGAGTTCCGATCGTGTCAGGAGTTGCACGGCATGCCGTAGCTGCGTCTGAATTAAGCGATTTAACg                                    >  1:450855/1‑68 (MQ=255)
                                         gtCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGcc                    <  1:354877/68‑1 (MQ=255)
                                                      tattcttcagcAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGc       <  2:57416‑M1/57‑1 (MQ=255)
                                                      cGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGc       <  2:287320/68‑1 (MQ=255)
                                                           tGACGTAACTTCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAc  <  2:109393/68‑1 (MQ=255)
                                                           tGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAc  <  2:602638/68‑1 (MQ=255)
                                                                 |                                                             
CTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAC  >  NC_000913/670648‑670774

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