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.

AGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAA  >  NC_000913/3910184‑3910460                                                                                                                                   |                                                                                                                                                   aGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTAc                                                                                                                                  >  7:134783/1‑149 (MQ=255)    gcCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAg                                                                                                                               <  8:134783/149‑1 (MQ=255)                             cccACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTAc                                                                                                                                  >  1:209107/1‑121 (MQ=255)                             cccACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTAc                                                                                                                                  <  2:209107/121‑1 (MQ=255)                                      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                                                       <  2:32503/135‑1 (MQ=255)                                                                                          gatCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTActc                                                       >  1:32503/1‑135 (MQ=255)                                                                                                     cAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGt                               >  1:66617/1‑148 (MQ=255)                                                                                                      aataatGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGt                                        >  1:125523/1‑138 (MQ=255)                                                                                                      aataatGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGt                                        <  2:125523/138‑1 (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)                                                                                                                                   |                                                                                                                                                   AGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAA  >  NC_000913/3910184‑3910460

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