New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 58089659 (0.450)9 (0.070) 7/166 NT 16.5% coding (63/546 nt) nohD DLP12 prophage; DNA packaging protein
?NC_000913 1636284 = 36 (0.290)pseudogene (84/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

CTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCC  >  NC_000913/1636185‑1636383
                                                                                                   |                                                                                                   
cTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACg                                                                                                     <  2:228092/100‑1 (MQ=38)
 ttttCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACgn                                                                                                    <  1:3637522/100‑2 (MQ=38)
   ttagTTCTCAATTTCAGCATCCCTTTCGTCATATTATCTTATAACCACGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAATgg                                                                                                    >  1:1030881/4‑98 (MQ=11)
   ttCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGc                                                                                                  <  1:1517300/100‑1 (MQ=38)
       ttCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTc                                                                                              >  2:486564/1‑100 (MQ=38)
         ctcAATTTCATCGTCCCTTTCGGCATACCATCTTATAACCGCAACAAAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTccc                                                                                            >  1:1659264/1‑100 (MQ=11)
         ctcAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTccc                                                                                            >  1:372130/1‑100 (MQ=38)
           cAATTTCATCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGCGCCTTCCCTg                                                                                          >  1:1683200/1‑100 (MQ=255)
           cAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGTGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTg                                                                                          >  2:117627/1‑100 (MQ=255)
           cAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATACAGCACCTCATTACACTTGCAACCGCCTCGCAGAACGGGCATCCCCt                                                                                           >  1:1289129/1‑99 (MQ=11)
           cAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTg                                                                                          >  2:1694074/1‑100 (MQ=38)
               tttttgtccgtggaatgtgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacTGGCATTCCCTGTTcc                                                                                      <  1:798965‑M2/16‑1 (MQ=255)
                 cagcagCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTg                                                                                    <  2:3565194/100‑1 (MQ=38)
                 cAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTg                                                                                    >  2:1236283/1‑100 (MQ=38)
                 cAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTg                                                                                    >  1:2554260/1‑100 (MQ=38)
                 cAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTg                                                                                    <  1:1606396/100‑1 (MQ=38)
                  aGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGc                                                                                   >  1:3350300/1‑100 (MQ=38)
                  aGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAACACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGc                                                                                   <  1:1508419/100‑1 (MQ=25)
                   gCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGcc                                                                                  >  2:2043677/1‑100 (MQ=38)
                     aTCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAg                                                                                <  2:1138598/100‑1 (MQ=38)
                      tCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCCTTCCCTGTTCCTGCCAGt                                                                               >  1:491312/1‑100 (MQ=38)
                        gtggaatgtgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTc                                                                             <  1:3281365‑M2/25‑1 (MQ=255)
                        ccTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTc                                                                             >  2:2944545/1‑100 (MQ=38)
                         cTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCt                                                                            <  2:1390663/100‑1 (MQ=35)
                         cTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCt                                                                            <  1:3380775/100‑1 (MQ=38)
                          tttCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTg                                                                           <  2:3371543/100‑1 (MQ=35)
                          ggaatgtgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCt                                                                            <  2:242981‑M2/26‑1 (MQ=255)
                            aatgtgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCTGaa                                                                         <  2:2870920‑M2/29‑1 (MQ=255)
                              tgtgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCTGAATg                                                                       <  2:1967649‑M2/31‑1 (MQ=255)
                               gtgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCTGAATgg                                                                      <  1:846027‑M2/32‑1 (MQ=255)
                               gaacaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCTGAATgg                                                                      <  2:2316892‑M2/32‑1 (MQ=255)
                               gaacaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCTGAATgg                                                                      <  2:2158804‑M2/32‑1 (MQ=255)
                               gCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATgg                                                                      >  1:1888738/1‑100 (MQ=38)
                                tgcaatggaagtcaacaaaaagcagctggctgacattttcggtgcgagtatccgtaccattcagaacGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGt                                                                     <  1:1416037‑M2/33‑1 (MQ=255)
                                  tACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTAc                                                                   <  1:513571/100‑1 (MQ=38)
                                  tACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTAc                                                                   <  2:2396034/100‑1 (MQ=38)
                                  tACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTAc                                                                   <  2:24004/100‑1 (MQ=35)
                                   aCCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACg                                                                  >  2:1273265/1‑100 (MQ=38)
                                     cATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGa                                                                <  1:251269/100‑1 (MQ=38)
                                     cATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGa                                                                <  2:3210599/100‑1 (MQ=38)
                                          nagaaCGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCGAACGCCTCGCCGAACGCGCATTCACTGTTACTCACAATTCTGAATGGTCCGGCTAc                                                             >  1:1606994/4‑98 (MQ=11)
                                           ataACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCg                                                          >  1:726599/1‑100 (MQ=32)
                                           ataACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCg                                                          >  1:723623/1‑100 (MQ=32)
                                           ataACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCg                                                          >  1:417808/1‑100 (MQ=32)
                                            taACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGc                                                         <  2:32008/100‑1 (MQ=31)
                                                                           aCCCTTGCCACCGCATCGCAGAACGGGCATTCCCTGTTACTGCCAGTTCTGCATGATACGAATGCTCGCGCCGACAATGTCATGCAAACGCGTTttgttg                          >  1:1333342/1‑100 (MQ=0)
                                                                                     cccccTCGCAGAGCGGCCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCAt                  <  1:707894/95‑1 (MQ=0)
                                                                                      cGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGt               <  1:1009892/100‑1 (MQ=11)
                                                                                      cGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGt               <  1:3040430/100‑1 (MQ=11)
                                                                                       gCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGtt              >  2:3224713/1‑100 (MQ=31)
                                                                                         cTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGtt              <  1:1883135/98‑1 (MQ=31)
                                                                                          tCGCAGAACGGGGATTCCCTGTTCCTGCCAGTTCTGACTGGTACGGATGCTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTc             <  1:3513179/98‑1 (MQ=255)
                                                                                          tCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCAt           <  2:193391/100‑1 (MQ=35)
                                                                                           cGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCAtt          >  1:1343880/1‑100 (MQ=35)
                                                                                           cGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCAtt          >  2:2639935/1‑100 (MQ=35)
                                                                                            gCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTc         <  2:1914132/100‑1 (MQ=38)
                                                                                            gCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTc         >  2:313922/1‑100 (MQ=38)
                                                                                            gCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTc         <  1:2932046/100‑1 (MQ=38)
                                                                                            gCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTc         >  1:2308082/1‑100 (MQ=38)
                                                                                             cagaacGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTcc        <  2:945887‑M2/94‑1 (MQ=255)
                                                                                                aacGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACg     >  2:1210583‑M2/4‑100 (MQ=255)
                                                                                                 acGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACgg    >  2:436449‑M2/3‑100 (MQ=255)
                                                                                                  cGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGc   >  2:1499350‑M2/2‑100 (MQ=255)
                                                                                                   gggCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTTTTCATTCCACGGcc  >  1:144356/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCC  >  NC_000913/1636185‑1636383

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.