New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5618172	42 (0.760)	4 (0.070)	4/248	NT	8.3%	coding (480/1176 nt)	wzy	O‑antigen polymerase
	?	NC_002947	= 5618204	48 (0.890)					coding (448/1176 nt)	wzy	O‑antigen polymerase

GCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAAC  >  NC_002947/5618059‑5618352                                                                                                                                                  |                                                                                                                                                     gCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcct                                                                                                                                                   <  7:82732‑M2/149‑4 (MQ=255)         ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCa                                                                                                                                           <  6:91301/149‑1 (MQ=255)                    ccAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGc                                                                                                                                >  4:137821/1‑149 (MQ=255)                             ccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                                                                                                                       <  3:16920/149‑1 (MQ=255)                                   acaacaACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAAc                                                                                                                 >  4:337003/1‑149 (MQ=255)                                   acaacaACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAAc                                                                                                                 >  4:203655/1‑149 (MQ=255)                                    caacaaCTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGAAATAAGCCCCCAGGATGGGGTGGGATATCTCGCCCATACACGCCAGCCTGAACAGCAAACTCGCACGCTGCTCACCATAGAACt                                                                                                                >  1:303638/1‑149 (MQ=255)                                       caacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtga                                                                                                             >  3:136777/1‑149 (MQ=255)                                        aacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAgg                                                                                                                                       >  1:156426/1‑122 (MQ=255)                                        aacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAgg                                                                                                                                       <  2:156426/122‑1 (MQ=255)                                             cAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGatgatg                                                                                                        <  8:82732/148‑1 (MQ=255)                                              aaGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGga                                                                                                      <  3:203655/149‑1 (MQ=255)                                                          cGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGc                                                                                          >  3:235355/1‑149 (MQ=255)                                                                   tCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCcag                                                                                 >  5:290056/1‑149 (MQ=255)                                                                   tCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCACGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCcag                                                                                 >  8:325585/1‑149 (MQ=255)                                                                      tAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCcagcag                                                                              >  2:67836/1‑149 (MQ=255)                                                                             tCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTg                                                                       <  4:235355/149‑1 (MQ=255)                                                                              cAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGc                                                                      >  1:318839/1‑149 (MQ=255)                                                                                agagCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGcc                                                                     >  8:63013/1‑148 (MQ=255)                                                                                 gagCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCaa                                                                   <  1:151535/149‑1 (MQ=255)                                                                                   gCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAagc                                                                 <  3:337003/149‑1 (MQ=255)                                                                                    cAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAagca                                                                <  1:67836/149‑1 (MQ=255)                                                                                    cAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAag                                                                  >  6:29293/1‑147 (MQ=255)                                                                                    cAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCCCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAagca                                                                >  5:96660/1‑149 (MQ=255)                                                                                       ggCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAgcagca                                                             <  3:137821/149‑1 (MQ=255)                                                                                          cGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAAt                                                           >  4:215042/1‑148 (MQ=255)                                                                                            aaCCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAgcgcg                                                                          >  3:352055/1‑131 (MQ=255)                                                                                            aaCCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAgcgcg                                                                          <  4:352055/131‑1 (MQ=255)                                                                                                aaTGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGAc                                                    >  3:5325/1‑149 (MQ=255)                                                                                                aaTGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGAc                                                    >  5:188475/1‑149 (MQ=255)                                                                                                        aaGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAAt                                                           >  5:181143/1‑134 (MQ=255)                                                                                                        aaGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAAt                                                           <  6:181143/134‑1 (MQ=255)                                                                                                            ccccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcctgggggcttatgtcattggttcggccttgctcttgatgctctacgagccacctcggcgacgtggcttgcagttgttgtggctggcggcactggctt                                                    <  4:5589‑M2/137‑99 (MQ=255)                                                                                                            ccccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcctgggggcttatgtcattggttcggccttgctcttgatgctctacgagccacctcggcgacgtggcttgcagttgttgtggctggcggcactggctt                                                    >  3:5589‑M2/1‑39 (MQ=255)                                                                                                            ccccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCtt                                        >  3:23331/1‑149 (MQ=255)                                                                                                              ccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcctgggggcttatgtcattggttcggccttgctcttgatgctctacgagccacctcggcgacgtggcttgcagttgttgtggctggcggcactgg                                                       >  5:245635‑M2/1‑37 (MQ=255)                                                                                                              ccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcctgggggcttatgtcattggttcggccttgctcttgatgctctacgagccacctcggcgacgtggcttgcagttgttgtggctggcggcactgg                                                       <  6:245635‑M2/132‑96 (MQ=255)                                                                                                              ccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTgg                                      <  3:215042/149‑1 (MQ=255)                                                                                                                  gATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCtt                                        <  5:220178/143‑1 (MQ=255)                                                                                                                  gATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCtt                                        >  6:220178/1‑143 (MQ=255)                                                                                                                   aTGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAg                                 <  7:325585/149‑1 (MQ=255)                                                                                                                        gTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCTCACGATGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTTGCCTAGCTggg                            <  6:96660/149‑1 (MQ=255)                                                                                                                            gATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGAACAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCaa                        >  5:13189/1‑149 (MQ=255)                                                                                                                                  tcGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCgg                                            <  8:233762/123‑1 (MQ=255)                                                                                                                                  tcGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCgg                                            >  7:233762/1‑123 (MQ=255)                                                                                                                                  tcGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTg                              <  2:289298/137‑1 (MQ=255)                                                                                                                                  tcGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCg                   >  6:110129/1‑148 (MQ=255)                                                                                                                                  tcGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCATGGCCAAGCTg                              >  1:289298/1‑137 (MQ=255)                                                                                                                                      cGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAAc              <  2:285782/149‑1 (MQ=255)                                                                                                                                          aCCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCca          <  2:318839/149‑1 (MQ=255)                                                                                                                                           cccGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCca          <  7:63013/148‑1 (MQ=255)                                                                                                                                           cccGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCca          <  4:136777/148‑1 (MQ=255)                                                                                                                                            ccGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCcag         >  2:314812/1‑148 (MQ=255)                                                                                                                                                  gCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAac  <  1:314812/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAAC  >  NC_002947/5618059‑5618352

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