New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3910314 =	32 (0.770)	5 (0.130)	5/246	4.2	13.1%	coding (85/960 nt)	pstC	phosphate ABC transporter permease
	?	NC_000913	3910344 =	36 (0.910)					coding (55/960 nt)	pstC	phosphate ABC transporter permease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTA  >  NC_000913/3910221‑3910490                                                                                                                            |                                                                                                                                                   tGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTAccc                                                                                                                              <  8:119810/146‑1 (MQ=255)   gTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTgg                                                                                                                         >  5:95030/1‑149 (MQ=255)   gTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAACAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTgg                                                                                                                         <  1:33603/149‑1 (MQ=255)     ccATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGt                                                                                                                       >  1:140610/1‑149 (MQ=255)      cATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGt                                                                                                                       <  6:95030/148‑1 (MQ=255)      cATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGt                                                                                                                       >  5:109316/1‑148 (MQ=255)              aGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGGGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAg                                                                                                              >  2:31882/1‑149 (MQ=255)                   gACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTAt                                                                                                                                              <  6:123491/112‑1 (MQ=255)                   gACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTAt                                                                                                                                              >  5:123491/1‑112 (MQ=255)                    aCCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTgg                                                                                                        <  6:55525/149‑1 (MQ=255)                        aaTTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTg                                                                                                     >  8:74713/1‑148 (MQ=255)                                aaTGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCAt                                                                                              >  3:49146/1‑147 (MQ=255)                                     tCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTAc                                                                                       >  5:26249/1‑149 (MQ=255)                                            ggaggaGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAg                                                                                <  4:49146/149‑1 (MQ=255)                                                     gatCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTActc                                                                                     >  1:32503/1‑135 (MQ=255)                                                     gatCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTActc                                                                                     <  2:32503/135‑1 (MQ=255)                                                                cAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGt                                                             >  1:66617/1‑148 (MQ=255)                                                                 aataatGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGt                                                                      <  2:125523/138‑1 (MQ=255)                                                                 aataatGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGt                                                                      >  1:125523/1‑138 (MQ=255)                                                                 aataatGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCa                                                           >  1:152347/1‑149 (MQ=255)                                                                     atGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCAtttt                                                       <  2:140610/149‑1 (MQ=255)                                                                        ccACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTa                                                    <  6:109316/149‑1 (MQ=255)                                                                          aCCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCAtttggatcg                                                  >  1:156784/1‑143 (MQ=255)                                                                                  tCAATAGCACAATCAGCGCCGCCAGTTTTACAAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATg                                          >  6:65340/1‑149 (MQ=255)                                                                                   cAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATgc                                         <  7:74713/149‑1 (MQ=255)                                                                                     aTAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATgcgg                                       >  2:135327/1‑149 (MQ=255)                                                                                     aTAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATgcgg                                       >  2:25450/1‑149 (MQ=255)                                                                                            aaTCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTaaaa                                >  8:70831/1‑149 (MQ=255)                                                                                                     cgccAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCtt                                                                                   >  8:53012/1‑89 (MQ=255)                                                                                                     cgccAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCtt                                                                                   <  7:53012/89‑1 (MQ=255)                                                                                                       ccAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCg                     >  4:154654/1‑149 (MQ=255)                                                                                                         aGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGc                   <  2:152347/149‑1 (MQ=255)                                                                                                          gTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGcc                  <  1:31882/149‑1 (MQ=255)                                                                                                           ttttACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCt                 <  5:65340/149‑1 (MQ=255)                                                                                                              tACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTgg               <  6:26249/148‑1 (MQ=255)                                                                                                                 cAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGtt           <  2:66617/149‑1 (MQ=255)                                                                                                                      cGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttatt      <  6:117070/149‑1 (MQ=255)                                                                                                                       gCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttatta     <  3:154654/149‑1 (MQ=255)                                                                                                                          gAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTa  >  1:153938/1‑149 (MQ=255)                                                                                                                           aaaaTTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTa  <  1:135327/148‑1 (MQ=255)                                                                                                                            |                                                                                                                                                   TGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTA  >  NC_000913/3910221‑3910490

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