New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4608119 =	27 (0.780)	4 (0.120)	4/258	4.0	12.7%	coding (317/414 nt)	holD	DNA polymerase III, psi subunit
	?	NC_000913	4608161 =	29 (0.870)					coding (359/414 nt)	holD	DNA polymerase III, psi subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGA  >  NC_000913/4607972‑4608266                                                                                                                                                    |                                                                                                                                                    ttCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATGACCg                                                                                                                                                    >  5:131108/1‑149 (MQ=255)    tGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGgcg                                                                                                                                                 >  5:135115/1‑149 (MQ=255)                         ccGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAAc                                                                                                                            <  6:131108/149‑1 (MQ=255)                               aGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACg                                                                                                                      <  2:164554/149‑1 (MQ=255)                                           tGACGCCAGAAAAAATCGCGATGCTGCCACAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCcgcg                                                                                                          >  4:39230/1‑149 (MQ=255)                                                ccAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATg                                                                                                     <  3:46545/149‑1 (MQ=255)                                                    aaaaaaTCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGcaa                                                                                                 <  1:13725/149‑1 (MQ=255)                                                            gcgATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCGCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAAttt                                                                                          <  8:84979/148‑1 (MQ=255)                                                                       gCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAac                                                                              <  6:135115/149‑1 (MQ=255)                                                                                  aCTGCAACAGTTGGCGTTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTAACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGATATGGCAACAAATTTGAACATATGAACACGCTttcttc                                                                   >  4:118568/1‑149 (MQ=255)                                                                                       aaCAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTc                                                               >  4:106773/1‑148 (MQ=255)                                                                                       aaCAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCg                                                              <  3:106773/149‑1 (MQ=255)                                                                                                      ggTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATgg                                                                                                    >  4:140442/1‑96 (MQ=255)                                                                                                      ggTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATgg                                                                                                    <  3:140442/96‑1 (MQ=255)                                                                                                                        tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATttctt                                                                    <  6:51895/110‑1 (MQ=255)                                                                                                                        tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATttctt                                                                    >  5:51895/1‑110 (MQ=255)                                                                                                                        tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          <  5:112675/120‑1 (MQ=255)                                                                                                                        tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          >  6:112675/1‑120 (MQ=255)                                                                                                                          aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          <  4:26563/118‑1 (MQ=255)                                                                                                                          aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          >  3:26563/1‑118 (MQ=255)                                                                                                                          aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAt                               <  5:39317/145‑1 (MQ=255)                                                                                                                          aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAt                               >  6:39317/1‑145 (MQ=255)                                                                                                                            tGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAtt                              >  1:74717/1‑144 (MQ=255)                                                                                                                            tGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAtt                              <  2:74717/144‑1 (MQ=255)                                                                                                                                    gcTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGccc                 >  5:145830/1‑149 (MQ=255)                                                                                                                                                   ccGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACTACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGa  <  3:39230/149‑1 (MQ=255)                                                                                                                                                   ccGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGa  <  6:145830/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGA  >  NC_000913/4607972‑4608266

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.