New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780563 (0.900)14 (0.240) 14/166 5.0 22.1% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 46 (0.790)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

CTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGT  >  NC_000913/1207708‑1207903
                                                                                                 |                                                                                                  
cTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                  >  2:1703677/1‑100 (MQ=255)
 ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                 <  1:876617/100‑1 (MQ=255)
 ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                 <  1:179492/100‑1 (MQ=255)
    gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                                                                              >  2:813009/1‑100 (MQ=255)
     cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgca                                                                                             <  1:192668/100‑8 (MQ=255)
       cgcTGCTGTGAAGCCTCAAGAATCCGTTTTTCGATATCAAATCTTATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACt                                                                                           <  2:992811/100‑1 (MQ=18)
         ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                                                          <  2:927497/99‑1 (MQ=255)
         ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                         >  2:153548/1‑100 (MQ=255)
         ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                         >  1:796577/1‑100 (MQ=255)
           gatgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                                       <  1:454642/98‑1 (MQ=255)
             tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                     <  2:1141200/100‑1 (MQ=255)
              gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                    >  1:1166701/1‑100 (MQ=255)
                    ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGTAACg                                                                              >  2:46044/1‑100 (MQ=255)
                    ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                                                              <  1:1592276/100‑1 (MQ=255)
                     cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                             >  2:1383428/1‑100 (MQ=255)
                     cGCAAGAATCCGTTTGGCAATATCAAATAGGATTGCACTGTGCAAACGGCTCTCGAAAACCTTGGGGTGCGAGCAGGAGCAAAACTCAATGCAGCAACgg                                                                             >  1:83719/1‑100 (MQ=21)
                          gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                                                        >  1:136749/1‑100 (MQ=255)
                           aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                       >  2:1116566/1‑100 (MQ=255)
                           aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                       <  2:1482963/100‑1 (MQ=255)
                            aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                                      >  2:786629/1‑100 (MQ=255)
                             tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                     <  1:294877/100‑1 (MQ=255)
                               ccaccggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                                                   >  1:1675858/6‑100 (MQ=255)
                                gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACTGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                                                  <  1:782311/100‑1 (MQ=255)
                                 tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                                                 <  2:394434/100‑1 (MQ=255)
                                    ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                              >  1:51007/1‑100 (MQ=255)
                                      cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                                            >  1:66173/1‑100 (MQ=255)
                                      cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                                            >  1:874988/1‑100 (MQ=255)
                                       gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAccc                                                           <  1:872278/100‑2 (MQ=255)
                                            cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                                      <  1:458368/100‑1 (MQ=255)
                                            cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                                      <  1:229090/100‑1 (MQ=255)
                                                 cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                 <  2:236651/100‑1 (MQ=255)
                                                       gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc                                           >  1:1283759/1‑100 (MQ=255)
                                                           tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt                                       >  2:270728/1‑100 (MQ=255)
                                                               cAAACGGCTGTCGAAAAGCTTTGTTGGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGAACCGGTTTCATACCTATACAGGTTATGa                                   >  1:313045/1‑100 (MQ=255)
                                                               cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAACGGTTATGa                                   >  1:1298615/1‑100 (MQ=255)
                                                                 aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                 <  2:1283765/100‑1 (MQ=255)
                                                                   tcattttcataaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                               >  2:1075947/11‑100 (MQ=38)
                                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                               >  1:707255/1‑100 (MQ=255)
                                                                    cattttcataaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                              <  2:1332437/91‑1 (MQ=38)
                                                                     gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtn                             <  2:1049276/100‑2 (MQ=255)
                                                                     gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                             <  1:813007/100‑1 (MQ=255)
                                                                     gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                             >  2:765966/1‑100 (MQ=255)
                                                                     gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                             <  1:277776/100‑1 (MQ=255)
                                                                      ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCCACGGCTACATTAGGACGAACCGGTTTCATAGCTATACCGGTTATGATtggtg                             <  1:762612/99‑1 (MQ=255)
                                                                          cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                        <  1:1703669/100‑1 (MQ=255)
                                                                           gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                        >  1:293698/1‑99 (MQ=255)
                                                                           gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAAGATTGGTGGTATTg                       >  1:1159959/1‑100 (MQ=255)
                                                                            taaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                      >  2:1370404/2‑100 (MQ=255)
                                                                              aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGc                    <  1:1456484/100‑1 (MQ=255)
                                                                                cgttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                  >  2:132223/3‑100 (MQ=255)
                                                                                ccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                  >  1:778628/1‑100 (MQ=255)
                                                                                ccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                  >  2:1371542/1‑100 (MQ=255)
                                                                                 gttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                 >  2:1573112/2‑100 (MQ=255)
                                                                                 cTTGGTTTGGGAGAAGGTGCAAAACTCAATGCAGCAACGGCTACCTTAGGACGCACCGGTTTCATAGCGAGACCGGTTATGATTGGTGGTATTGAGCaa                  >  1:1023374/1‑99 (MQ=255)
                                                                                    ggTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg              >  1:171784/1‑100 (MQ=255)
                                                                                     ntttGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAATAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt             >  1:929369/2‑100 (MQ=255)
                                                                                      tttGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt             <  1:270727/99‑1 (MQ=255)
                                                                                      tttGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTa            >  1:1567393/1‑100 (MQ=255)
                                                                                      tttGGGAGAAAGAGCAAAACTCAATGCAGCAACGGCTACATTAGGCCGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTAAGCAATCAGTa            <  1:440267/100‑1 (MQ=255)
                                                                                        tGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt          <  2:484526/100‑1 (MQ=255)
                                                                                         gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGATATGATTGGTGGTATTGAGCAATCAGTAAt          >  1:660913/1‑99 (MQ=255)
                                                                                              aaGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCa    <  2:708990/100‑1 (MQ=255)
                                                                                                ggAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt  >  2:1281345/1‑100 (MQ=255)
                                                                                                ggAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt  >  2:691676/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGT  >  NC_000913/1207708‑1207903

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.