New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 388340538 (1.110)6 (0.180) 6/242 3.2 14.6% coding (325/1404 nt) dnaA chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
?NC_000913 = 3883413 34 (1.050)coding (317/1404 nt) dnaA chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTAT  >  NC_000913/3883270‑3883560
                                                                                                                                               |                                                                                                                                                   
ccgccgCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGc                                                                                                                                                <  3:143360/149‑1 (MQ=255)
           gcCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCAc                                                                                                                                     <  6:117730/149‑1 (MQ=255)
                              cTTCAACGAAGTTATCAAAAGTGTTTTTTACGTTTACGTTAGAAAGATAGGTCGGTTCTGACGGGGCCCGGACGATATCCCAACCTGAGCGCGAAGAAGGCGAAGCACGTTGCGGCTGCGTTTCAGCAACCTCTGAAGTTCCCGCGAc                                                                                                                   >  2:52984/1‑148 (MQ=255)
                                 cAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGtt                                                                                                                <  2:17513/148‑1 (MQ=255)
                                  aaCGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGc                                                                                                                               <  6:64252/132‑1 (MQ=255)
                                  aaCGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGc                                                                                                                               >  5:64252/1‑132 (MQ=255)
                                            tCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCt                                                                                                                                   <  3:91413/118‑1 (MQ=255)
                                            tCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCt                                                                                                                                   >  4:91413/1‑118 (MQ=255)
                                                    gtgtTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGc                                                                                                                                                >  4:15984/1‑97 (MQ=255)
                                                    gtgtTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGc                                                                                                                                                <  3:15984/97‑1 (MQ=255)
                                                     tgtTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTgg                                                                                           >  4:96984/1‑149 (MQ=255)
                                                       tttGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGTGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCg                                                                                         >  8:47748/1‑149 (MQ=255)
                                                           aCGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTg                                                                                     >  6:92654/1‑149 (MQ=255)
                                                              tttACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGt                                                                                  <  3:96984/149‑1 (MQ=255)
                                                                  cGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCAcc                                                                              >  5:106332/1‑149 (MQ=255)
                                                                   gTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                            <  3:129286/102‑1 (MQ=255)
                                                                   gTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                            >  4:129286/1‑102 (MQ=255)
                                                                    ttAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCgg                                                                            >  4:102616/1‑149 (MQ=255)
                                                                       gAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGtt                                                                          <  5:92654/148‑1 (MQ=255)
                                                                                    gTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                            >  4:26038/1‑85 (MQ=255)
                                                                                    gTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                            <  3:26038/85‑1 (MQ=255)
                                                                                          ccggggccgggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCaaaa                                                         <  6:106332/146‑1 (MQ=255)
                                                                                           cggggccgggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGTCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCa                                                     <  6:23936/149‑1 (MQ=255)
                                                                                                  gggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCgg                                                                            <  3:68297/119‑1 (MQ=255)
                                                                                                  gggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCgg                                                                            >  4:68297/1‑119 (MQ=255)
                                                                                                  gggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTgg                                              <  3:102616/149‑1 (MQ=255)
                                                                                                             ccAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgctc                                   >  2:141504/1‑149 (MQ=255)
                                                                                                              cAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgctcc                                  >  3:51525/1‑149 (MQ=255)
                                                                                                                  cTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTggg                                             <  4:141599/134‑1 (MQ=255)
                                                                                                                  cTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTggg                                             >  3:141599/1‑134 (MQ=255)
                                                                                                                      gcgcGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAg                                                    >  1:124606/1‑123 (MQ=255)
                                                                                                                      gcgcGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAg                                                    <  2:124606/123‑1 (MQ=255)
                                                                                                                          gTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACtt                                                              <  7:17055/109‑1 (MQ=255)
                                                                                                                          gTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACtt                                                              >  8:17055/1‑109 (MQ=255)
                                                                                                                                     agcaCGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAt           <  1:52984/149‑1 (MQ=255)
                                                                                                                                          gTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGata      <  1:141504/149‑1 (MQ=255)
                                                                                                                                               ggCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGAtattat  >  2:159893/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTAT  >  NC_000913/3883270‑3883560

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