New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 621688	22 (0.710)	5 (0.170)	3/356	NT	17.0%	coding (502/1005 nt)	fepD	ferric enterobactin ABC transporter membrane subunit FebD
	?	NC_000913	623383 =	28 (0.960)					coding (1084/1251 nt)	entS	enterobactin exporter EntS

AGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAAC  >  NC_000913/623131‑623614                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        aGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAggcg                                                                                                                                                                                                                                   <  1:330829/259‑1 (MQ=255)                            ccAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACAACgg                                                                                                                                                                                            >  1:360114/1‑271 (MQ=255)                                   tATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCgctgct                                                                                                                                                                                                                         <  1:79729/235‑1 (MQ=255)                                                              cGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATGGGCGGGGCGCTGCTGGGTGGTTTGGGCGGGAGGATGACACGGGTTGCTTCCGCCAGCGCGGGCGGCTGTGGTttgttg                                                                                                                                                         >  1:13022/1‑272 (MQ=255)                                                               gTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGa                                                                                                                                                        <  1:157366/272‑1 (MQ=255)                                                                   cTCGCCATTGTTTTTTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTa                                                                                                                                                     <  1:392280/271‑1 (MQ=255)                                                                       ccATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCGGGGGGGTTTGGGCGCGATGGTGACACCGGTTGGTTCCGCCAGCGGGGGGGGGTTTGtgttgtggctttagggg                                                                                                                                                >  2:392280/1‑261 (MQ=255)                                                                                                        agcgcgcgactcagtaatagcgcagttgctccggcgatcagcaccgggatcagcaccacttttaaggtatgtagattgcgaatatccagcgaaccggcttgccagaaacgcaactgatcgtagacgtcaggattaagcagggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGa                                                                                                               >  1:365893‑M2/150‑272 (MQ=255)                                                                                                                                     ggCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGGTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcg                                                                                                                                                                             >  1:112369/1‑181 (MQ=255)                                                                                                                                      gCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcgc                                                                                                                                                                            <  2:112369/181‑1 (MQ=255)                                                                                                                                             gAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGCAGCGATGTTAGGGCGGATTAACGGTTTGTGGGCGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGCTGACACCGGTTGCTTCCGCAAGCGCGGGCGGGTTTGGGTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGGGACATTTTCGTCAGACGCCGCCGCAGGTGAAAGAGTCCg                                                                          >  1:461242/1‑272 (MQ=255)                                                                                                                                                     tCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcgc                                                                                                                                                                            <  1:452482/166‑1 (MQ=255)                                                                                                                                                        gCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGc                                                                                                                                                                        >  1:321949/1‑167 (MQ=255)                                                                                                                                                         cTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCg                                                                                                                                                                       <  2:321949/167‑1 (MQ=255)                                                                                                                                                              accacttttaaggtatgtagattgcgaatatccagcgaaccggcttgccagaaacgcaactgatcgtagacgtcaggattaagcagggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcg                                                                                                                                                                             >  1:60504‑M2/96‑156 (MQ=255)                                                                                                                                                              accacttttaaggtatgtagattgcgaatatccagcgaaccggcttgccagaaacgcaactgatcgtagacgtcaggattaagcagggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACGGGGTCCGACAGTTAATGCTTaaa                                                          >  1:255839‑M2/96‑271 (MQ=255)                                                                                                                                                              accacttttaaggtatgtagattgcgaatatccagcgaaccggcttgccagaaacgcaactgatcgtagacgtcaggattaagcagggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGATTaa                                                           >  1:5362‑M2/96‑270 (MQ=255)                                                                                                                                                               ccacttttaaggtatgtagattgcgaatatccagcgaaccggcttgccagaaacgcaactgatcgtagacgtcaggattaagcagggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcgc                                                                                                                                                                            <  2:60504‑M2/62‑1 (MQ=255)                                                                                                                                                                    aaTACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAgcg                                                    <  1:444183/271‑1 (MQ=255)                                                                                                                                                                           aaTGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTc                                                                                                                                                                                     >  1:484898/1‑135 (MQ=255)                                                                                                                                                                            aTGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTcc                                                                                                                                                                                    <  2:484898/135‑1 (MQ=255)                                                                                                                                                                                tgcAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATa                                                                                                                                                                                                                                       >  2:397747/1‑80 (MQ=255)                                                                                                                                                                                 gcaaacgcaaacCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAg                                                                                                                                                                                                                                      <  1:397747/80‑1 (MQ=255)                                                                                                                                                                                     acgcaaacCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATcc                                      >  1:64143/1‑268 (MQ=255)                                                                                                                                                                                          aacCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATa                                                                                                                                                                                                                                       >  2:378700/1‑70 (MQ=255)                                                                                                                                                                                           acCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAg                                                                                                                                                                                                                                      <  1:378700/70‑1 (MQ=255)                                                                                                                                                                                              ccGGATGCGTTGTTGGGGCGGTTTAACGGTTGGTGGACGGGGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATg                          <  2:64143/271‑1 (MQ=255)                                                                                                                                                                                                                 ggATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTaaa                                                          >  1:223978/1‑220 (MQ=255)                                                                                                                                                                                                                     tAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcg                                                                                                                                                                             >  1:355900/1‑101 (MQ=255)                                                                                                                                                                                                                     tAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTttgttg                                                                                                                                                         >  2:136402/1‑121 (MQ=255)                                                                                                                                                                                                                     tAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCAGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCAAGCACTTGCATGGCGCTGTAGTACTCCAGACGGAAc  >  1:56182/1‑272 (MQ=255)                                                                                                                                                                                                                      aaCGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAgcgc                                                                                                                                                                            <  2:355900/101‑1 (MQ=255)                                                                                                                                                                                                                      aaCGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGa                                                                                                                                                        <  1:136402/121‑1 (MQ=255)                                                                                                                                                                                                                                ggACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGc                                                                                                                                                >  1:427178/1‑119 (MQ=255)                                                                                                                                                                                                                                 gACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCg                                                                                                                                               <  2:427178/119‑1 (MQ=255)                                                                                                                                                                                                                                       gcAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTaaa                                                          >  1:442613/1‑198 (MQ=255)                                                                                                                                                                                                                                        cAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTaaaa                                                         <  2:442613/198‑1 (MQ=255)                                                                                                                                                                                                                                                   gggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCtt                                                >  2:280433‑M2/11‑196 (MQ=255)                                                                                                                                                                                                                                                    ggcgatgccATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTa                                               <  1:280433‑M2/187‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        AGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAAC  >  NC_000913/623131‑623614

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