New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	202637 =	26 (0.760)	12 (0.360)	12/270	1.7	26.2%	coding (78/789 nt)	lpxA	UDP‑N‑acetylglucosamine acetyltransferase
	?	NC_000913	= 4542986	42 (1.250)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

GCTCTGGTCGATGGTAAAGTAGTTTGCGAAGCAACGATGATGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTCGCCTCCATCGGCGAAGTTAACCAGGATC  >  NC_000913/202491‑202782                                                                                                                                                   |                                                                                                                                                  gCTCTGGTCGATGGTAAAGTAGTTTGCGAAGCAACGATGATGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTg                                                                                                                                                  <  7:994/148‑1 (MQ=255)         gacgcatcttcctcattcttctctccaaaaaccacctcatgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCttttt                                                                                                                                         >  8:7450‑M1/139‑149 (MQ=255)         cGATGGTAAAGTAGTTTGCGAAGCAACGATGATGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCttttt                                                                                                                                         <  4:93570/149‑1 (MQ=255)                 ttcctcattcttctctccaaaaaccacctcatgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGtt                                                                                                                                 <  5:29959‑M1/19‑1 (MQ=255)                  tcctcattcttctctccaaaaaccacctcatgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTg                                                                                                                                <  5:93902‑M1/20‑1 (MQ=255)                  agtagtTTGCGAAGCAACGATGATGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTg                                                                                                                                <  5:19326/149‑1 (MQ=255)                                 aaCGATGATGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCCACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGaaa                                                                                                                 <  8:106967/149‑1 (MQ=255)                                  caaaaaccacctcatgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAAt                                                                                                                <  5:61202‑M1/36‑1 (MQ=255)                                    gatgatGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCTCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGCCCCCATGTCGAAATTg                                                                                                              >  8:26566/1‑149 (MQ=255)                                         cacctcatgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAg                                                                                                         >  8:113646‑M1/107‑149 (MQ=255)                                          acctcatgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAgg                                                                                                        >  6:93902‑M1/106‑149 (MQ=255)                                           tgtgCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAgg                                                                                                        <  5:118359/148‑1 (MQ=255)                                           tgtgCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCATTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAgg                                                                                                        <  7:127102/148‑1 (MQ=255)                                               atgcaatataaacatctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTAcc                                                                                                   >  6:29959‑M1/101‑149 (MQ=255)                                                          gAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAAtc                                                                                        <  5:57680/149‑1 (MQ=255)                                                          gAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAAtc                                                                                        <  5:55389/149‑1 (MQ=255)                                                            ggCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAAtctc                                                                                      <  2:148362/149‑1 (MQ=255)                                                            atctataaataaagataacaatagaatattaagccaacaaataaactgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAAtctc                                                                                      >  6:61202‑M1/88‑149 (MQ=255)                                                                                 aTCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTTGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATa                                                                 <  6:123975/149‑1 (MQ=255)                                                                                  tCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATAc                                                                >  3:118625/1‑149 (MQ=255)                                                                                  tCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATAc                                                                >  8:16362/1‑149 (MQ=255)                                                                                        tttGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAAt                                                          >  3:118092/1‑149 (MQ=255)                                                                                               aTCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCcgc                                                   >  1:137419/1‑149 (MQ=255)                                                                                                       gCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAAt                                             >  2:67075/1‑147 (MQ=255)                                                                                                         cATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAAt                                             >  5:136648/1‑145 (MQ=255)                                                                                                          tgaaaaagtttgtccgcgatgctttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGAt                                        <  4:75141‑M1/108‑1 (MQ=255)                                                                                                          aTTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGAt                                        >  4:36766/1‑149 (MQ=255)                                                                                                            aaaaagtttgtccgcgatgctttcctctatgattcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTAACGTACTGAAATCTCACGTTGTCGTGAAAGGTCATACTAAAATTGGCCGCGATAATGAGAtt                                       <  5:89626‑M1/109‑1 (MQ=255)                                                                                                                gAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTAt                                    >  6:73529/1‑147 (MQ=255)                                                                                                                gAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCa                                  <  3:36766/149‑1 (MQ=255)                                                                                                                    atggCGCGTCAATTGGCGCGAACGCAGAGATTGGTCCTTTTTGTATCGTTGGACACCATGTCGAAATTGGTGAGGGTACCGTACTTAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTc                              <  7:26566/147‑1 (MQ=255)                                                                                                                    agGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTc                              <  6:136648/149‑1 (MQ=255)                                                                                                                                tttcctctatgagtcaaaaTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTCGCCTCCATCGGc                  <  2:63898‑M1/130‑1 (MQ=255)                                                                                                                                  ggCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGa                                                                                            >  2:143444/1‑73 (MQ=255)                                                                                                                                   gcgcGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTCGCCTCCATCGGCGaa               <  4:118625/149‑1 (MQ=255)                                                                                                                                          cGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTCGCCTCCATCGGCGAAGTTAAcc        <  7:16362/149‑1 (MQ=255)                                                                                                                                             tcacaTTGGTCCTTTTTGTATCGTTGGACCCCATTTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAATTTGGCCGCGATAATGAGAt                                        <  5:100952/113‑1 (MQ=255)                                                                                                                                                taaTGTTCCTTTTTGTATCGTTGGACCCCATGTCGAATTTGGTGAGGGTACCGTACTGAAATCTCTCTTTGTCGTTAATTGTCATACTAAAATTGGCAGCGATAATGAGATTTATCATTTCGCCTCCATCGTCGAAGTTAACCAGGATc  >  6:89626‑M1/4‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCTCTGGTCGATGGTAAAGTAGTTTGCGAAGCAACGATGATGTGTGCTCGTAGCCGGGAGGCCTGATACGTGATTGATAAATCCGCCTTTGTGCATCCAACCGCCATTGTGGAAGAGGGCGCGTCAATTGGCGCGAACGCACACATTGGTCCTTTTTGTATCGTTGGACCCCATGTCGAAATTGGTGAGGGTACCGTACTGAAATCTCACGTTGTCGTGAATGGTCATACTAAAATTGGCCGCGATAATGAGATTTATCAGTTCGCCTCCATCGGCGAAGTTAACCAGGATC  >  NC_000913/202491‑202782

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