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.

TTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGA  >  NC_000913/4607972‑4608266
                                                                                                                                                   |                                                                                                                                                   
ttCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATGACCg                                                                                                                                                    >  3:131108/1‑149 (MQ=255)
   tGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGgcg                                                                                                                                                 >  3:135115/1‑149 (MQ=255)
                        ccGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAAc                                                                                                                            <  4:131108/149‑1 (MQ=255)
                           aaacgttttttcactccacggaaaggcgtgggcgcgttgttcaatgtggtaagccgccggtaaatcagtcgtttcgagggaagaaatcgtgttcatatgtgcaaatttgttgccataacgCGGCGCTCATCGATTTACGGGCAAAccc                                                                                                                          <  5:91754‑M1/28‑1 (MQ=255)
                                       actccacggaaaggcgtgggcgcgttgttcaatgtggtaagccgccggtaaatcagtcgtttcgagggaagaaatcgtgttcatatgtgcaaatttgttgccataacgCGGCGCTCATCGATTTACGGGCAAACCCAACGGCACGcgcc                                                                                                             >  6:91754‑M1/109‑149 (MQ=255)
                                          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                                                                                                           <  3:107416‑M1/43‑1 (MQ=255)
                                                                                  gccggtaaatcagtcgtttcgagggaagaaatcgtgttcatatgtgcaaatttgttgccataacgCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGcgccgc                                                                                                           >  4:107416‑M1/66‑108 (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                                                                    >  3:51895/1‑110 (MQ=255)
                                                                                                                       tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATttctt                                                                    <  4:51895/110‑1 (MQ=255)
                                                                                                                       tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          >  4:112675/1‑120 (MQ=255)
                                                                                                                       tCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          <  3:112675/120‑1 (MQ=255)
                                                                                                                         aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          >  1:26563/1‑118 (MQ=255)
                                                                                                                         aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAAc                                                          <  2:26563/118‑1 (MQ=255)
                                                                                                                         aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAt                               <  3:39317/145‑1 (MQ=255)
                                                                                                                         aCTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACAt                               >  4:39317/1‑145 (MQ=255)
                                                                                                                                   gcTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGccc                 >  3:145830/1‑149 (MQ=255)
                                                                                                                                                  ccGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACTACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGa  <  1:39230/149‑1 (MQ=255)
                                                                                                                                                  ccGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGa  <  4: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.