New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 670713 =34 (0.760)6 (0.150) 5/276 4.0 15.9% coding (893/1032 nt) holA DNA polymerase III, delta subunit
?NC_000913 670903 = 32 (0.770)coding (703/1032 nt) holA DNA polymerase III, delta subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAGCCAGCGCCAGCAGGTTACCTTCATAACAGTAGCAGAGCACCTGATTTGCCGCGTCATCCAGTTCTAAGTTGA  >  NC_000913/670687‑671133
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
gtgtCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATgcgcgc                                                                                                                                                                                                                  <  1:124007/239‑1 (MQ=255)
                   cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATcaaca                                                                                                                                                                                               <  1:449283/239‑1 (MQ=255)
                                                                     cccGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCt                                                                                                                                                                                                                                      >  1:609747/1‑150 (MQ=255)
                                                                     cccGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCt                                                                                                                                                                                                                                      <  2:609747/150‑1 (MQ=255)
                                                                                    aTACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAAAGCGCGCTTACTTTTTTCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTTACCGCCTGTTCAACGCCCggg                                                                                                                              >  1:717996/1‑238 (MQ=255)
                                                                                                  tCAAACAACGCACGCAGTGGCGTATGGGCAGTCTGGCGTTTCAGGTTAACCAGTAACAACAGTTCAGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATg                                                                                                                                                                                                                       <  1:544350/135‑1 (MQ=255)
                                                                                                  tCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCAGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATg                                                                                                                                                                                                                       >  2:544350/1‑135 (MQ=255)
                                                                                                          cgcacgcaGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCt                                                                                                                                                                                                                                      >  2:465225/1‑113 (MQ=255)
                                                                                                          cgcacgcaGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCt                                                                                                                                                                                                                                      <  1:465225/113‑1 (MQ=255)
                                                                                                          cgcacgcaGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGcat                                                                                                                                                        >  2:235903/1‑191 (MQ=255)
                                                                                                          cgcacgcaGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGcat                                                                                                                                                        <  1:235903/191‑1 (MQ=255)
                                                                                                               gcaGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAAcg                                                                                                                                    >  1:387776/1‑206 (MQ=255)
                                                                                                               gcaGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAAcg                                                                                                                                    <  2:387776/206‑1 (MQ=255)
                                                                                                                  agtcagtgtGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGagcagc                                                                                                >  2:312809/9‑239 (MQ=255)
                                                                                                                    ggCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTga                                                                                                                                                                                                                                    >  1:195605/1‑105 (MQ=255)
                                                                                                                    ggCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTga                                                                                                                                                                                                                                    <  2:195605/105‑1 (MQ=255)
                                                                                                                           gggCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCt                                                                                                                                                                                                                                      <  1:302479/96‑1 (MQ=255)
                                                                                                                           gggCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCt                                                                                                                                                                                                                                      >  2:302479/1‑96 (MQ=255)
                                                                                                                             gCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATcaaca                                                                                                                                                                                               <  2:578954/133‑1 (MQ=255)
                                                                                                                             gCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATcaaca                                                                                                                                                                                               >  1:578954/1‑133 (MQ=255)
                                                                                                                             gCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCt                                                                                      >  2:281581/1‑238 (MQ=255)
                                                                                                                             gCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCAATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGCAAATTTGCCGTCTGGCCAGAGCAACGATAAACGCt                                                                                      >  2:277203/1‑238 (MQ=255)
                                                                                                                                      cGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGACTATGCAATGCGCGCTTCCTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCCCGGCATCATTCACCGCCTGTTCACCGAGCGGTAATGTCAATTTCACGTCTGGCCCAGGCATCGGTAAACCGTCCAGTGccgt                                                                            >  2:17688/1‑237 (MQ=255)
                                                                                                                                               ttAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAgcca                                                                       <  1:312809/235‑1 (MQ=255)
                                                                                                                                                                ttCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCa                                                                                                                                                                                  <  1:387844/111‑1 (MQ=255)
                                                                                                                                                                ttCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCa                                                                                                                                                                                  >  2:387844/1‑111 (MQ=255)
                                                                                                                                                                 tCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCTATAAAAGCACCCGAGCCCGAGTCACCACCTGCCGGttcactt                                                     >  2:398751/1‑230 (MQ=255)
                                                                                                                                                                        gTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAGCCAGCGCCAGCAGGTTACCTTCATAACAGTAg                                          >  2:301667/1‑239 (MQ=255)
                                                                                                                                                                                cgcAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGca                                                                                                                                                         >  2:308363/1‑120 (MQ=255)
                                                                                                                                                                                cgcAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGca                                                                                                                                                         <  1:308363/120‑1 (MQ=255)
                                                                                                                                                                                                  tCGCTCCCTTCCAGACGCATTTGCGGAAAATTATGCAATGCGCGTTTATTTTTCCCCATCACCAAAGCATCAACCAAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAGCCAGCGCCAGCAGGTTACCTTCATAACAGTAGCAGAGCACCTGATTTGCCGCGTCATc                <  2:717996/239‑1 (MQ=255)
                                                                                                                                                                                                         cTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAGCCAGCGCCAGCAGGTTACCTTCATAACAGTAGCAGAGCACCTGATTTGCCGCGTCATCCAGTTCt         <  1:160798/239‑1 (MQ=255)
                                                                                                                                                                                                                  gCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAAc                                                                                                                                                                                     <  2:537700/58‑1 (MQ=255)
                                                                                                                                                                                                                  gCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAAc                                                                                                                                                                                     >  1:537700/1‑58 (MQ=255)
                                                                                                                                                                                                                    aGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATcaacaa                                                                                                                                                                                              <  1:683368/47‑1 (MQ=255)
                                                                                                                                                                                                                    aGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATcaacaa                                                                                                                                                                                              >  2:683368/1‑47 (MQ=255)
                                                                                                                                                                                                                    aGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAGCCAGCGCCAGCAGGTTACCTTCATAACAGTAGCAGAGCACCTGATTTGCCGCGTCATCCAGTTCTAAGTTGa  >  1:636014/1‑235 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
GTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTTCGCTGCCTTCCAGACGCAGTTGCTGAAGAATATGCAATGCGCGCTTACTTTTTCCCATCAACAAAGCATCAACCCAATGAAAAGGGGTGAAATGCGCGGCATCATTCACCGCCTGTTCAACGCGCGGTAATGTCAATTTGCCGTCTGGCCAGAGCAGCGATAAACGCTCCAGTGCCTGAGCCAGCGCCAGCAGGTTACCTTCATAACAGTAGCAGAGCACCTGATTTGCCGCGTCATCCAGTTCTAAGTTGA  >  NC_000913/670687‑671133

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