New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4576863 =45 (1.620)4 (0.160) 4/246 2.9 8.8% intergenic (+8/+49) yjiV/mcrC pseudogene; conserved hypothetical protein/5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
?NC_000913 4576890 = 41 (1.600)intergenic (+35/+22) yjiV/mcrC pseudogene; conserved hypothetical protein/5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Rejected: Frequency below/above cutoff threshold.

CTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGC  >  NC_000913/4576752‑4577034
                                                                                                                                          |                                                                                                                                                
atattggcttgtgtaccgtcaatttaggtcaggaatggccgtgtatacatcaagaattactcgacattttcgatgaatatctcaaataagcggtaagttttatcgaaaataCCGCTCCCGATGTGGCGTCGGGAGCGGTATTTTCtata                                                                                                                                        >  4:82266‑M1/112‑149 (MQ=255)
      aGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGATGTGATCGCCTCCTGATGATGAGCCGATCCCGATGTGGTTTCGGGAGCGGTATTTTATATAAAACt                                                                                                                                   >  6:108309/1‑148 (MQ=255)
            tCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCt                                                                                                                            >  1:125766/1‑149 (MQ=255)
              aaaaTGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTa                                                                                                                          >  3:26253/1‑149 (MQ=255)
                 aTGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt                                                                                                                       <  2:125766/149‑1 (MQ=255)
                    gCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgag                                                                                                                    <  5:108309/149‑1 (MQ=255)
                    gCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTga                                                                                                                     <  4:4531/148‑1 (MQ=255)
                        aGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGat                                                                                                                  <  4:26253/147‑1 (MQ=255)
                        aGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATAtt                                                                                                               >  6:21732/1‑149 (MQ=255)
                            cTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCAt                                                                                                            <  1:130446/149‑1 (MQ=255)
                            cTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTGATTTGAGATATTCAt                                                                                                            >  1:90707/1‑149 (MQ=255)
                                 aggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaaa                                                                                                      <  5:21732/149‑1 (MQ=255)
                                  ggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaaa                                                                                                      <  5:14749/149‑1 (MQ=255)
                                   gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTc                                                                                                              >  5:85101/1‑140 (MQ=255)
                                   gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTc                                                                                                              <  6:85101/140‑1 (MQ=255)
                                     aggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGt                                                                                                   <  8:81670/149‑1 (MQ=255)
                                     aggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGt                                                                                                   >  7:81670/1‑149 (MQ=255)
                                       ggAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                                                                  >  3:17840/1‑148 (MQ=255)
                                           tatacatcaataattactcgacattttcgatgaatatctcaaataagcggtaagttttatagaaaataCCGCTCCCGATGTGGTGTCTGGAGCGGTATTTTc                                                                                                                                            >  6:102834‑M1/69‑102 (MQ=255)
                                             tCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTTGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAtt                                                                                                                        <  6:117928/120‑1 (MQ=255)
                                             tCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAtt                                                                                                                        >  5:117928/1‑120 (MQ=255)
                                                 cAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt                                                                                       <  2:90707/149‑1 (MQ=255)
                                                  aaCAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt                                                                                       >  1:79557/1‑148 (MQ=255)
                                                   aCAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa                                                                                     <  2:79557/149‑1 (MQ=255)
                                                         aaGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                                                                  >  2:98876/1‑130 (MQ=255)
                                                         aaGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                                                                  <  1:98876/130‑1 (MQ=255)
                                                           ctcgacattttcgatgaatatctcaaataagcggtaagttttatagaaaataCCGCTCCCGATGTGGTGTCGGGAGCGGTAtt                                                                                                                                               <  1:114115‑M1/31‑1 (MQ=255)
                                                           ctcgacattttcgatgaatatctcaaataagcggtaagttttatagaaaataCCGCTCCCGATGTGGTGTCGGGAGCGGTAtt                                                                                                                                               >  2:114115‑M1/53‑83 (MQ=255)
                                                                ttGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAAt                                                                                                     <  8:2946/120‑1 (MQ=255)
                                                                ttGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtat                                                                                <  2:20021/141‑1 (MQ=255)
                                                                ttGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc                                                                         >  5:126655/1‑148 (MQ=255)
                                                                ttGGGCAACGGATGTGCTGGTGGAGGTGATCCCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtat                                                                                >  1:20021/1‑141 (MQ=255)
                                                                    gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAg                                                                                               >  8:24075/1‑122 (MQ=255)
                                                                    gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAg                                                                                               <  7:24075/122‑1 (MQ=255)
                                                                    gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTc                                                                    <  4:17840/149‑1 (MQ=255)
                                                                              tgCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAg                                                                                               >  4:64601/1‑112 (MQ=255)
                                                                                            aTCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTac                                                >  6:78250/1‑145 (MQ=255)
                                                                                            aTCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTac                                                <  5:78250/145‑1 (MQ=255)
                                                                                                    tatagaaaataccgctcccgacgccacatcgggagcggTATTTTCGATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAatat                                    >  4:82266‑M2/39‑149 (MQ=255)
                                                                                                      gatgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                                                                  >  2:47384/1‑85 (MQ=255)
                                                                                                      gatgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                                                                  <  1:47384/85‑1 (MQ=255)
                                                                                                        gaaaataccgctccagacaccacatcgggagcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTATTGATGtata                                                                               >  6:102834‑M2/35‑102 (MQ=255)
                                                                                                         gatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACgg                                                                          <  8:63636/106‑1 (MQ=255)
                                                                                                         gatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACgg                                                                          >  7:63636/1‑106 (MQ=255)
                                                                                                          atgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGc                              >  3:88160/1‑149 (MQ=255)
                                                                                                           tgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGc                              >  1:23393/1‑148 (MQ=255)
                                                                                                           aataccgctcccgacaccacatcgggagcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAg                                                                                               <  1:114115‑M2/52‑1 (MQ=255)
                                                                                                           aataccgctcccgacaccacatcgggagcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAg                                                                                               >  2:114115‑M2/32‑83 (MQ=255)
                                                                                                                     CCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAAtttt                    >  2:23393‑M1/1‑148 (MQ=255)
                                                                                                                       GATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGTAGCCATTAATTTtat                  <  2:101264‑M1/148‑1 (MQ=255)
                                                                                                                            ggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCa                                                                       >  8:105419/1‑90 (MQ=255)
                                                                                                                                  gggagcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCa      <  6:126655‑M2/141‑1 (MQ=255)
                                                                                                                                      gcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGc  <  4:88160‑M2/145‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGC  >  NC_000913/4576752‑4577034

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.