New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4576350 =	36 (0.570)	3 (0.050)	3/252	NT	8.8%	coding (747/795 nt)	ybhP‑II	phosphohydrolase
	?	NC_002947	4576394 =	28 (0.470)					coding (791/795 nt)	ybhP‑II	phosphohydrolase

GCAACGCGACTTGCACGAAGCCTTCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGTTGAGCGTTACGAC  >  NC_002947/4576203‑4576490                                                                                                                                                    |                                                                                                                                             gagaggtcgcggcttgcgtcaagaaagcgtgaaggaagtgttgagggtgcgtcgtaacgctcaacctgttgttttataggtaaaaacagcagggtatagcggggatcgctatttgccgtgccgtggctatgaagaatgattgctaaaGt                                                                                                                                             <  2:124258‑M1/2‑1 (MQ=255)  cAACGCGACTTGCACGAAGCCTTCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTgg                                                                                                                                              <  5:124271/147‑1 (MQ=255)   aaCGCGACTTGCACGAAGCCTTCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATGCCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCg                                                                                                                                           >  5:43961/1‑149 (MQ=255)    aggtcgcggcttgcgtcaagaaagcgtgaaggaagtgttgagggtgcgtcgtaacgctcaacctgttgttttataggtaaaaacagcagggtatagcggggatcgctatttgccgtgccgtggctatgaagaatgattgctaaaGTCGc                                                                                                                                          >  1:52231‑M1/145‑149 (MQ=255)                cGAAGCCTTCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGAt                                                                                                                               >  7:85008/1‑148 (MQ=255)                  aaGCCTTCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCAc                                                                                                                            <  8:104971/149‑1 (MQ=255)                       ttCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGcc                                                                                                                       >  1:191543/1‑149 (MQ=255)                       ttCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGcc                                                                                                                       >  8:119466/1‑149 (MQ=255)                           aGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTg                                                                                                                   >  8:224051/1‑149 (MQ=255)                            gcgcCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATGCCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTgg                                                                                                                  <  6:43961/149‑1 (MQ=255)                              gcCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGAATCTGGGGCACAAACCTTGGTCGCACCTGTCAGATAACCTGCCGCAGGcc                                                                                                                >  4:67506/1‑149 (MQ=255)                                cacctcGGCCACCTCGCGCGCACGTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGGACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGt                                                                                                              <  3:67506/149‑1 (MQ=255)                                   cacGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTg                                                                                                            >  7:3616/1‑148 (MQ=255)                                      ggCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGt                                                                                                        <  2:191543/149‑1 (MQ=255)                                          aCCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACCAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGGTGAAGTGCGa                                                                                                    >  5:174666/1‑149 (MQ=255)                                          aCCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGa                                                                                                    >  7:294259/1‑149 (MQ=255)                                             gggtgcgtcgtaacgctcaacctgttgttttataggtaaaaacagcagggtatagcggggatcgctatttgccgtgccgtggctatgaagaatgattgctaaaGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACtt                                                                                                 >  1:393344‑M1/104‑149 (MQ=255)                                                       cTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCTTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCAt                                                                                       <  4:382782/149‑1 (MQ=255)                                                          cccGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCAtt                                                                                      <  7:224051/147‑1 (MQ=255)                                                           ccGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAg                                                                                              <  6:18214/138‑1 (MQ=255)                                                           ccGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAg                                                                                              >  5:18214/1‑138 (MQ=255)                                                            ccgCCCCCCAGCCGCTGCTACGCCTGGACCGCAAATACCTTCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCAttcttc                                                                                  <  4:361774/147‑1 (MQ=255)                                                            cGGCCCGCCTGCCGCTGCTACGCCAGGACCTCATCTACCTGCGTAACGCCGAAAGCCAGGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTCTCAGATCACCTGCCGCGGGCCGCTGAAGTGCGACTTTAGCAATCAttcttc                                                                                  <  6:384985/149‑1 (MQ=255)                                                             ggCCCGCATGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGACGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTTGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCa                                                                                 <  3:273700/149‑1 (MQ=255)                                                               cccGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATa                                                                               <  7:17494/149‑1 (MQ=255)                                                                   tgttgttttataggtaaaaacagcagggtatagcggggatcgctatttgccgtgccgtggctatgaagaatgattgctaaaGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGt                                                                                                        <  4:76374‑M1/39‑1 (MQ=255)                                                                   tgttgttttataggtaaaaacagcagggtatagcggggatcgctatttgccgtgccgtggctatgaagaatgattgctaaaGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGt                                                                                                        >  3:76374‑M1/82‑120 (MQ=255)                                                                          ctgctACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCcacggcacg                                                                    <  6:298207/149‑1 (MQ=255)                                                                              tACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACGTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCAcggcacggca                                                                 >  1:124258/1‑148 (MQ=255)                                                                                           cATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCg                                                   >  5:25552/1‑149 (MQ=255)                                                                                           cATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCg                                                   >  3:237644/1‑149 (MQ=255)                                                                                            atatACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAg                                                            >  4:111711/1‑139 (MQ=255)                                                                                            atatACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAg                                                            <  3:111711/139‑1 (MQ=255)                                                                                                    gCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCct                                          >  2:336291/1‑149 (MQ=255)                                                                                                           gCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGttt                                   <  4:237644/149‑1 (MQ=255)                                                                                                           gCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGttt                                   <  8:85008/149‑1 (MQ=255)                                                                                                           gCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTAGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGttt                                   >  5:438/1‑149 (MQ=255)                                                                                                                      tGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATaaa                        <  6:25552/149‑1 (MQ=255)                                                                                                                               gATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGt               <  1:336291/149‑1 (MQ=255)                                                                                                                               gATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGt               <  3:135422/149‑1 (MQ=255)                                                                                                                                      gggCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGTTGAGc          <  4:125877/147‑1 (MQ=255)                                                                                                                                      gggCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGTTGAGCGt        >  8:347813/1‑149 (MQ=255)                                                                                                                                            aaaCCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGTTGAGCGTTacgac  <  8:3616/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             GCAACGCGACTTGCACGAAGCCTTCGAGCGCCACCACGGCCACCTCGCGCGCACCTACCCGGCCCGCCTGCCGCTGCTACGCCTGGACCGCATATACCTGCGTAACGCCGAAAGCCATGGGCCGCGGATTCTGGGGCACAAACCTTGGTCGCACCTGTCAGATCACCTGCCGCTGGCCGTTGAAGTGCGACTTTAGCAATCATTCTTCATAGCCACGGCACGGCAAATAGCGATCCCCGCTATACCCTGCTGTTTTTACCTATAAAACAACAGGTTGAGCGTTACGAC  >  NC_002947/4576203‑4576490

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