New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3907801	35 (1.060)	4 (0.130)	4/264	3.9	11.1%	coding (566/774 nt)	pstB	phosphate ABC transporter ATPase
	?	NC_000913	= 3907849	31 (0.990)					coding (518/774 nt)	pstB	phosphate ABC transporter ATPase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGG  >  NC_000913/3907707‑3907996                                                                                                                                               |                                                                                                                                                    gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCggg                                                                                                                                                <  5:119778/148‑1 (MQ=255) gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGc                                                                                                                                               >  5:134306/1‑149 (MQ=255)     gAACAACGCGCAGGCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGTTAGAGATAGGGTCGAGCGCCGATCACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAAt                                                                                                                                           <  7:60155/149‑1 (MQ=255)     gAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAAt                                                                                                                                           <  8:92814/149‑1 (MQ=255)        caacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTAGTCGAGCAGCAGCACTTCCGGGCGAATGGc                                                                                                                                        >  1:88260/1‑149 (MQ=255)             cgcAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCg                                                                                                                                       >  8:129376/1‑145 (MQ=255)             cgcAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCg                                                                                                                                       <  7:129376/145‑1 (MQ=255)              gcAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATAcc                                                                                                                                  >  5:125836/1‑149 (MQ=255)                  ccTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCCGCACTTCCGGGCGAATGGCGATACCAcgc                                                                                                                              >  3:75207/1‑149 (MQ=255)                         cATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGGAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATa                                                                                                                        >  1:85885/1‑148 (MQ=255)                          aTGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTAGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATaca                                                                                                                      <  1:56627/149‑1 (MQ=255)                                    tGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATacac                                                                                                                     >  7:90191/1‑140 (MQ=255)                                    tGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATacac                                                                                                                     <  8:90191/140‑1 (MQ=255)                                             ccacGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACACTGTTGCTGacca                                                                                                   >  6:62234/1‑149 (MQ=255)                                                        tCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCagagaga                                                                                          >  8:130732/1‑147 (MQ=255)                                                        tCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGGCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCagagagaga                                                                                        <  2:33857/149‑1 (MQ=255)                                                        tCCTGCTTCAGATGTGTGATCATCTATTCAATACGCCCGGTAGAGATAGGGTCGAGAGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGTCGCTGTTGCTGACCACCagagagaga                                                                                        <  2:108173/149‑1 (MQ=255)                                                             cTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAAc                                                                                   >  5:68422/1‑149 (MQ=255)                                                                   ttCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGctct                                                                             <  6:68422/149‑1 (MQ=255)                                                                           tCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAAt                                                                     <  6:134306/149‑1 (MQ=255)                                                                             aGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAAtt                                                                    >  3:84827/1‑148 (MQ=255)                                                                                  ttCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCggg                                                                                                                                                <  5:139719/67‑1 (MQ=255)                                                                                  ttCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCggg                                                                                                                                                >  6:139719/1‑67 (MQ=255)                                                                                         cGCCCGGTAGAGACAGGGTCGAGCTCCGAACACGGGTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTc                                                       <  2:88260/149‑1 (MQ=255)                                                                                           cccGGTAGAGATAGGGTCGAGCGCCGAACACGGTGCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGt                                                     <  2:45491/149‑1 (MQ=255)                                                                                               gTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTCCGTTCca                                                 >  2:84397/1‑149 (MQ=255)                                                                                                   agaTAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCCCAAt                                             >  5:133684/1‑149 (MQ=255)                                                                                                    gaTAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATg                                            <  4:84827/149‑1 (MQ=255)                                                                                                      tAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGctct                                                                             <  8:33204/114‑1 (MQ=255)                                                                                                      tAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGctct                                                                             >  7:33204/1‑114 (MQ=255)                                                                                                             gAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGTCGATACCACGCGCAATACACAGACTCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTg                                    <  4:30471/148‑1 (MQ=255)                                                                                                             gAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACACTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTgg                                   <  5:62234/149‑1 (MQ=255)                                                                                                                  ccGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTATCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAAt                              <  7:130732/149‑1 (MQ=255)                                                                                                                           gTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTg                      >  1:23994/1‑148 (MQ=255)                                                                                                                                              ttCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCgg  >  4:117730/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                    GTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGG  >  NC_000913/3907707‑3907996

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