New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	559162 =	35 (0.940)	5 (0.140)	5/250	4.1	12.8%	coding (189/693 nt)	sfmC	putative periplasmic pilus chaperone
	?	NC_000913	559214 =	35 (0.980)					coding (241/693 nt)	sfmC	putative periplasmic pilus chaperone
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAA  >  NC_000913/559019‑559308                                                                                                                                                |                                                                                                                                                   aTCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAAcag                                                                                                                                               >  3:55527/1‑149 (MQ=255)              gcgcCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTaaaa                                                                                                                                  <  4:55527/149‑1 (MQ=255)                           cGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtcattca                                                                                                                     <  2:80979/149‑1 (MQ=255)                             gAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATtcatc                                                                                                                   <  3:48935/149‑1 (MQ=255)                                                   gcgTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtca                                                                                                                         >  4:4237/1‑121 (MQ=255)                                                   gcgTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtca                                                                                                                         <  3:4237/121‑1 (MQ=255)                                                      tattatTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGt                                                                                           >  2:74097/1‑148 (MQ=255)                                                       attattTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTg                                                                                          >  3:148700/1‑148 (MQ=255)                                                           tttatCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAAcc                                                                                      >  4:85389/1‑148 (MQ=255)                                                              atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCaa                                                                                   >  5:131231/1‑148 (MQ=255)                                                                  cGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTac                                                                                                              >  4:57605/1‑117 (MQ=255)                                                                  cGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTac                                                                                                              <  3:57605/117‑1 (MQ=255)                                                                             aaCAGACTTCGGTATGGATTAGAAATAGCCAAACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACAGCCACtgttt                                                                                                 >  7:98047/1‑119 (MQ=255)                                                                                aGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTa                                                                                             <  6:152731/120‑1 (MQ=255)                                                                                aGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTa                                                                                             >  5:152731/1‑120 (MQ=255)                                                                                 gACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGta                                                               <  1:74097/149‑1 (MQ=255)                                                                                       ggTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTAttt                                                         <  3:85389/149‑1 (MQ=255)                                                                                        gTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTa                                                        >  1:105221/1‑149 (MQ=255)                                                                                           tGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTacac                                                     <  6:131231/149‑1 (MQ=255)                                                                                                  gAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATtcatca                                                                                                                  >  8:115455/1‑81 (MQ=255)                                                                                                  gAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATtcatca                                                                                                                  <  7:115455/81‑1 (MQ=255)                                                                                                     aTAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGc                                                                  >  2:51828/1‑126 (MQ=255)                                                                                                     aTAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGc                                                                  <  1:51828/126‑1 (MQ=255)                                                                                                     aTAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCg                                           >  4:58414/1‑149 (MQ=255)                                                                                                          cATACCAATGAGCGCATTCTGGTCAATTCTTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTgg                                       >  6:105334/1‑148 (MQ=255)                                                                                                            tACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGaaaa                                                                                                                                                  <  3:97725/39‑1 (MQ=255)                                                                                                            tACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGaaaa                                                                                                                                                  >  4:97725/1‑39 (MQ=255)                                                                                                                aaTGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACAcc                                                    >  8:123933/1‑129 (MQ=255)                                                                                                                aaTGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACAcc                                                    <  7:123933/129‑1 (MQ=255)                                                                                                                               gTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGtt                 >  5:55176/1‑149 (MQ=255)                                                                                                                                          gATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAAt       <  2:169126/148‑1 (MQ=255)                                                                                                                                          gATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAAt       >  1:169126/1‑148 (MQ=255)                                                                                                                                          gATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATg      >  8:66353/1‑149 (MQ=255)                                                                                                                                               caaaCAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTgg                                       <  3:42558/110‑1 (MQ=255)                                                                                                                                               aaaaCAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTgg                                       >  4:42558/1‑111 (MQ=255)                                                                                                                                               aaaaCAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTaa  <  3:58414/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   ATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAA  >  NC_000913/559019‑559308

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.