New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1619859 =	33 (0.810)	15 (0.370)	15/298	1.1	30.9%	coding (286/384 nt)	marA	multiple antibiotic resistance transcriptional regulator
	?	NC_000913	= 1623531	34 (0.830)					coding (320/891 nt)	dgcZ	diguanylate cyclase, zinc‑sensing

GCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCC  >  NC_000913/1623385‑1623681                                                                                                                                                   |                                                                                                                                                       gCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGaaa                                                                                                                                                    >  1:206152/1‑151 (MQ=255)  cTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGaaa                                                                                                                                                    <  2:258511/150‑1 (MQ=255)  cTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGaaa                                                                                                                                                    <  2:299677/150‑1 (MQ=255)      gCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaat                                                                                                                                               <  1:231855‑M2/151‑10 (MQ=255)         ttGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                                                                                                                                            >  2:292855/1‑151 (MQ=255)         ttGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                                                                                                                                            >  1:120028/1‑151 (MQ=255)         ttGCGTAACTGATGATCAAAGGATTCATCAAGAACCCAACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                                                                                                                                            >  1:513841/1‑151 (MQ=255)          tGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaa                                                                                                                                           <  2:582537/151‑1 (MQ=255)           gCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactt                                                                                                                                          <  1:176271‑M2/151‑15 (MQ=255)                aCTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcg                                                                                                                                      >  2:232170/1‑150 (MQ=255)                      gatCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccg                                                                                                                                <  1:303625‑M2/150‑25 (MQ=255)                       atCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTg                                                                                                                               >  1:518808/1‑150 (MQ=255)                        tCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGc                                                                                                                              >  2:629800/1‑150 (MQ=255)                                 tcatcaAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaatac                                                                                                                    >  2:580640‑M2/1‑115 (MQ=255)                                  catcaAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGAtttt                                                                                                                    >  2:221235/1‑150 (MQ=255)                                  catcaAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTc                                                                                                                   <  1:36571/151‑1 (MQ=255)                                   atcaAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccg                                                                                                                  <  1:246341‑M2/151‑39 (MQ=255)                                   atcaAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTca                                                                                                                  >  2:96469/1‑151 (MQ=255)                                    tcaAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTcaa                                                                                                                 >  1:197324/1‑151 (MQ=255)                                     caAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccgga                                                                                                                >  1:176497‑M2/1‑111 (MQ=255)                                          aCCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAg                                                                                                           <  1:138781/151‑1 (MQ=255)                                           cccGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgt                                                                                                                                    <  1:120034‑M2/125‑21 (MQ=255)                                           cccGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgt                                                                                                                                    >  2:120031‑M2/1‑105 (MQ=255)                                                  acCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgc                                                                                                    >  2:303620‑M2/1‑98 (MQ=255)                                                                   aaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGtt                                                                                  >  2:348844/1‑151 (MQ=255)                                                                     aCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgcagggcgaatcgcgcttttta                                                                                <  1:521993‑M2/151‑73 (MQ=255)                                                                       aTCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGc                                                                              <  2:618122/151‑1 (MQ=255)                                                                         ccATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCat                                                                            <  2:475727/151‑1 (MQ=255)                                                                         ccATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCat                                                                            >  2:431537/1‑151 (MQ=255)                                                                           atatTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgcagggcgaatcgcgctttttacatcca                                                                          <  1:400331‑M2/151‑79 (MQ=255)                                                                             atTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGt                                                                        <  1:390899/151‑1 (MQ=255)                                                                              ttGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGtt                                                                       >  2:78459/1‑151 (MQ=255)                                                                              ttGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGtt                                                                       <  2:579759/151‑1 (MQ=255)                                                                                      gATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgcagggcgaatcgcgctttttacatccattaaatcatta                                                               <  2:552312‑M2/151‑90 (MQ=255)                                                                                               cAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATcc                                                         <  1:246515/148‑1 (MQ=255)                                                                                               cAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCAtt                                                      <  2:118923/151‑1 (MQ=255)                                                                                               cAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCAtt                                                      <  2:460583/151‑1 (MQ=255)                                                                                               cAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCAtt                                                      <  2:522428/151‑1 (MQ=255)                                                                                                   taaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGcc                                                  <  2:96080/151‑1 (MQ=255)                                                                                                    aaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCg                                                 <  1:414774/151‑1 (MQ=255)                                                                                                                cGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgcagggcgaatcgcgctttttacatccattaaatcattacaacagctagttgaaaacgtgacaa                                      <  2:199730‑M2/150‑115 (MQ=255)                                                                                                                      aTGCCGCAGTAAAAGAAAGCAACCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgcagggcgaatcgcgctttttacatccattaaatcattacaacagctagttgaaaacgtgacaacgtca                                 <  1:365628‑M2/149‑120 (MQ=255)                                                                                                                            cAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTAt                          >  2:254368/1‑150 (MQ=255)                                                                                                                            cAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTAt                          >  2:257366/1‑150 (MQ=255)                                                                                                                               tAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAg                      <  2:603535/151‑1 (MQ=255)                                                                                                                                           aCCCCTCCTttcaaaaattactttgatgttccgccgcataaataccggatgaccaatatgcagggcgaatcgcgctttttacatccattaaatcattacaacagctagttgaaaacgtgacaacgtcactgaggcaatcatgaaaccactt          >  2:176268‑M2/1‑9 (MQ=255)                                                                                                                                              cctcctGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATc       <  1:78460/151‑1 (MQ=255)                                                                                                                                               ctcctGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCa      >  2:409094/1‑151 (MQ=255)                                                                                                                                                   tGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATcc  <  1:348849/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGATTTTCAACAATAGCCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCC  >  NC_000913/1623385‑1623681

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