New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =6 (0.190)7 (0.270) 7/128 1.9 58.9% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 5 (0.200)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207711‑1207868
                                                                               |                                                                              
cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                               >  2:237208‑M1/1‑81 (MQ=255)
  cgcgcTGCTGTTAAGCCGCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                                             <  2:237176‑M1/81‑1 (MQ=255)
   gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                            <  2:237177‑M1/81‑1 (MQ=255)
    cgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTT                                                                           <  2:237178‑M1/81‑1 (MQ=255)
    cgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTT                                                                           >  1:237386‑M1/1‑81 (MQ=255)
        gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                        >  2:237391‑M1/1‑80 (MQ=255)
         ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                                       <  2:237180‑M1/80‑1 (MQ=255)
         ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                                       <  2:237179‑M1/80‑1 (MQ=255)
          catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                                        <  1:237434‑M1/9‑1 (MQ=255)
            ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                    <  1:237182‑M1/80‑1 (MQ=255)
            ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                    <  1:237181‑M1/80‑1 (MQ=255)
              ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                 >  2:237491‑M1/66‑81 (MQ=255)
                 aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                              >  1:237464‑M1/63‑81 (MQ=255)
                             gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgt                                                  <  2:237183‑M1/81‑16 (MQ=255)
                             gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                  <  2:237184/81‑1 (MQ=255)
                                   aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                            >  1:237462‑M1/45‑81 (MQ=255)
                                         cAAATCGGATGGCACTGTGCAGACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                       >  2:237215/1‑80 (MQ=255)
                                           aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                    >  2:237431‑M1/1‑52 (MQ=255)
                                            aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                   <  2:237185/81‑1 (MQ=255)
                                              catttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                                 >  1:237459‑M1/34‑81 (MQ=255)
                                              cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                                 >  2:237206/1‑81 (MQ=255)
                                               atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                >  1:237465‑M1/33‑81 (MQ=255)
                                                  tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcc                             <  1:237186‑M1/81‑37 (MQ=255)
                                                     ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                          <  2:237210‑M1/55‑1 (MQ=255)
                                                         gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                       >  1:237430‑M1/1‑38 (MQ=255)
                                                         gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                      >  2:237204/1‑81 (MQ=255)
                                                          tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                      <  1:237187‑M1/80‑44 (MQ=255)
                                                          tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                      <  1:237188‑M1/80‑44 (MQ=255)
                                                             aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                  >  2:237429‑M1/1‑34 (MQ=255)
                                                              cttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                  <  1:237209‑M1/63‑1 (MQ=255)
                                                              aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc                 >  1:237427‑M1/1‑33 (MQ=255)
                                                                 ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc               >  1:237426‑M1/1‑30 (MQ=255)
                                                                    tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa            >  2:237425‑M1/1‑27 (MQ=255)
                                                                      tcGAAAACCTTGGTTTGGGAGAAGGctctgcactacctgttggtgtcccggttccgtggccttcagccactccgccaac           <  2:237189‑M1/79‑55 (MQ=255)
                                                                       cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc         >  1:237457‑M1/9‑80 (MQ=255)
                                                                        gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc       >  2:237423‑M1/1‑23 (MQ=255)
                                                                        gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc        <  1:237190/80‑1 (MQ=255)
                                                                         aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct      <  2:237422‑M1/81‑60 (MQ=255)
                                                                          aaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctg     >  1:237421‑M1/1‑21 (MQ=255)
                                                                             ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct  <  2:237420‑M1/81‑64 (MQ=255)
                                                                               |                                                                              
CGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207711‑1207868

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