New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260479	3422061 =	18 (0.670)	7 (0.290)	5/188	1.7	28.4%	coding (609/933 nt)	h16_A3162	Periplasmic protein TonB, links inner and outer membranes
	?	AM260479	3422097 =	19 (0.790)					coding (573/933 nt)	h16_A3162	Periplasmic protein TonB, links inner and outer membranes

CGTTGTAGCCGTCGCGGTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCGGCCTCCAGC  >  AM260479/3421928‑3422175                                                                                                                                      |                                                                                                                   cGTTGTAGCCGTCGCGGTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGgcgc                                                                                                                    <  2:977216/134‑1 (MQ=255) cGTTGTAGCCGTCGCGGTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGgcgc                                                                                                                    <  2:1075181/134‑1 (MQ=255) cGTTGTAGCCGTCGCGGTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGgcgc                                                                                                                    <  1:852341/134‑1 (MQ=255)    caccccgcagcgcggccaggacgagcgcacctcgaccgacgagatggccaggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAgg                                                                                                                 >  2:746981‑M1/131‑134 (MQ=255)    caccccgcagcgcggccaggacgagcgcacctcgaccgacgagatggccaggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAgg                                                                                                                 >  2:10018‑M1/131‑134 (MQ=255)         cgcagcgcggccaggacgagcgcacctcgaccgacgagatggccaggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            >  1:505240‑M1/126‑134 (MQ=255)          cGTCGCGGTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCg                                                                                                           >  2:38180/1‑134 (MQ=255)                 gTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCgt                                                                                                          >  2:1018536/1‑128 (MQ=255)                 gTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCgt                                                                                                          <  1:1018536/128‑1 (MQ=255)                    gTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTgc                                                                                                  <  2:329365/133‑1 (MQ=255)                                 cctcgaccgacgagatggccaggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGc                                                                                    <  2:125880‑M1/33‑1 (MQ=255)                                     gaccgacgagatggccaggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCtgc                                                                                >  1:125880‑M1/98‑133 (MQ=255)                                                     aggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAgg                                                                                                                 <  2:887889‑M1/4‑1 (MQ=255)                                                     aggctggaggccgagatcggccgcaacctggaccactacgccaagcgacccaagcgcttccagctgaccgccaccagtgccCAgg                                                                                                                 >  1:887889‑M1/82‑85 (MQ=255)                                                                                  tGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATa                                                                                             >  1:1187910/1‑76 (MQ=255)                                                                                  tGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATa                                                                                             <  2:1187910/76‑1 (MQ=255)                                                                                    ccGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGgcgc                                                                                                                    <  1:305872/51‑1 (MQ=255)                                                                                    ccGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGgcgc                                                                                                                    >  2:305872/1‑51 (MQ=255)                                                                                             ccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            <  2:191564‑M1/9‑1 (MQ=255)                                                                                             ccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            <  2:968302‑M1/9‑1 (MQ=255)                                                                                             ccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            <  2:422610‑M1/9‑1 (MQ=255)                                                                                             ccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            >  1:191564‑M1/42‑50 (MQ=255)                                                                                             ccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            >  1:422610‑M1/42‑50 (MQ=255)                                                                                             ccaagcgacccaagcgcttccagctgaccgccaccagtgccCAGGCGGTc                                                                                                            >  1:968302‑M1/42‑50 (MQ=255)                                                                                               cgGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTg                       <  1:38180/133‑1 (MQ=255)                                                                                                          tggtgCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGc                                       <  2:675070/106‑1 (MQ=255)                                                                                                          tggtgCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGc                                       >  1:675070/1‑106 (MQ=255)                                                                                                          gcgcttccagctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCg           <  2:754317‑M1/106‑1 (MQ=255)                                                                                                          gcgcttccagctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCg           >  1:754317‑M1/29‑134 (MQ=255)                                                                                                           ggtgCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCgg          >  2:248818/1‑134 (MQ=255)                                                                                                             tgCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTgg                                        >  2:725186/1‑102 (MQ=255)                                                                                                             tgCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTgg                                        <  1:725186/102‑1 (MQ=255)                                                                                                               ccGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCGGCCTc      <  1:248818/134‑1 (MQ=255)                                                                                                                   gctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGa                <  2:2477‑M1/101‑1 (MQ=255)                                                                                                                   gctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGa                >  1:2477‑M1/20‑120 (MQ=255)                                                                                                                   gctgaccgccaccagtgccCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCGGCCTCCAGc  <  1:846870‑M1/115‑1 (MQ=255)                                                                                                                                      |                                                                                                                   CGTTGTAGCCGTCGCGGTTGTAGCCGAGCTTGCCGAGCCGGTTGACGTTGATCACCAGGATCAGCTGGCCGTACAGCGGCTTGCCGTTGCGTTGCGGGAAGTCAGTGGTGCCGCGCGCCTCGATCTTGTGGCGCAGGCGGTCGTAGTACTGCGCATAGTCGACCGCCTGGGCACTGGTGGCGGTCAGCTGGAAGCGCTTGGGTCGCTTGGCGTAGTGGTCCAGGTTGCGGCCGATCTCGGCCTCCAGC  >  AM260479/3421928‑3422175

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