New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 580918 =38 (0.410)7 (0.080) 3/166 NT 16.6% coding (85/546 nt) nohD DLP12 prophage; DNA packaging protein
?NC_000913 = 1636304 35 (0.400)pseudogene (64/528 nt) nohQ pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin

TCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGC  >  NC_000913/580819‑581017
                                                                                                   |                                                                                                   
tCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGCGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCg                                                                                                     <  1:183624/100‑1 (MQ=38)
atccctttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccag                                                                                                     <  1:1878410‑M1/1‑1 (MQ=255)
 tccctttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccaGt                                                                                                    >  1:1427063‑M1/99‑100 (MQ=255)
  cGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGtt                                                                                                   >  2:41606/1‑100 (MQ=38)
   cctttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccaGTTc                                                                                                  >  1:603082‑M1/97‑100 (MQ=255)
   cctttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccaGTTc                                                                                                  <  2:2549814‑M1/4‑1 (MQ=255)
    ctttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccaGTTCt                                                                                                 >  1:1112285‑M1/96‑100 (MQ=255)
    ctttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccaGTTCt                                                                                                 <  1:783584‑M1/5‑1 (MQ=255)
     tttcggcataccatcttataacggcggcagagtcataaagcacctcattacccttgccaccgcctcgcagaacgggcattccctgttcctgccaGTTCTg                                                                                                >  2:1498846‑M1/95‑100 (MQ=255)
      gAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGc                                                                                               <  1:538163/100‑1 (MQ=38)
      gAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGc                                                                                               <  1:2355595/100‑1 (MQ=38)
          gtgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAgg                                                                                           <  2:722256/100‑1 (MQ=32)
           tgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGa                                                                                             >  2:1049820/1‑97 (MQ=31)
           tgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGc                                                                                          >  1:2378029/1‑100 (MQ=31)
            gCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCg                                                                                         >  1:1780606/1‑100 (MQ=11)
                                                     tGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGTGTAATGAGGTGCTTTATGACTCTGCCGCCGTc                                                <  1:2471878/100‑1 (MQ=11)
                                                         aGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATaa                                            <  1:1842613/100‑1 (MQ=32)
                                                           tATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATaaaa                                          >  2:2280243/1‑100 (MQ=35)
                                                             tCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATg                                        >  1:1631829/1‑100 (MQ=38)
                                                              ccGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATgg                                       >  2:2348320/1‑100 (MQ=38)
                                                              ccGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATgg                                       <  2:2188653/100‑1 (MQ=35)
                                                              ccGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATgg                                       <  2:2188357/100‑1 (MQ=38)
                                                               cGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGt                                      >  2:526700/1‑100 (MQ=38)
                                                                gTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTACTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTa                                     >  2:1391207/1‑100 (MQ=25)
                                                                 tACGATTCAGAACTGGCAGGAACAGGTAATGACCGTGTTGCGAGGCGGTCGCTAGGGTAATGAGGTGAATTAAGACATTGCCAACGTCATAAAATGGTa                                     <  2:2441583/99‑1 (MQ=2)
                                                                  gcctcgcagaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATg                                   <  1:1864711‑M1/67‑1 (MQ=255)
                                                                  gcctcgcagaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATg                                   <  1:1578834‑M1/67‑1 (MQ=255)
                                                                   ccATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGc                                  <  1:1420837/100‑1 (MQ=38)
                                                                    ctcgcagaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTTATAAGATGGTATGcc                                 <  1:578518‑M1/69‑1 (MQ=255)
                                                                     aTTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCg                                >  2:972401/1‑100 (MQ=38)
                                                                     aTTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCg                                <  1:45133/100‑1 (MQ=35)
                                                                       tCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGaa                              >  2:1546212/1‑100 (MQ=35)
                                                                       gcagaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTTATAAGATGGTATGCCGaa                              <  1:87952‑M1/72‑1 (MQ=255)
                                                                       gcagaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGaa                              <  1:2066369‑M1/72‑1 (MQ=255)
                                                                        cAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGCAATGCGGTGCTTTATGACTCTGCCGCCGTCATACAATGGTAt                                    >  2:2235760/1‑93 (MQ=255)
                                                                          gaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAgg                           <  1:1139490‑M1/75‑1 (MQ=255)
                                                                          gaacgggcattccctgttcctgccaGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAgg                           <  1:1143045‑M1/75‑1 (MQ=255)
                                                                          gAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAgg                           >  1:130413/1‑100 (MQ=38)
                                                                               ggCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGc                      >  1:1183555/1‑100 (MQ=38)
                                                                               ggCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGc                      >  1:1247634/1‑100 (MQ=38)
                                                                               ggCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGc                      <  1:2280256/100‑1 (MQ=38)
                                                                                   ggAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGaa                  >  2:95503/1‑100 (MQ=35)
                                                                                   ggAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGaa                  >  2:1653390/1‑100 (MQ=38)
                                                                                    gAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGaaa                 >  1:953115/1‑100 (MQ=38)
                                                                                    gAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGaaa                 >  2:2011204/1‑100 (MQ=35)
                                                                                      aCAGGGAATGCCCGTTCTGAGAGGCGGTGGCAAGAGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAAtt               <  2:519179/100‑1 (MQ=255)
                                                                                        aGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTga             <  2:1308732/100‑1 (MQ=38)
                                                                                         gggAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTgag            <  2:616163/100‑1 (MQ=35)
                                                                                         gggAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTgag            >  2:524466/1‑100 (MQ=38)
                                                                                         gggAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTgag            <  2:616225/100‑1 (MQ=38)
                                                                                         gggAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTgag            <  1:1049825/100‑1 (MQ=35)
                                                                                           gAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTACGCCGAAAGGGATGCTGAAATTGAGaa          >  2:705372/1‑100 (MQ=255)
                                                                                            aaTGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAAc         <  1:970290/100‑1 (MQ=38)
                                                                                             aTGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACg        <  1:747816/100‑1 (MQ=38)
                                                                                              tGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGa       >  2:2073317/1‑100 (MQ=38)
                                                                                                 ccGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGaaaa    <  1:19728/100‑1 (MQ=38)
                                                                                                  cGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAg   >  1:945904/1‑100 (MQ=38)
                                                                                                  cGTTCTGCGAGGCGGTGGCAAGGGAAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAg   >  2:264813/1‑100 (MQ=25)
                                                                                                  cGTTCTGCGAGGCGGTGGCAAGGGAAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCAGAAAGGGATGCTGAAATTGAGAACGAAAAg   >  2:2054022/1‑100 (MQ=255)
                                                                                                   GTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGc  >  2:1865928‑M1/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGC  >  NC_000913/580819‑581017

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