New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 671548 =57 (1.100)4 (0.080) 3/260 7.6 7.1% coding (58/1032 nt) holA DNA polymerase III, delta subunit
?NC_000913 671580 = 52 (1.090)coding (26/1032 nt) holA DNA polymerase III, delta subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGGCAGCTTACGAA  >  NC_000913/671438‑671720
                                                                                                                                              |                                                                                                                                            
gtCAATGGTAAAAGTGTGGTGTTCTGCGAATCCTTGAGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAAGAACAGAGGATCGTTACCAAGTCAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGa                                                                                                                                              <  1:382973/142‑1 (MQ=255)
aTCAATGGAAAAAGTGTGGTTTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGa                                                                                                                                              <  1:704386/143‑1 (MQ=255)
aTCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGa                                                                                                                                              <  1:589768/143‑1 (MQ=255)
     tGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTAGCAAGTAAAAGATACGCCGCGCGCAGCCGGTCATTGAGCTGCGCGCGGAgtggc                                                                                                                                         >  2:950782/1‑140 (MQ=255)
            aGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCATTCCTTGAGCTGCGCGCGGAGTGGTTCCGGGt                                                                                                                                  >  2:351760/1‑143 (MQ=255)
              tgtgGTGTTCTTCAAATCCTTGTGCCGCAACAACCCGACGAACAACGCCCCGGCTGCCCCGCAACAACACAGCAGCCTGACCAAGAAAAACATACGCCGCCCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTAc                                                                                                                                <  1:726376/143‑1 (MQ=255)
                 ggtgTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACacc                                                                                                                             >  2:19007/1‑141 (MQ=255)
                 ggtgTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAAc                                                                                                                             <  2:812037/143‑1 (MQ=255)
                   tgtTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTCAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGCGCTGTTCCGGGTACCACCg                                                                                                                           >  2:515765/1‑143 (MQ=255)
                    gtTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCCCACCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACACCCGa                                                                                                                          >  1:328800/1‑143 (MQ=255)
                       cttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAatc                                                                                                                       >  1:154751/1‑143 (MQ=255)
                        ttcGAATCCTTGTGCCGCAGCTACCTGCCGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCCTCATTGAGCTGCGCGCGGAGCTGTCCCCGGTACAACCAActca                                                                                                                      >  2:484621/1‑143 (MQ=255)
                         tcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAatcgt                                                                                                                     >  1:76690/1‑141 (MQ=255)
                            aaTCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGGACAACCGAATCATCAg                                                                                                                  >  1:543073/1‑143 (MQ=255)
                                   gtgCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCa                                                                                                           >  1:579509/1‑143 (MQ=255)
                                        gCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCACGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGGTCTTCCGGGTCCAACCGAATCATCAGTTACCCCGCgtg                                                                                                      >  1:502928/1‑143 (MQ=255)
                                        gCAGCTACATGTCGAACTGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGTGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCgtg                                                                                                      <  2:444871/143‑1 (MQ=255)
                                             tACCTGACGAACAGCGTCCTGGCTTTCCGGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGga                                                                                                 <  2:950892/143‑1 (MQ=255)
                                             tACCTGACGAACAGCGTCCGGGCTGTCCTGCAATACCAGAGTATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGga                                                                                                 <  1:873071/143‑1 (MQ=255)
                                              aCCTGACGGACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGGGTTGTTCCGGGTACACCCGAATCATCAGCTGCCCAGCGgcgtgggg                                                                                                >  2:102880/1‑141 (MQ=255)
                                               ccTGACGAACAGCGTCCTGGTTTTACTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGgaga                                                                                               <  2:489126/143‑1 (MQ=255)
                                                   aCGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGAcgcg                                                                                           >  1:452509/1‑143 (MQ=255)
                                                          gCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGc                                                                                    <  2:452509/143‑1 (MQ=255)
                                                                     ttCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGGGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCGGGAGt                                                                         >  1:317078/1‑143 (MQ=255)
                                                                         tGCAATAACAGAGGATCGTTACCAAGTAAAAGGTACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGGGGAGACCCGTGCAGGCGTTGCCGGAgcaacc                                                                     >  2:500221/1‑138 (MQ=255)
                                                                         tGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGAGCTGGCAGTCATCGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATcc                                                                     <  1:947436/143‑1 (MQ=255)
                                                                            aaTAACGGGGGATCGTTACCAGGTAAGGGGTACGCCGCGCGCAGCCAGTCTGTGAGCTGCGCGCGGAGTGGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGtt                                                                  <  2:154751/143‑1 (MQ=255)
                                                                            aaTAACGGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGTGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGtt                                                                  <  2:258306/143‑1 (MQ=255)
                                                                            aaTAACACAGGATCGTTACCTAGTAAAAGATACGCCGCTCGCATCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGGGACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGtt                                                                  <  1:12699/143‑1 (MQ=255)
                                                                             aTAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTtg                                                                 <  1:502724/143‑1 (MQ=255)
                                                                               caCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGTGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTtgtg                                                               <  1:515765/142‑1 (MQ=255)
                                                                                aCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTGCCCAGCGTGGTGGAGACGCGTGCAGGCGCTGCCGGAGGATCCGTTGTgg                                                              >  1:613622/1‑143 (MQ=255)
                                                                                 cAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGCGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGt                                                             >  2:303783/1‑143 (MQ=255)
                                                                                 cAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTGGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGt                                                             <  1:97574/143‑1 (MQ=255)
                                                                                  cggggATCGTTACCAAGTAAAAGGTACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTc                                                            <  2:798993/140‑1 (MQ=255)
                                                                                      gATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACg                                                        >  1:139127/1‑143 (MQ=255)
                                                                                       aTCGTTATCAAGGAAAAGAAACGGCGCGCGCAGCCCTTCATTGAGCGGCGCGCGGAGTTGTTCCGGGTACAACCGAATCGTCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGt                                                       <  1:402394/143‑1 (MQ=255)
                                                                                             gccAAGTCATAGAGACGCGGCGCGCAGCCCTGCATGGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttc                                                 <  2:751537/142‑1 (MQ=255)
                                                                                             acgaaatAAAATCTACGCCGTGCGCAGCCCTTCATTGAGCTGCGCGCGCAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttc                                                 <  2:547385/137‑1 (MQ=255)
                                                                                             aCCGAGTAAAAGCTACGCCGCGCGCAGCCCCTCATTGAGCTGCGCGCGGCGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGTCGCGCTCTGCGGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttc                                                 <  2:886737/143‑1 (MQ=255)
                                                                                               cAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGGACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTCTCCGTTGTGGTCGGCGTCTGttctt                                               >  1:478095/1‑143 (MQ=255)
                                                                                                  ggaaaaGAGACGCCGCGCGCAGCCCTGCATTGAGCTGCGCGCGGAGTTGTGCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgt                                            <  1:229391/141‑1 (MQ=255)
                                                                                                   tgactgATACGCCGCGCGCTGGGCTTCATTGGGCTGCGCGCGGAGTTGTTCCGGGTACAACCGCATCATCAGTTACCGCGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgtc                                           <  2:574224/138‑1 (MQ=255)
                                                                                                    aaaaGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgtcg                                          <  1:575259/143‑1 (MQ=255)
                                                                                                       cgacaCGCCGCGCGCACCGCTTCAGTGAGCGGCGCGCGGAGTTGTTCCGGCTACAACCGGATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGaa                                       <  1:165925/139‑1 (MQ=255)
                                                                                                        gctACGCCGCGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGTTCCGCGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAAc                                      <  2:608264/141‑1 (MQ=255)
                                                                                                                   cgcAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCgca                           <  1:608144/143‑1 (MQ=255)
                                                                                                                      aGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACCACCGAATCATCAGTTACCCAGCGTGGAGGATACGCGGGGAGGCGCTCGCGGCGGAGCCGGTGCGGCCGACGTGTGGTCGGCGTCTGGACGACTATCGggagca                        >  2:85479/1‑143 (MQ=255)
                                                                                                                        cccTTCATTGAGCTGCGCGCGGAATTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCCGGCGTTGCCGGAGTCCCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGACTATCCGCAGCACg                      >  2:308925/1‑143 (MQ=255)
                                                                                                                         ccTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACgg                     <  2:129452/143‑1 (MQ=255)
                                                                                                                                  aGCTGCGCGGGGGGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGGTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTGGGTCGGAATATCCGCAGCACGGATGCTTGGc            <  2:582325/143‑1 (MQ=255)
                                                                                                                                  aGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGGGGTGGAGACGCGTGCAGGCGTTGCCGGAGTACCCGTTGTGGTCGCCGTCGGTTCTTCGTCGGAACGTATATCCGCAGCACGGATGCTTGGc            >  2:605983/1‑143 (MQ=255)
                                                                                                                                       cgcgcgGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGACTATCCGCAGCACGGATGCTTGCCAGccc       >  2:597585/1‑141 (MQ=255)
                                                                                                                                        gcgcgTAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGGCAGCTTa      <  2:76690/143‑1 (MQ=255)
                                                                                                                                        gcgcgGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGCGTATCCGTTGTGGTCGACGTCTGTTCTTCGGCGGAACGCATCACCGGAGCCCGCATGCTTGGCAGCtgg      >  2:712748/1‑141 (MQ=255)
                                                                                                                                            ggAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTCTCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGGCAGCTTACGaa  >  1:151152/1‑143 (MQ=255)
                                                                                                                                            ggAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGGCAGCTTACgca  >  1:953993/1‑141 (MQ=255)
                                                                                                                                              |                                                                                                                                            
ATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGGCAGCTTACGAA  >  NC_000913/671438‑671720

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.