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.

CTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCA  >  NC_000913/3907671‑3907949                                                                                                                                   |                                                                                                                                                     cTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGTTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAgcg                                                                                                                                    >  3:30471/1‑149 (MQ=255)      tcaatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGaa                                                                                                                               >  1:33857/1‑149 (MQ=255)               ccAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcg                                                                                                                      >  1:45491/1‑149 (MQ=255)                 aGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgtc                                                                                                                    >  6:119778/1‑149 (MQ=255)                       aTAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgtc                                                                                                                    >  4:70783/1‑143 (MQ=255)                       aTAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgcc                                                                                                                    <  3:70783/143‑3 (MQ=255)                          aaCGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagcagca                                                                                                            <  1:99132/148‑1 (MQ=255)                                     gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCggg                                                                                                 <  5:119778/148‑1 (MQ=255)                                     gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGc                                                                                                >  5:134306/1‑149 (MQ=255)                                        ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTc                                                                                                                                         >  2:89894/1‑105 (MQ=255)                                        ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTc                                                                                                                                         <  1:89894/105‑1 (MQ=255)                                         gAACAACGCGCAGGCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGTTAGAGATAGGGTCGAGCGCCGATCACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAAt                                                                                            <  7:60155/149‑1 (MQ=255)                                         gAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAAt                                                                                            <  8:92814/149‑1 (MQ=255)                                            caacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTAGTCGAGCAGCAGCACTTCCGGGCGAATGGc                                                                                         >  1:88260/1‑149 (MQ=255)                                                 cgcAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCg                                                                                        <  7:129376/145‑1 (MQ=255)                                                 cgcAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCg                                                                                        >  8:129376/1‑145 (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                                                                      <  8:90191/140‑1 (MQ=255)                                                                        tGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATacac                                                                      >  7:90191/1‑140 (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                                                                                                 >  6:139719/1‑67 (MQ=255)                                                                                                                      ttCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCggg                                                                                                 <  5:139719/67‑1 (MQ=255)                                                                                                                             cGCCCGGTAGAGACAGGGTCGAGCTCCGAACACGGGTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTc        <  2:88260/149‑1 (MQ=255)                                                                                                                               cccGGTAGAGATAGGGTCGAGCGCCGAACACGGTGCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGt      <  2:45491/149‑1 (MQ=255)                                                                                                                                   gTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTCCGTTCca  >  2:84397/1‑149 (MQ=255)                                                                                                                                   |                                                                                                                                                     CTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCA  >  NC_000913/3907671‑3907949

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.