New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 4038798	31 (0.790)	4 (0.110)	4/268	4.7	12.4%	noncoding (1280/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
	?	NC_000913_3_pae_tpiA	= 4038805	27 (0.720)					noncoding (1287/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAG  >  NC_000913_3_pae_tpiA/4038651‑4038872                                                                                                                                                    |                                                                           aaCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGt                                                                           <  2:469085/149‑1 (MQ=255)      ggCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGaaaa                                                                      <  2:53280/149‑1 (MQ=40)         tAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGcc                                                                   <  2:95606/149‑1 (MQ=255)         tAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGcc                                                                   >  1:459761/1‑149 (MQ=255)                tGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGc                                                            >  1:441453/1‑149 (MQ=255)                    ccGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                        >  1:432749/1‑149 (MQ=255)                                 aGCGACACTAGGAGGTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTc                                           <  2:38929/149‑1 (MQ=255)                                     acacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGGCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGt                                       >  1:691360/1‑149 (MQ=17)                                      aacttagtgtTGTTGGGTAGGGGAGGGTTCGTTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATACGTAACGACAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAATTGTTCCTGTc                                      <  2:588553/143‑1 (MQ=11)                                        cTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGGGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATCAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGATCCttcgca                                    >  2:29888/1‑144 (MQ=255)                                               ttgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGACATACTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAAACCAAGGGTTCCGGTCCAACGTTaa                             >  2:466412/1‑149 (MQ=255)                                                 gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTAAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                       >  2:364537/1‑121 (MQ=18)                                                 gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTAAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                       <  1:364537/121‑1 (MQ=17)                                                   tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgg                         <  2:441453/149‑1 (MQ=32)                                                   tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAAGGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGCCCAAGGGTTCCTGTCCAACGTTAATCgg                         <  1:699326/149‑1 (MQ=17)                                                    gggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                        <  1:271484/149‑1 (MQ=32)                                                    gggCAGTGGAGCGTTGTGTATGCATGTGAAGGTGCGCTGCGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAACTAACGATAAAGCGGGTGAAAAGCGCGCTCGCCGGAAGACCAGGGTTTCCGGGCCAACGTTGATCggg                        <  2:574097/149‑1 (MQ=11)                                                     ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgggg                       >  2:11411/1‑149 (MQ=32)                                                     ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgggg                       >  2:47818/1‑149 (MQ=32)                                                       tAGGGGAGCGTTCCGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAGGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAGCAAGGGTTCCTGTCCAACGTTAATCGGGGCa                     >  2:280405/1‑149 (MQ=12)                                                        aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAg                    <  1:116432/149‑1 (MQ=32)                                                         ggggAGCGTTCTGTAAGCCTGTGAAGGTGTGCAGTGAGGCATACTGGAGGCATCAGAAGTGCGAACGCTGACATAGGTAACGATAAAGCGGGTGAAAATCCCGCTAGCCGGAAGACCGAGGGTTCCTGTCCAACGTTAATCGGGTCAgg                   >  2:662609/1‑149 (MQ=11)                                                           ggAGTGTTCTGTAAGGCTGTGAAGGGTTGCGGTGAGGCATGTGGGAGGTAGCAGAAGTGCTAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGt                 <  2:487298/149‑1 (MQ=11)                                                           ggAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGt                 <  2:579757/149‑1 (MQ=32)                                                            gAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg                >  1:228268/1‑149 (MQ=32)                                                            gAGCGTTCTGTAAGCCTGTGAAGGTGTGCCGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg                >  1:592555/1‑149 (MQ=17)                                                                  tCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGa          <  1:47699/149‑1 (MQ=32)                                                                    tGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGACGACGCAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc        >  1:306674/1‑149 (MQ=12)                                                                    tGTAAGCCTGTGAAGGTGTGCTGTGAGCCCTGCGTCGGTTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc        >  2:502790/1‑149 (MQ=11)                                                                     gTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAccc       <  2:459761/149‑1 (MQ=32)                                                                          ccTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAg  <  1:11411/149‑1 (MQ=32)                                                                                                                                                    |                                                                           AACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAG  >  NC_000913_3_pae_tpiA/4038651‑4038872

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