New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3907905	26 (0.700)	4 (0.120)	4/238	4.3	14.4%	coding (462/774 nt)	pstB	phosphate ABC transporter ATPase
	?	NC_000913	= 3907917	23 (0.660)					coding (450/774 nt)	pstB	phosphate ABC transporter ATPase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAAC  >  NC_000913/3907777‑3908061                                                                                                                                             |                                                                                                                                                 gTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTg                                                                                                                                          >  2:132416/1‑149 (MQ=255)                aTACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTg                                                                                                                                             >  4:93702/1‑131 (MQ=255)                aTACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTg                                                                                                                                             <  3:93702/131‑1 (MQ=255)                                  gggTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCgg                                                                                                          <  2:44633/148‑1 (MQ=255)                                    gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCaa                                                                                                                                       <  4:137516/117‑1 (MQ=255)                                    gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCaa                                                                                                                                       >  3:137516/1‑117 (MQ=255)                                    gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCtt                                                                                                                               >  8:96173/1‑125 (MQ=255)                                    gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCtt                                                                                                                               <  7:96173/125‑1 (MQ=255)                                    gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCtt                                                                                                       <  1:132416/149‑1 (MQ=255)                                        aGCGCCGATCACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGt                                                                                                   <  4:79956/149‑1 (MQ=255)                                                      ttCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCa                                                                                     >  3:131101/1‑149 (MQ=255)                                                         gtcgAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCAcgc                                                                                  >  3:4789/1‑149 (MQ=255)                                                              gcagcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCg                                                                             >  3:72066/1‑149 (MQ=255)                                                                agcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTgg                                                                                                    >  5:160667/1‑124 (MQ=255)                                                                agcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTgg                                                                                                    <  6:160667/124‑1 (MQ=255)                                                                 gcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTc                                                                           >  2:93706/1‑148 (MQ=255)                                                                      aCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTc                                                                     <  4:4789/149‑1 (MQ=255)                                                                       cTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCg                                                                    >  8:106707/1‑149 (MQ=255)                                                                       cTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCg                                                                    <  4:131101/149‑1 (MQ=255)                                                                                  aTGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGagag                                                         <  4:72066/149‑1 (MQ=255)                                                                                                              gCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGga                                                            <  4:21615/118‑1 (MQ=255)                                                                                                              gCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGga                                                            >  3:21615/1‑118 (MQ=255)                                                                                                                          caccaGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCg                    >  7:23117/1‑146 (MQ=255)                                                                                                                                   gagagTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGTCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATc        <  7:106707/149‑1 (MQ=255)                                                                                                                                      agTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTACACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCgg      >  4:75408/1‑148 (MQ=255)                                                                                                                                          aCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAAc  <  1:93706/148‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                 GTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAAC  >  NC_000913/3907777‑3908061

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