New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3883405	38 (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.

CTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAT  >  NC_000913/3883262‑3883551                                                                                                                                                |                                                                                                                                                   cTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGtt                                                                                                                                               <  6:45213/149‑1 (MQ=255) cTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGGCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGtt                                                                                                                                               <  6:42936/149‑1 (MQ=255)  tGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTg                                                                                                                                              >  5:117730/1‑149 (MQ=255)         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                                                                                                                          >  4:91413/1‑118 (MQ=255)                                                     tCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCt                                                                                                                          <  3:91413/118‑1 (MQ=255)                                                             gtgtTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGc                                                                                                                                       <  3:15984/97‑1 (MQ=255)                                                             gtgtTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGc                                                                                                                                       >  4:15984/1‑97 (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                                                                                                                   >  4:129286/1‑102 (MQ=255)                                                                            gTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                   <  3:129286/102‑1 (MQ=255)                                                                             ttAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCgg                                                                   >  4:102616/1‑149 (MQ=255)                                                                                gAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGtt                                                                 <  5:92654/148‑1 (MQ=255)                                                                                             gTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                   <  3:26038/85‑1 (MQ=255)                                                                                             gTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAgg                                                                                                                   >  4:26038/1‑85 (MQ=255)                                                                                                   ccggggccgggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCaaaa                                                <  6:106332/146‑1 (MQ=255)                                                                                                    cggggccgggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGTCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCa                                            <  6:23936/149‑1 (MQ=255)                                                                                                           gggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCgg                                                                   >  4:68297/1‑119 (MQ=255)                                                                                                           gggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCgg                                                                   <  3:68297/119‑1 (MQ=255)                                                                                                           gggACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTgg                                     <  3:102616/149‑1 (MQ=255)                                                                                                                      ccAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgctc                          >  2:141504/1‑149 (MQ=255)                                                                                                                       cAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCAtccgctcc                         >  3:51525/1‑149 (MQ=255)                                                                                                                           cTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTggg                                    >  3:141599/1‑134 (MQ=255)                                                                                                                           cTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTggg                                    <  4:141599/134‑1 (MQ=255)                                                                                                                               gcgcGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAg                                           <  2:124606/123‑1 (MQ=255)                                                                                                                               gcgcGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAg                                           >  1:124606/1‑123 (MQ=255)                                                                                                                                   gTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACtt                                                     <  7:17055/109‑1 (MQ=255)                                                                                                                                   gTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACtt                                                     >  8:17055/1‑109 (MQ=255)                                                                                                                                              agcaCGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAt  <  1:52984/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCGCAGCACGTTGCGGCTGCGTTTGCGCCACCTGTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCAT  >  NC_000913/3883262‑3883551

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