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.

ACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTC  >  NC_000913/2258589‑2258861
                                                                                                                                                   |                                                                                                                             
aCCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCgg                                                                                                                               <  5:122475/148‑1 (MQ=255)
 gataatgccgccaacgtacccgtatttgtcgatccggtttcggcatggaaaggtgtcaaagtgcgcgaccggctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGtt                                                                                                                             >  7:31143‑M1/103‑149 (MQ=255)
       atatCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCtt                                                                                                                       <  6:83365/149‑1 (MQ=255)
       atatCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCtt                                                                                                                       <  6:118646/149‑1 (MQ=255)
                   taaACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGAt                                                                                                           <  4:122234/149‑1 (MQ=255)
                                atccggtttcggcatggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCgg                                                                                                                               >  3:108051‑M1/72‑116 (MQ=255)
                                atccggtttcggcatggaaatgtgtcaaagtgcgcgaccgtctaaatcagatccacactctcaagccaaacCGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCgg                                                                                                                               <  4:108051‑M1/45‑1 (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                                                                                                                            >  8:79124/1‑83 (MQ=255)
                                                                    gATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGttt                                                                                                                            <  7:79124/83‑1 (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                                                         <  5:105408/133‑1 (MQ=255)
                                                                                     aCTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGAt                                                         >  6:105408/1‑133 (MQ=255)
                                                                                             cacaTCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                                                            >  5:117596/1‑90 (MQ=255)
                                                                                             cacaTCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                                                            <  6:117596/90‑1 (MQ=255)
                                                                                                   aaaccgccttgaagcggaaaccctgagtgggattgcgctgtcagggcgGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTcc                                                                        <  1:113587‑M2/56‑1 (MQ=255)
                                                                                                   aaaccgccttgaagcggaaaccctgagtgggattgcgctgtcagggcgGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTcc                                                                        >  2:113587‑M2/49‑104 (MQ=255)
                                                                                                    aaccgccttgaagcggaaaccctgagtgggattgcgctgtcagggcgGTTTGGCTTGAGAGTGTGGATCTGATTTAGCCGGTCGCGCACTTTGACACCTTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATc                          >  7:31143‑M2/48‑149 (MQ=255)
                                                                                                      ccgccttgaagcggaaaccctgagtgggattgcgctgtcagggcgGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGAt                                                         <  4:108051‑M2/71‑1 (MQ=255)
                                                                                                      ccgccttgaagcggaaaccctgagtgggattgcgctgtcagggcgGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGAt                                                         >  3:108051‑M2/46‑116 (MQ=255)
                                                                                                       CGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                                                            <  3:43547‑M1/80‑1 (MQ=255)
                                                                                                       CGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTc                                                                                            >  4:43547‑M1/1‑80 (MQ=255)
                                                                                                        GCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACaa                                                   >  8:54178‑M1/1‑120 (MQ=255)
                                                                                                        GCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACaa                                                   <  7:54178‑M1/120‑1 (MQ=255)
                                                                                                                gcggaaaccctgagtgggattgcgctgtcagggcgGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGc              <  8:31143‑M2/114‑1 (MQ=255)
                                                                                                                   gcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACgg                                             <  5:4859/115‑1 (MQ=255)
                                                                                                                   gcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACgg                                             >  6:4859/1‑115 (MQ=255)
                                                                                                                      aTCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGc        <  2:150722/149‑1 (MQ=255)
                                                                                                                           aCTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCtt   <  5:70350/149‑1 (MQ=255)
                                                                                                                                       cGCTTCAAGTCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACTCATTTCCATGCCGATACGGGATCGACAATTACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTc  <  3:66082/138‑1 (MQ=255)
                                                                                                                                        gCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCa                        <  1:119035/115‑1 (MQ=255)
                                                                                                                                        gCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCa                        >  2:119035/1‑115 (MQ=255)
                                                                                                                                                   gTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTAc                                         <  6:93301/87‑1 (MQ=255)
                                                                                                                                                   gTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTAc                                         >  5:93301/1‑87 (MQ=255)
                                                                                                                                                   |                                                                                                                             
ACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTC  >  NC_000913/2258589‑2258861

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.