New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3463082 =	12 (0.670)	4 (0.230)	4/334	NT	22.4%	coding (159/984 nt)	gspK	Type II secretion system protein GspK
	?	NC_000913	3463296 =	16 (0.920)					coding (373/984 nt)	gspK	Type II secretion system protein GspK

AGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAGCACTGTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGA  >  NC_000913/3462825‑3463351                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              aGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCCGGGAAACCATCATCTGCCACGGCAATATGAATATGAAATGGCgcaacaaccttgg                                                                                                                                                                                                                                                                >  2:24466/1‑261 (MQ=255)    cTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTg                                                                                                                                                                                                                                                            <  1:24466/274‑1 (MQ=255)                                                                                                  cactatgccggggcggcgtctggctttggtaaaagctatcttctgcgccgtgaaagcgtggtgagtcatcgacatcaatatagtcagcgattgattgcacgatctcatcggtgtttccgcgatcgataccggcgttaatcagcagcgcactgaatacctgGGCGGAGAAACTTgcg                                                                                                                                                                                                                                                                <  2:64392‑M1/16‑1 (MQ=255)                                                                                                     tatgccggggcggcgtctggctttggtaaaagctatcttctgcgccgtgaaagcgtggtgagtcatcgacatcaatatagtcagcgattgattgcacgatctcatcggtgtttccgcgatcgataccggcgttaatcagcagcgcactgaatacctgGGCGGAGAAACTTGCGCtt                                                                                                                                                                                                                                                             >  1:64392‑M1/158‑176 (MQ=255)                                                                                                                                                gTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTgcgc                                                                                                                                                                                                                                                               >  2:240488/1‑131 (MQ=255)                                                                                                                                                  cgtgaaagcgtggtgggtcatcgacatcaatatagtcagcgattgattgcacgatctcatcggtgtttccgcgatcgataccggcgttaatcagcagcgcactgaatacctgGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTgg                                                                                                                                                                        <  2:257345‑M1/104‑1 (MQ=255)                                                                                                                                                   aCTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTg                                                                                                                                                                                                                                                            <  1:240488/131‑1 (MQ=255)                                                                                                                                                     gaaagcgtggtgggtcatcgacatcaatatagtcagcgattgattgcacgatctcatcggtgtttccgcgatcgataccggcgttaatcagcagcgcactgaatacctgGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAg                                                                                                                                                                     >  1:257345‑M1/110‑216 (MQ=255)                                                                                                                                                                        ggCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTAcgc                                                                                                                                        <  1:21639/226‑1 (MQ=255)                                                                                                                                                                                        cTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAg                                                                                                                               >  1:232401/1‑219 (MQ=255)                                                                                                                                                                                           aGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAGCAc                                                                                                                            <  2:232401/219‑1 (MQ=255)                                                                                                                                                                                                        agcTTCAGCGAACCCGCCAGGTAAACCCTCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAGCACTGTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCa                                                        <  1:41916/274‑1 (MQ=255)                                                                                                                                                                                                               gCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGtt                                                                                                                                                                                                              >  1:182357/1‑117 (MQ=255)                                                                                                                                                                                                                  aaCCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTgca                                                                                                                                                                                                           <  2:182357/117‑1 (MQ=255)                                                                                                                                                                                                                     ccGTCAGGTAAACCATCATCTGCAACTGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAGCACTGTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGt                                              <  1:82885/271‑1 (MQ=255)                                                                                                                                                                                                                     ccGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGcc                                                                                                                                                                                         >  1:206811/1‑132 (MQ=255)                                                                                                                                                                                                                        ccAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCcaca                                                                                                                                                                                      <  2:206811/132‑1 (MQ=255)                                                                                                                                                                                                                                                               aCTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAGCACTGTTTTAATCTCAATGCGTTAGCAACAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATga  >  1:214051/1‑273 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              AGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCGCTTGCCAGCCTCACCCAGGACGTTAAAGATAACGACAGACAGACCACGTTGCAGCAGTACTGGGCACAGCCACAACAGTTGCAGCTGGAAGATGGCAATACAGTGAAGTGGCAGTTACGCGATGCCCAGCACTGTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGA  >  NC_000913/3462825‑3463351

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