New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3495085	109 (1.940)	3 (0.050)	3/488	3.9	4.3%	coding (1075/2544 nt)	nirB	nitrite reductase (NADH) large subunit
	?	NC_000913	= 3803485	26 (0.460)					coding (653/1080 nt)	waaB	UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑1, 6‑D‑galactosyltransferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTAATCATTTGCTGTTTAATGCCAGAGCTGATTGCT  >  NC_000913/3803256‑3803730                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      aatGTCCGCAGGCCCGGAGACAAAATCGGCACTGATACAAGGACTCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGTTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCAAATACATgg                                                                                                                                                                                                                                   >  2:313136/1‑250 (MQ=255)             cccGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTAc                                                                                                                                                                                                                      >  2:91011/1‑251 (MQ=255)                  agaCACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGAACATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTc                                                                                                                                                                                                                 <  1:313136/251‑1 (MQ=255)                      acaATCGGCACTGATACACGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGCAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCtt                                                                                                                                                                                                             >  2:54397/1‑251 (MQ=255)                             gCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATaa                                                                                                                                                                                                      <  1:127388/251‑1 (MQ=255)                                                                    aGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCaa                                                                                                                                                                <  2:377854/250‑1 (MQ=255)                                                                              ggAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCACTATCTTGCTGGACCAATTACCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGc                                                                                                                                                     >  1:68909/1‑251 (MQ=255)                                                                              ggAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGCTATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCCTTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGc                                                                                                                                                     >  2:189029/1‑251 (MQ=255)                                                                                      tCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCACTATCTTGCTGCCCCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATACGCCTGGCATTTTTCGAAATCGGAGCCATAACCCAATACATGGAGTTTCAAGTTACCTTTAGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACat                                                                                                                                              >  1:426947/1‑250 (MQ=255)                                                                                          aaGATGACGTGCGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACataaat                                                                                                                                          >  1:351973/1‑250 (MQ=255)                                                                                          aaGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCACTATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCACCataaat                                                                                                                                          >  1:331853/1‑250 (MQ=255)                                                                                          aaGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACataaat                                                                                                                                          >  2:13667/1‑250 (MQ=255)                                                                                                      gCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGtt                                                                                                                             <  1:13667/251‑1 (MQ=255)                                                                                                                actcactTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCATCATAAATAAACGTTGCCGTTTCACCTTCCt                                                                                                                   <  1:56172/251‑1 (MQ=255)                                                                                                                 ctcactTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGc                                                                                                                                 >  2:30881/1‑236 (MQ=255)                                                                                                                 ctcactTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGc                                                                                                                                 <  1:30881/236‑1 (MQ=255)                                                                                                                                                       agagttgtgCCAACCAGACCAGTAAATTCGATCATCAATATTTAATGCTCTACCATAAGCCTGGCAGTTGTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATGTCATACGCCCAACATAAATAAACGTTGCCGTGTCACCTTCCTCGGGCGCGGGAATGACTGAGTTTTTAGTTGCGACGGGAt                                                                            <  2:85838/244‑1 (MQ=255)                                                                                                                                                                    ccATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTAACTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACAtt                                                               >  1:320734/1‑251 (MQ=255)                                                                                                                                                                        aCCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTAtc                                                           <  2:286976/251‑1 (MQ=255)                                                                                                                                                                        aCCAGACAATTCGATCATCAATATTTAAGTCGCTACCATAAGCCTTGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATGTCATACGCCCAACAGAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACTTGATTAAAAATAAAATTTAtc                                                           <  1:400055/251‑1 (MQ=255)                                                                                                                                                                                     atcatcAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTATCGTCGATTCTGCt                                               >  1:263555/1‑250 (MQ=255)                                                                                                                                                                                                       tctACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTACTCAttt                            >  2:98857/1‑251 (MQ=255)                                                                                                                                                                                                                aGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAACACATGGAGTTTCCAGTTAACTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTCATCATTTGCTGTTTa                    >  2:40266/1‑250 (MQ=255)                                                                                                                                                                                                                           tttCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTAATCATTTGCTGTTTAATGCCAGAGCt         <  1:18602/250‑1 (MQ=255)                                                                                                                                                                                                                             tCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGCCGGGCTTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTAATCATTTGCTGTTTAATGCCAGAGCTGAt      >  1:314774/1‑251 (MQ=255)                                                                                                                                                                                                                                 aaTCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGATCGTGATAAGCCATCAAGCCGGTCTTTAACTCTTTTCTGGCCTTAAAATTTCATACGCCCAACATAAATAAAAGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGCAGGGATTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTAATCATTTGCTGTTTAATGCCAGAGCTGATTGCt  >  1:423053/1‑251 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      AATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCAAGCAGGTCTTTAACTCTTTTCTGGCCTTCAAATTTCATACGCCCAACATAAATAAACGTTGCCGTTTCACCTTCCTCGGGCGCGGGAATGACTGAGTCTTTAGTTTCGACGGGATTAAAAATAACATTTATCGTCGATTCTGCTACACCTCGATTAATCATTTGCTGTTTAATGCCAGAGCTGATTGCT  >  NC_000913/3803256‑3803730

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