New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	905931 =	34 (0.870)	6 (0.160)	6/268	3.9	15.9%	coding (823/1014 nt)	ybjS	putative NAD(P)H‑dependent oxidoreductase
	?	NC_000913	= 912978	31 (0.830)					coding (837/1653 nt)	hcp	hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATCGTCATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCT  >  NC_000913/912830‑913117                                                                                                                                                     |                                                                                                                                            gATCGTCATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAg                                                                                                                                             <  1:23196/149‑1 (MQ=255) gATCGTCATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTAACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAg                                                                                                                                             >  8:51555/1‑149 (MQ=255)      tCATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATg                                                                                                                                        >  1:175439/1‑149 (MQ=255)       cATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATgc                                                                                                                                       >  2:33318/1‑149 (MQ=255)       cATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATgc                                                                                                                                       >  6:39627/1‑149 (MQ=255)        aTAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATgcg                                                                                                                                      >  4:49553/1‑149 (MQ=255)         tAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATgcg                                                                                                                                      >  5:52128/1‑148 (MQ=255)             cgcCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCgcaggc                                                                                                                                 <  2:175439/149‑1 (MQ=255)                          ggcgctcatgCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTGCCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGt                                                                                                                    <  8:39689/142‑1 (MQ=255)                                    cAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACCCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGAc                                                                                                          <  5:39627/149‑1 (MQ=255)                                      gTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACAt                                                                                                        >  4:132253/1‑149 (MQ=255)                                          gAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAAc                                                                                                    <  6:52128/149‑1 (MQ=255)                                                     gATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCtt                                                                                         >  8:125723/1‑149 (MQ=255)                                                         gggCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGt                                                                                     <  8:104499/149‑1 (MQ=255)                                                               gggAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTc                                                                               <  3:49553/149‑1 (MQ=255)                                                                    aCGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCcagca                                                                          <  1:61285/149‑1 (MQ=255)                                                                                cACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGAAGCTTGAATTTACGCAGCTCCGGGAAGACATGCTCAGGCAGCAATACGCCGTGGGTGTAGCCATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGgtagta                                                                   >  6:33889/1‑144 (MQ=255)                                                                                   ttGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGCCATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGATATCtt                                                           >  3:90377/1‑149 (MQ=255)                                                                                   ttGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCtt                                                           >  6:115430/1‑149 (MQ=255)                                                                                   ttGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCtt                                                           >  6:114551/1‑149 (MQ=255)                                                                                    tGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCttt                                                          >  5:85316/1‑149 (MQ=255)                                                                                      cTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTa                                                                   >  7:12620/1‑138 (MQ=255)                                                                                      cTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTa                                                                   <  8:12620/138‑1 (MQ=255)                                                                                      cTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTga                                                        <  6:85316/149‑1 (MQ=255)                                                                                      cTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTga                                                        <  7:125723/149‑1 (MQ=255)                                                                                        gATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCttt                                                          <  2:61816/145‑1 (MQ=255)                                                                                        gATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCttt                                                          >  1:61816/1‑145 (MQ=255)                                                                                               ccAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGt                                                  <  3:132253/146‑1 (MQ=255)                                                                                                                    gACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTAc                          >  2:96053/1‑149 (MQ=255)                                                                                                                          aTGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACccgccg                    >  4:67623/1‑149 (MQ=255)                                                                                                                              ttGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGc                >  6:44949/1‑149 (MQ=255)                                                                                                                                   tttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGa           <  5:114551/149‑1 (MQ=255)                                                                                                                                   tttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGa           >  7:20244/1‑149 (MQ=255)                                                                                                                                   tttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGa           >  2:81650/1‑149 (MQ=255)                                                                                                                                        gCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGtt       >  2:13324/1‑148 (MQ=255)                                                                                                                                            cTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCt  <  8:20244/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                            GATCGTCATAAGCGCCTACGGTTGGGTCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCT  >  NC_000913/912830‑913117

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.