New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4608119 =	23 (0.930)	4 (0.170)	4/258	2.7	14.5%	coding (317/414 nt)	holD	DNA polymerase III, psi subunit
	?	NC_000913	4608161 =	25 (1.050)					coding (359/414 nt)	holD	DNA polymerase III, psi subunit
Rejected: Frequency below/above cutoff threshold.

TGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACGTTTGCCAGC  >  NC_000913/4608014‑4608286                                                                                                                                                    |                                                                                                                              tGACGCCAGAAAAAATCGCGATGCTGCCACAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCcgcg                                                                                                                              >  2:39230/1‑149 (MQ=255)      ccAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATg                                                                                                                         <  1:46545/149‑1 (MQ=255)                  gcgATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCGCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAAttt                                                                                                              <  6:84979/148‑1 (MQ=255)                             gCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAac                                                                                                  <  4:135115/149‑1 (MQ=255)                                        aCTGCAACAGTTGGCGTTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTAACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGATATGGCAACAAATTTGAACATATGAACACGCTttcttc                                                                                       >  2:118568/1‑149 (MQ=255)                                         gccggtaaatcagtcgtttcgagggaagaaatcgtgttcatatgtgcaaatttgttgccataacgCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGcgccgc                                                                                                                               >  4:107416‑M1/66‑108 (MQ=255)                                         gccggtaaatcagtcgtttcgagggaagaaatcgtgttcatatgtgcaaatttgttgccataacgCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGcgccgc                                                                                                                               <  3:107416‑M1/43‑1 (MQ=255)                                             aaCAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTc                                                                                   >  2:106773/1‑148 (MQ=255)                                             aaCAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCg                                                                                  <  1:106773/149‑1 (MQ=255)                                                            ggTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATgg                                                                                                                        <  1:140442/96‑1 (MQ=255)                                                            ggTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATgg                                                                                                                        >  2:140442/1‑96 (MQ=255)                                                                              tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATttctt                                                                                        <  4:51895/110‑1 (MQ=255)                                                                              tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATttctt                                                                                        >  3:51895/1‑110 (MQ=255)                                                                              tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                                              <  3:112675/120‑1 (MQ=255)                                                                              tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                                              >  4:112675/1‑120 (MQ=255)                                                                                aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                                              <  2:26563/118‑1 (MQ=255)                                                                                aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                                              >  1:26563/1‑118 (MQ=255)                                                                                aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAt                                                   <  3:39317/145‑1 (MQ=255)                                                                                aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAt                                                   >  4:39317/1‑145 (MQ=255)                                                                                          gcTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGccc                                     >  3:145830/1‑149 (MQ=255)                                                                                                         gcggcgcgtgccgttgggtttgcccgtaaatcggtgagcgccgCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACcggc                                                               <  3:107416‑M2/65‑1 (MQ=255)                                                                                                         gcggcgcgtgccgttgggtttgcccgtaaatcggtgagcgccgCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACcggc                                                               >  4:107416‑M2/44‑108 (MQ=255)                                                                                                         ccGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACTACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGa                      <  1:39230/149‑1 (MQ=255)                                                                                                         ccGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGa                      <  4:145830/149‑1 (MQ=255)                                                                                                           ggcgcgtgccgttgggtttgcccgtaaatcgatgagcgccgCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGt                    >  6:91754‑M2/42‑149 (MQ=255)                                                                                                           GGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGt                    <  6:88379‑M1/149‑1 (MQ=255)                                                                                                                tCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGaaa                >  4:46600/1‑148 (MQ=255)                                                                                                                       tttACTGGCAAACCCATCGGCACGCGCCGCGTTATTGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGAGGCTTACCAGATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACGttt        <  1:118568/149‑1 (MQ=255)                                                                                                                        gggtttgcccgtaaatcgatgagcgccgCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACGttt        <  5:91754‑M2/120‑1 (MQ=255)                                                                                                                            gggCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGAt                                                                       >  4:75967/1‑81 (MQ=255)                                                                                                                            gggCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGAt                                                                       <  3:75967/81‑1 (MQ=255)                                                                                                                              gCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACGTTTGCCAGc  >  4:81981/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                              TGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACGTTTGCCAGC  >  NC_000913/4608014‑4608286

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.