New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	3908495 =	43 (0.980)	6 (0.140)	6/272	5.0	13.7%	intergenic (‑129/+54)	pstB/pstA	phosphate ABC transporter ATPase/phosphate ABC transporter permease
	?	NC_000913_3_pae_tpiA	3908518 =	34 (0.800)					intergenic (‑152/+31)	pstB/pstA	phosphate ABC transporter ATPase/phosphate ABC transporter permease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCA  >  NC_000913_3_pae_tpiA/3908349‑3908640                                                                                                                                                   |                                                                                                                                                   aGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTcccc                                                                                                                                                  <  1:202203/149‑1 (MQ=255)     tCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACGTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑c                                                                                                                                             <  2:189342/149‑1 (MQ=255)       aaCCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTCCCATAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACGGTGTACAGAGTCAGAGAGGCGTTCGGTTGGTGCCCTCG‑gcc                                                                                                                                           >  1:575065/1‑149 (MQ=255)          cATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAACTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGCCGGCCCCCCCG‑CCCctc                                                                                                                                        >  2:137345/1‑149 (MQ=255)                 aTTGCAATCTTTTCCTCATGGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGgag                                                                                                                                 >  2:657895/1‑149 (MQ=255)                 aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCAGCATCCGGCAAAGTGCTGACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCGCTCG‑CCCCTCTGGGgag                                                                                                                                 >  2:658905/1‑149 (MQ=255)                   tGCAATCCTGTCCTCATCGGGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGGTCGAGTATACTTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGagag                                                                                                                               <  1:763819/149‑1 (MQ=255)                    gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAACGGGGAGGCAGGTCTTTCCCCTCG‑CCCCTCTGGGGAGAgg                                                                                                                              <  1:188854/149‑1 (MQ=38)                     cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAggg                                                                                                                             >  1:108220/1‑149 (MQ=255)                     cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACAGTAGGCGAGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGgggg                                                                                                                             >  1:613750/1‑149 (MQ=255)                               cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                   <  2:108220/149‑1 (MQ=255)                                  aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACTGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                >  2:617845/1‑149 (MQ=255)                                                  ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAc                                                                                                           <  2:166649/138‑1 (MQ=255)                                                  ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAc                                                                                                           >  1:166649/1‑138 (MQ=255)                                                         aTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                         >  2:293658/1‑149 (MQ=255)                                                             ggCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGAGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGGCGGTCCCCTCG‑CCACCCACAAGACATTCTTCGGGGCAGGGGAGAACCGTGTCAGCAATATCAACCgtg                                                                                     >  1:760585/1‑149 (MQ=25)                                                                     gTAGTCCGGATCAGGGGTTGCGGTCGCTTCCGGCAAGGTGTTAACGGTGTACAGGGTCAAAAGGGCGCTCGGTCGGCCCCCTCG‑CCCCTCTGGGGGGGGGGTTAGGGTGAGGGGAAAACTGTGTCAGCAATATCAACCGTGTTTAttct                                                                             <  1:657895/149‑1 (MQ=25)                                                                        ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAgg                                                                                                                              >  2:510141/1‑97 (MQ=255)                                                                        ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAgg                                                                                                                              <  1:510141/97‑1 (MQ=255)                                                                              tacgggCGTTCCCGCCTCTCCCGGAAAAGGGTTAACCGTGTACAGAGCCAAAAAGGCGCTGGGTCGGTCCCCCCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                    <  1:564217/145‑1 (MQ=37)                                                                                  ggCTTTCTCGCGGCACCCGGCAAAGTGTTACCCGTGCACAGAGTCAAAAAGGCGCGCGCTCGGCCCCCTCG‑CCCCTCGGGGGGGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                <  1:559797/149‑1 (MQ=37)                                                                                           gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAACACAACGCGCGCCac                                                       >  2:9873/1‑148 (MQ=255)                                                                                           gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                       >  1:255386/1‑149 (MQ=255)                                                                                            ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGACGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGGATCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACACGACCCAGa                                                      >  2:641696/1‑149 (MQ=38)                                                                                               cttCCGGCAAAGTGTGAACTGTGTCGAGAGTAAAAAAGTCGCTCGGTCGGCCCCCTCG‑CCCCTCGGGGGAGGGGGTTGGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATg                                                   <  1:570549/147‑1 (MQ=255)                                                                                                 gccGACAAAGTGTAAACCGTGTACAGAGCCAAAACGGCGCCCGGTCGGTCCCCTCC‑CCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                 <  1:534191/148‑1 (MQ=38)                                                                                                  ccGGCAAAGTGTTAACGGTGTACAGGGTCAAAAAGGCGCCGGGCCGGTCCCCGCC‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                <  2:343952/149‑1 (MQ=38)                                                                                                                   gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAg                               >  1:424167/1‑149 (MQ=255)                                                                                                                    tttACAGAGTAAAAAAAGGGCATGGGCGGCCCCCACG‑CCCCTAAGGGAAAGGGGTTAGGGTAAGGGGAAAACCGTATCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                              <  2:575065/147‑1 (MQ=255)                                                                                                                    tgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                              >  1:63541/1‑149 (MQ=255)                                                                                                                            gTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                      >  1:106585/1‑149 (MQ=255)                                                                                                                             tCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGGGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACgcgcgc                                                          >  2:490714/1‑112 (MQ=255)                                                                                                                             tCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGGGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACgcgcgc                                                          <  1:490714/112‑1 (MQ=255)                                                                                                                               aaaaaGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                   <  1:293658/149‑1 (MQ=255)                                                                                                                               aaaaaGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCACCCGTGTTTATTCTTCGCAAAAACACCGCGCGCCAGAATGTTAAGCCGCAGTACGCCCAGGGGAAGCATCaa                   >  2:40523/1‑149 (MQ=37)                                                                                                                                  aaGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAc                >  1:391069/1‑149 (MQ=255)                                                                                                                                        cTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGGTAGGGTGAGGAGAAAACAGTGACAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCACTACGCACAGGGTAATGATCAATACCCCgac          >  2:341640/1‑147 (MQ=39)                                                                                                                                         tcggtcggTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc         >  1:632993/1‑149 (MQ=255)                                                                                                                                           ggtcggtcCCCTCG‑CCCCTCTGGGGGGAGGGTTAGGGTGAGGCGAAAACCGTGACAGCAATAGCAACGGTGTTTATTGTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCa       >  2:631082/1‑149 (MQ=38)                                                                                                                                             tcggtcCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACCCACAGGGTAATGATCAATACCCCGGCCCAgg     >  2:131046/1‑149 (MQ=255)                                                                                                                                             tcggt‑cCCTCGCCCCCCCTGGGCGGAGGGTTAGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg     <  1:127530/149‑1 (MQ=255)                                                                                                                                               ctaacccTCT‑CCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc   >  2:754435/5‑149 (MQ=255)                                                                                                                                                gtcCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCa  >  1:369482/1‑149 (MQ=255)                                                                                                                                                gtcCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCa  >  1:369600/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCG‑CCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCA  >  NC_000913_3_pae_tpiA/3908349‑3908640

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