New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1774019 =	45 (1.180)	4 (0.110)	4/268	4.7	8.9%	coding (231/867 nt)	ydiB	quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
	?	NC_000913	1774066 =	39 (1.080)					coding (278/867 nt)	ydiB	quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATAGCTTTCCTGGAGCAATTGAAGGATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTG  >  NC_000913/1773927‑1774206                                                                                                                                            |                                                                                                                                             gATAGCTTTCCTGGAACAAATGAAGGATTAAAAGCCCTCAAAATACGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCAt                                                                                                                                           >  2:52272/1‑143 (MQ=255)          ccgtgggCAATTGAAGGACTAAAAGCCCTGAAAACGCGCGGAACTGGTGGATCGATGCCGAACAAACGCCTGGCGCATGAATATGTTGTTGAATTAACACAAGCGGCCAAACTGGTGGGGGCCATCAACACCATCGTTAatga                                                                                                                                  <  1:655058/137‑1 (MQ=255)           cTGGAGCAATTGAAGGATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACCCCAGCTGCCAAGCTGGGGGGGTCCATCCACACCATCGTTaataat                                                                                                                                 >  1:122032/1‑143 (MQ=255)            tGGAGCAATTGAAGGATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACACACAACTGGCGTGTGAATATGTTGCTGACTTAACACCAGCTGCCAAACTGATGGGGACCATCAACAACGTCGTTAatgatg                                                                                                                                >  1:396918/1‑143 (MQ=255)                 cAATTGAAGGATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCCGCTGCCAAGCCGGTGGGGGCCATCAACACCATCGTTAATGATGATGGc                                                                                                                           >  1:227812/1‑143 (MQ=255)                    ttGAAGGATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATAGCTAt                                                                                                                        >  1:203068/1‑143 (MQ=255)                         ggATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCg                                                                                                                   <  2:455678/143‑1 (MQ=255)                         ggATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGCTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCg                                                                                                                   >  1:532348/1‑143 (MQ=255)                           aTTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTg                                                                                                                 <  1:371505/143‑1 (MQ=255)                            ttAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTgg                                                                                                                >  2:578284/1‑143 (MQ=255)                                aaGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCtat                                                                                                            <  1:395351/143‑1 (MQ=255)                                     cTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGGATTATCACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAacac                                                                                                       >  1:518517/1‑143 (MQ=255)                                     cTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGCCCCTCAACACCAACGTTAATGATCATGGCTATCTGCGTGGCCATAacac                                                                                                       >  2:588281/1‑143 (MQ=255)                                      tCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGGGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACAcc                                                                                                      <  2:459812/143‑1 (MQ=255)                                      tCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGGGTGAATATGAGGATGAATTAACACCAGCTGCCAAACTGGGGGGGGCAATCAACACCAGCGATAATGATGATGGCCATCTGCGAGGCGATAACAcc                                                                                                      >  2:72196/1‑143 (MQ=255)                                           aTGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATCATGATGGCTAACCGCGAGGCTAcaccacgcccgg                                                                                                 >  2:600398/1‑137 (MQ=255)                                                gagAACTGGTCTAGCGCTGCCCAACACACAACTGGCGTGTGAATATGGTCAAAAATCACTACCATCTCCCAAACTCGTTGCGGCCATCAACACCATCGTTAATGATTATGGCTATCTGCGTTGCTATAACACCGACGGCACgg                                                                                            <  2:433409/141‑1 (MQ=255)                                                   aaCTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACACACGCCACGGGcc                                                                                         >  2:687222/1‑143 (MQ=255)                                                      tGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCACACTGGTGGGGGCCATCACCACCATCGTTAATGATGATGGCTATCAGCGTGGCTAAAACACCGACAGCACGGGACata                                                                                      >  1:576967/1‑143 (MQ=255)                                                      tGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGAGGGCTATCTGCGTGGCTATAGCACCGACGGCATGGGCCata                                                                                      >  1:576805/1‑143 (MQ=255)                                                           tATCGATCCTGAACAATCAACTGTCGGGTGATTATTTTTATTAATTAACACCATCTGCCAAACTGGTGGGTGCCATCAACACCATCGTTATTGATGATGGCTATCTGCGTGTCTATAACACCGACGGCACGGGCCATATTcgc                                                                                 <  2:61581/143‑1 (MQ=255)                                                              cgaggacGAACAAATAACGGGGGTGTGAATATGTTAATGAATTAAAACCAGGTGCCAAACTGGTGGGGGGCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGcc                                                                              <  2:227189/137‑1 (MQ=255)                                                                         aaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAagagag                                                                   >  1:135132/1‑143 (MQ=255)                                                                         aaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATCACagagcg                                                                   >  1:500718/1‑141 (MQ=255)                                                                           tcgacgtgCGTGTGAATATGTTGATGAATTAACACCATCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCg                                                                 <  2:242884/136‑1 (MQ=255)                                                                               gtGGAGTGTGAATATGTAGATGAATTAATACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCCGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGttt                                                             <  1:573835/142‑1 (MQ=255)                                                                                tGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAAGGATGATGGCTATCTGCGTGGCTATAACACCGACGTCACGGGGAATATTCGGCCCATAAACGAAACCGGtttt                                                            >  1:501851/1‑143 (MQ=255)                                                                                  gCGTGTGAATATGTTGATGAATTAGCACCAGCTGCCGATCTGGTGGGGGCCATCAACACCATCGTTAATGATAATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGTCATTAAAGAGAGCGGTTTTGa                                                          <  2:220032/143‑1 (MQ=255)                                                                                    gtgtGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATCTTCGCGCCCTTAAGGAGAGCGGTTTTGata                                                        >  2:231076/1‑143 (MQ=255)                                                                                         aaTATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACCGCACGGGCCATATTTGGGCGATTAAAGAGAGCGGTATTGATATCaaa                                                   >  1:22520/1‑143 (MQ=255)                                                                                                 gatgaATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATCTTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAAc                                           >  1:626458/1‑143 (MQ=255)                                                                                                   tgaATTATCACCAGCTGCCAAACTGTTGGGGCCCATCAACACCATCGTTAGTGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGa                                         <  1:425209/143‑1 (MQ=255)                                                                                                     tatcagcacGAGTTGTCAAACTCGTGTGTGCCATCAACACCATCGTTAATGATGATGCCTTTCTGCGTGGCCACAACTCCGACGGCACTGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATg                                       <  2:228928/137‑1 (MQ=255)                                                                                                        tAACACCAGCTGCCAAACTGGTGGGGGCCATAAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGCGTTGATATAAAAGGCAAAACGATGGTg                                    >  2:626481/1‑143 (MQ=255)                                                                                                          acacCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGATTTGATATCAAAGGCAAAACGATGGTCCt                                  >  2:70012/1‑143 (MQ=255)                                                                                                                    ccAAACTGGTGGGGGCCAAAAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGAGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGATAGGGGcc                        >  2:640340/1‑143 (MQ=255)                                                                                                                        gcTGGTGGGGGCCGTCAACACCATCGTTGATGATGGTGGCTATCTGTGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCggtg                    <  2:203068/142‑1 (MQ=255)                                                                                                                           gatggGTGCCATAAACATCATCGTTAATGATTATAGCAATCTGCGTGGCAATAACACCGACGGCACGGGCCATATGCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGgtggtg                 <  1:167059/141‑1 (MQ=255)                                                                                                                            gttggggCCATAGACGCCTTCGTTAAGGATGAGGGTTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGc                <  2:607499/140‑1 (MQ=255)                                                                                                                                  ttcATTAACACTGTAGTTTATGATGTTTGCTATCTGTGTGGTTATAACACCGACGGCACGGGCCAGATTCGCGGCATCAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAAc          <  1:216606/141‑1 (MQ=255)                                                                                                                                     aTCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGc       <  1:578284/143‑1 (MQ=255)                                                                                                                                      tCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCa      >  1:424968/1‑143 (MQ=255)                                                                                                                                       cAACACCATCGTTAATCATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATCTCATAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCaa     <  2:518517/143‑1 (MQ=255)                                                                                                                                          cacCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTGTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCAGGAGGTGACTCAACCGCAAGTg  >  2:199204/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                             GATAGCTTTCCTGGAGCAATTGAAGGATTAAAAGCCCTCAAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTG  >  NC_000913/1773927‑1774206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ 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.