New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078056 (0.190)15 (0.590) 14/128 0.8 75.4% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 5 (0.200)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGC  >  NC_000913/1207725‑1207885
                                                                                |                                                                                
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)
                                                                  TGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                <  2:237205‑M1/81‑1 (MQ=255)
                                                                   ggTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaa               >  1:237419‑M1/1‑14 (MQ=255)
                                                                   GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                >  2:237452‑M1/1‑80 (MQ=255)
                                                                    gTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaat              <  1:237461‑M1/81‑69 (MQ=255)
                                                                       tGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc            <  1:237418‑M1/80‑71 (MQ=255)
                                                                         ggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaat         <  2:237417‑M1/81‑74 (MQ=255)
                                                                         GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg          >  2:237191‑M1/1‑80 (MQ=255)
                                                                            gaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccaccccgccaacaggctggctgaaatgcaatggt      >  1:237416‑M1/1‑5 (MQ=255)
                                                                              aGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc    >  1:237415‑M1/1‑3 (MQ=255)
                                                                               GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg   <  1:237221‑M1/81‑1 (MQ=255)
                                                                                gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccg  <  1:237414‑M1/81‑81 (MQ=255)
                                                                                |                                                                                
AAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGC  >  NC_000913/1207725‑1207885

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