New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 2127296 =7 (0.400)4 (0.240) 3/292 NT 31.1% coding (1838/2163 nt) wzc colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
?CP009273 2129302 = 11 (0.660)coding (278/444 nt) wzb colanic acid production protein‑tyrosine‑phosphatase; Wzc‑P dephosphorylase

GGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGG  >  CP009273/2127036‑2127473
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                 
ggAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTGCTGTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACACGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACcgcgc                                                                                                                                                                     <  1:88103/275‑1 (MQ=255)
                                                                                                     gcCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACcgcgcg                                                                                                                                                                    >  2:225469/1‑175 (MQ=255)
                                                                                                         ccaTTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGAt                                                                                                                                                                <  1:225475/175‑1 (MQ=255)
                                                                                                               aCGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTNGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGTCATTCCCGTTATTTGTGCCCAACAGCTCg                                                      <  1:229583/275‑1 (MQ=255)
                                                                                                                                                                    ggcggcGTATCAATCAACACCAGGTCATATTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGTTTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTgg  <  1:162428/274‑1 (MQ=255)
                                                                                                                                                                                ccgtcaaatcagccgccgtctgtgtcgcaactacgacctgattttgaccatggaaaagcgccatatcgaacgcttatgcgagatGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTag                                                                                                                                <  1:288102‑M1/52‑1 (MQ=255)
                                                                                                                                                                                 cgtcaaattagccgccgtctgtgtcgcaactacgacctgattttgaccatggaaaagcgccatatcgaacgcttatgcgagatGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATATAGGTCGGTTTag                                                                                                                                <  1:299355‑M1/52‑1 (MQ=255)
                                                                                                                                                                                    caaatcagccgccgtctgtgtcgcaactacgacctgattttgaccatggaaaagcgccatatcgaacgcttatgcgagatGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTagcagc                                                                                                                            >  2:288091‑M1/81‑136 (MQ=255)
                                                                                                                                                                                     gaattagccgccgtctgagtcgcaactacgtcctgatttttaccatggcaaagggcaatatcgaacgcttatgcgtgatGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTagcagc                                                                                                                            >  2:299344‑M1/80‑135 (MQ=255)
                                                                                                                                                                                          ggTCATAGTTTTTACTCGCCCAGTTCACTAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACg                                                                           >  1:247987/1‑179 (MQ=255)
                                                                                                                                                                                              aTAGTTTTTACTCGCCCAGTTCACCCGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGttat                                                                       <  2:247979/179‑1 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                 
GGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGG  >  CP009273/2127036‑2127473

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ 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.