New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =11 (0.160)51 (0.880) 41/144 NT 87.9% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 5 (0.090)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAT  >  NC_000913/1209524‑1209717
                                                                                               |                                                                                                  
ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                 >  1:1513662‑M2/96‑99 (MQ=255)
ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                 <  1:1399846‑M2/4‑1 (MQ=255)
  caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                               <  2:658311‑M2/6‑1 (MQ=255)
  aaataccaccaatcataaccggtatagctatgaaaccggtgcatcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                               <  1:1941321‑M2/6‑1 (MQ=255)
     taccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                                                            >  1:1150340‑M2/91‑99 (MQ=255)
      accaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                                                           <  1:231049‑M2/10‑1 (MQ=255)
                   accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                              >  2:396002‑M2/77‑99 (MQ=255)
                   accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                              <  1:1527900‑M2/23‑1 (MQ=255)
                    ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                             >  1:1332655‑M2/76‑99 (MQ=255)
                     gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                            <  1:136499/99‑1 (MQ=255)
                      ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                           <  2:1853667‑M2/26‑1 (MQ=255)
                      ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                           <  1:1151184‑M2/26‑1 (MQ=255)
                       gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                          <  1:491278‑M2/27‑1 (MQ=255)
                       gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                          <  2:1069525‑M2/27‑1 (MQ=255)
                        tatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                                         <  2:657355‑M2/28‑1 (MQ=255)
                        tatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                                         <  2:1842105‑M2/28‑1 (MQ=255)
                        tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                                         >  2:316192/1‑99 (MQ=255)
                         atagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                                                        >  1:1137651‑M2/71‑99 (MQ=255)
                          tagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                                                       <  1:789114‑M2/30‑1 (MQ=255)
                              tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                   <  2:1146105‑M2/34‑1 (MQ=255)
                               atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                  <  1:1659981‑M2/35‑1 (MQ=255)
                                tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                 <  2:299591‑M2/36‑1 (MQ=255)
                                 gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                <  2:1867196‑M2/37‑1 (MQ=255)
                                  aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                               <  2:1297823‑M2/38‑1 (MQ=255)
                                  aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                               >  1:1497804‑M2/62‑99 (MQ=255)
                                   aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                              >  1:945176‑M2/61‑99 (MQ=255)
                                   aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                              <  1:866769‑M2/39‑1 (MQ=255)
                                    accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                             >  2:452928‑M2/60‑99 (MQ=255)
                                    accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                             <  1:1397907‑M2/40‑1 (MQ=255)
                                    accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                             >  1:583966‑M2/60‑99 (MQ=255)
                                     ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                            <  2:1064956‑M2/41‑1 (MQ=255)
                                      cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                           >  1:783305‑M2/58‑99 (MQ=255)
                                        gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                         <  2:1509358‑M2/44‑1 (MQ=255)
                                        gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                         >  2:1452838‑M2/56‑99 (MQ=255)
                                        gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                         <  1:165026‑M2/44‑1 (MQ=255)
                                        gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                         <  1:1714051‑M2/44‑1 (MQ=255)
                                         tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                        <  1:1681770‑M2/45‑1 (MQ=255)
                                            gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                     <  2:571688‑M2/48‑1 (MQ=255)
                                             tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                    >  2:201741‑M2/51‑99 (MQ=255)
                                                cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                 <  1:909617/99‑1 (MQ=255)
                                                 aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                <  1:278601‑M2/53‑1 (MQ=255)
                                                  atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                               <  1:1372068‑M2/54‑1 (MQ=255)
                                                  atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                               >  2:303922‑M2/46‑99 (MQ=255)
                                                      agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                           <  2:1233465‑M2/58‑1 (MQ=255)
                                                       gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                          >  2:964903/1‑99 (MQ=255)
                                                       gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                          >  2:1007119‑M2/41‑99 (MQ=255)
                                                       gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                          >  2:188832‑M2/41‑99 (MQ=255)
                                                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                         <  2:1189181‑M2/60‑1 (MQ=255)
                                                          gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                       <  1:1851793‑M2/62‑1 (MQ=255)
                                                           ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                      >  1:1887374‑M2/37‑99 (MQ=255)
                                                           ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                      >  1:1815744‑M2/37‑99 (MQ=255)
                                                             gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                    <  2:1240791‑M2/65‑1 (MQ=255)
                                                             ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                    <  2:1046826/99‑1 (MQ=255)
                                                              ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                   >  2:545292‑M2/34‑99 (MQ=255)
                                                               tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                  <  2:143080‑M2/67‑1 (MQ=255)
                                                               tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                  <  2:75990‑M2/67‑1 (MQ=255)
                                                                  attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                               <  1:1869944‑M2/70‑1 (MQ=255)
                                                                     gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                            <  1:616606‑M2/73‑1 (MQ=255)
                                                                          ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                       <  1:1496264‑M2/78‑1 (MQ=255)
                                                                          ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                       >  2:1931553‑M2/22‑99 (MQ=255)
                                                                             ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                    >  1:3465‑M2/19‑99 (MQ=255)
                                                                              tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATg                   <  2:1017601‑M2/82‑1 (MQ=255)
                                                                              tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATg                   <  1:957093‑M2/82‑1 (MQ=255)
                                                                               ccttctcccagaccaaCGTTTATGAAAATGAAGAAATAACGAGCATATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                  <  1:1063399‑M2/83‑1 (MQ=255)
                                                                                  tctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg               <  2:207325‑M2/86‑1 (MQ=255)
                                                                                        aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc         >  1:1022214‑M2/8‑99 (MQ=255)
                                                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa       >  2:1327692‑M2/6‑99 (MQ=255)
                                                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa       <  1:1626070‑M2/94‑1 (MQ=255)
                                                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa       <  1:1407368‑M2/94‑1 (MQ=255)
                                                                                           ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac      >  2:1217071‑M2/5‑99 (MQ=255)
                                                                                           ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac      >  1:985032‑M2/5‑99 (MQ=255)
                                                                                              aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa   <  1:961318‑M2/98‑1 (MQ=255)
                                                                                              aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa   >  1:396926‑M2/2‑99 (MQ=255)
                                                                                               cGTTTATGAAAATGAAGTAATCGCAAGCAGATGGCAACATTTGTGCTCTTACTAGGGGGATTTAACATGATTATTGGCTATGTACGCGTATCAAGAAAt  >  2:1384162/1‑99 (MQ=255)
                                                                                               cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCTTCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt  >  2:473466/1‑99 (MQ=255)
                                                                                               cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt  <  2:1619312/99‑1 (MQ=255)
                                                                                               |                                                                                                  
GGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAT  >  NC_000913/1209524‑1209717

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