New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4128795 =	46 (1.160)	5 (0.130)	5/258	4.0	10.0%	coding (124/1161 nt)	metB	cystathionine gamma‑synthase, PLP‑dependent
	?	NC_000913	4128835 =	46 (1.220)					coding (164/1161 nt)	metB	cystathionine gamma‑synthase, PLP‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTA  >  NC_000913/4128691‑4128981                                                                                                                                                 |                                                                                                                                                   ccaTCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGAtgtg                                                                                                                                                <  4:150217/149‑1 (MQ=255)  caTCGCATTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGGCGCGGCAACCCAACGCGCGATGTgg                                                                                                                                               <  5:163951/149‑1 (MQ=255)            cGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGCTCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTgc                                                                                                                                     >  7:152389/1‑149 (MQ=255)             gTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTgcg                                                                                                                                    >  6:1108/1‑149 (MQ=255)              tAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGc                                                                                                                                         >  1:57465/1‑143 (MQ=255)              tAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGc                                                                                                                                         <  2:57465/143‑1 (MQ=255)                         aTGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTg                                                                                                                        <  5:1108/149‑1 (MQ=255)                                   aCAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGgtggtg                                                                                                              >  5:120060/1‑149 (MQ=255)                                      gTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctg                                                                                                           <  6:120060/149‑1 (MQ=255)                                         tGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg                                                                                                        <  1:64572/149‑1 (MQ=255)                                            ttGCGTTGTCCCACCGATCCATCTTTCCAGCAACTATAACTATACCGGAATTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAg                                                                                                     >  2:52172/1‑149 (MQ=255)                                                gTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtact                                                                                                 >  7:93009/1‑149 (MQ=255)                                                 ttGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGc                                                                                                                                         >  2:18006/1‑108 (MQ=255)                                                 ttGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGc                                                                                                                                         <  1:18006/108‑1 (MQ=255)                                                   gTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtactta                                                                                               <  5:30173/148‑1 (MQ=255)                                                   gTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtactta                                                                                               <  8:152389/148‑1 (MQ=255)                                                    tCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCACCGCGCGATGTGGTTCa                                                                                                                                           >  2:101853/1‑103 (MQ=255)                                                    tCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCa                                                                                                                                           <  1:101853/103‑1 (MQ=255)                                                      ccACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGgtggtg                                                                                                              >  8:61809/1‑130 (MQ=255)                                                      ccACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGgtggtg                                                                                                              <  7:61809/130‑1 (MQ=255)                                                        aCCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATa                                                                                         >  7:155229/1‑149 (MQ=255)                                                             tCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCgg                                                                                     >  5:46148/1‑148 (MQ=255)                                                             tCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGc                                                                                    <  6:46148/149‑1 (MQ=255)                                                                aTCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgc                                                                                                             <  4:134244/121‑1 (MQ=255)                                                                aTCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgc                                                                                                             >  3:134244/1‑121 (MQ=255)                                                                aTCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtactta                                                                                               >  4:151313/1‑135 (MQ=255)                                                                aTCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtactta                                                                                               <  3:151313/135‑1 (MQ=255)                                                                aTCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATg                                                                                 >  5:177392/1‑149 (MQ=255)                                                                       cAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGa                                                                          <  8:155229/149‑1 (MQ=255)                                                                       cAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGa                                                                          >  1:80947/1‑149 (MQ=255)                                                                             cTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtacttact                                                                                             <  5:12389/124‑1 (MQ=255)                                                                             cTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtacttact                                                                                             >  6:12389/1‑124 (MQ=255)                                                                              tataACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCt                                                                   >  5:123412/1‑149 (MQ=255)                                                                               ataACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAggtg                                                                                                                 >  2:90927/1‑102 (MQ=255)                                                                               ataACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAggtg                                                                                                                 <  1:90927/102‑1 (MQ=255)                                                                                 aaCTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGGTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTc                                                                               <  8:69210/134‑1 (MQ=255)                                                                                 aaCTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTc                                                                               >  7:69210/1‑134 (MQ=255)                                                                                    tttACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAAc                                                             >  6:98855/1‑149 (MQ=255)                                                                                       aCCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGa                                                           <  6:123412/148‑1 (MQ=255)                                                                                                             cgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGAt                                    <  8:93009/149‑1 (MQ=255)                                                                                                              gcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATc                                   <  3:119022/149‑1 (MQ=255)                                                                                                                          gcgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGTAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTgcgc                       >  2:171100/1‑149 (MQ=255)                                                                                                                            gtcgcgGCAACCCAACGCGCGATGTGGTCCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCcgc                                                                            <  5:128203/94‑1 (MQ=255)                                                                                                                            gtcgcgGCAACCCAACGCGCGATGTGGTCCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCcgc                                                                            >  6:128203/1‑94 (MQ=255)                                                                                                                               gcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCAc                  >  7:11453/1‑149 (MQ=255)                                                                                                                                       ccAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTa          >  1:1237/1‑149 (MQ=255)                                                                                                                                           cgcgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGCCCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggc      >  7:123401/1‑149 (MQ=255)                                                                                                                                               cgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTa  <  2:1237/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTA  >  NC_000913/4128691‑4128981

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