New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2258692 =	28 (0.970)	6 (0.210)	6/244	2.6	18.7%	coding (605/942 nt)	psuK	pseudouridine kinase
	?	NC_000913	2258736 =	25 (0.890)					coding (561/942 nt)	psuK	pseudouridine kinase
Rejected: Frequency below/above cutoff threshold.

ATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGG  >  NC_000913/2258547‑2258819                                                                                                                                                  |                                                                                                                                aTTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGcc                                                                                                                              >  3:125467/1‑149 (MQ=255)         cATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAgc                                                                                                                      >  7:25401/1‑149 (MQ=255)          aTTGGTTTTGCTCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAgcg                                                                                                                     >  6:122475/1‑149 (MQ=255)           ttGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCAc                                                                                                                                 <  6:25501/136‑1 (MQ=255)           ttGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCAc                                                                                                                                 >  5:25501/1‑136 (MQ=255)                 ttGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATcc                                                                                                               <  8:25401/148‑1 (MQ=255)                                gcatggattctggataatgccgccaacgtacccgtatttgtcgatccggtttcggcatggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGc                                                                                               <  8:31143‑M1/35‑1 (MQ=255)                                cAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGc                                                                                               >  5:83365/1‑149 (MQ=255)                                   ccACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTc                                                                                            <  4:125467/149‑1 (MQ=255)                                           aCCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCgg                                                                                     <  5:122475/148‑1 (MQ=255)                                            gataatgccgccaacgtacccgtatttgtcgatccggtttcggcatggaaaggtgtcaaagtgcgcgaccggctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGtt                                                                                   >  7:31143‑M1/103‑149 (MQ=255)                                                  atatCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCtt                                                                             <  6:118646/149‑1 (MQ=255)                                                  atatCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCtt                                                                             <  6:83365/149‑1 (MQ=255)                                                              taaACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGAt                                                                 <  4:122234/149‑1 (MQ=255)                                                                           atccggtttcggcatggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCgg                                                                                     <  4:108051‑M1/45‑1 (MQ=255)                                                                           atccggtttcggcatggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCgg                                                                                     >  3:108051‑M1/72‑116 (MQ=255)                                                                                          ggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGttt                                                                                  <  1:113587‑M1/48‑1 (MQ=255)                                                                                          ggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGttt                                                                                  >  2:113587‑M1/57‑104 (MQ=255)                                                                                             gTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACAt                                  >  6:70350/1‑149 (MQ=255)                                                                                                ggTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTc                               <  3:122102/149‑1 (MQ=255)                                                                                                      tgaaaTGTTGATGGAACCAGGAAGAAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCGTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCACGCACTTTGACACATTTCCATGc                          <  4:66234/144‑1 (MQ=255)                                                                                                      ggCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCGTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGCTTTAGACGGTCACGCACTTTGACACATTTCCATGc                          >  3:66234/1‑148 (MQ=255)                                                                                                               gATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGttt                                                                                  <  7:79124/83‑1 (MQ=255)                                                                                                               gATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGttt                                                                                  >  8:79124/1‑83 (MQ=255)                                                                                                               gATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGa                >  1:150722/1‑149 (MQ=255)                                                                                                                          gcagcaACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACAtt                                 >  6:45571/1‑121 (MQ=255)                                                                                                                          gcagcaACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACAtt                                 <  5:45571/121‑1 (MQ=255)                                                                                                                             gcaACTTTTGCCACATCTTCACGCCCTGCCAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACggg  >  8:62012/1‑149 (MQ=255)                                                                                                                                aCTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGAt               >  6:105408/1‑133 (MQ=255)                                                                                                                                aCTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGAt               <  5:105408/133‑1 (MQ=255)                                                                                                                                       ccACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTc                                                                                                  <  7:61082/43‑1 (MQ=255)                                                                                                                                       ccACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTc                                                                                                  >  8:61082/1‑43 (MQ=255)                                                                                                                                        cacaTCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                  <  6:117596/90‑1 (MQ=255)                                                                                                                                        cacaTCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                  >  5:117596/1‑90 (MQ=255)                                                                                                                                                  CGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                  >  4:43547‑M1/1‑80 (MQ=255)                                                                                                                                                  CGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                  <  3:43547‑M1/80‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                ATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGG  >  NC_000913/2258547‑2258819

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