New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780510 (0.290)20 (0.720) 19/128 0.6 66.7% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 12 (0.430)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAG  >  NC_000913/1207725‑1207884
                                                                                |                                                                               
ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                                 >  1:266202‑M1/66‑81 (MQ=255)
  gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                               <  2:265833‑M1/81‑3 (MQ=255)
   aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                                               <  2:266182‑M1/18‑1 (MQ=255)
    atccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                                              <  2:266181‑M1/19‑1 (MQ=255)
     tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                                            >  2:266213‑M1/61‑81 (MQ=255)
         gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcat                                                                         >  1:266180‑M1/1‑72 (MQ=255)
         gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACt                                                                        <  1:265834/81‑1 (MQ=255)
         cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACt                                                                        <  1:265860‑M1/25‑1 (MQ=255)
           tggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTTGGAGAAGGAGCAAAACTc                                                                       <  1:265857‑M1/26‑1 (MQ=255)
             gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                    >  2:266220‑M1/53‑81 (MQ=255)
             gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                    <  1:265855‑M1/29‑1 (MQ=255)
              cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                                                   >  2:266179‑M1/1‑67 (MQ=255)
                    gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                                              <  1:265835‑M1/80‑20 (MQ=255)
                    gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                                              <  1:265836‑M1/80‑20 (MQ=255)
                     aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                                            <  1:265873‑M1/37‑1 (MQ=255)
                      gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                                           >  1:265862/1‑81 (MQ=255)
                       atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                                          <  1:265837‑M1/81‑24 (MQ=255)
                       atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                          >  2:265856/1‑81 (MQ=255)
                        tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                                         <  1:265838‑M1/81‑25 (MQ=255)
                        tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                         >  1:265859/1‑81 (MQ=255)
                         atCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                        >  2:266177‑M1/1‑56 (MQ=255)
                          tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                       <  1:265839/81‑1 (MQ=255)
                           gatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                       <  2:265864‑M1/42‑1 (MQ=255)
                            atgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                                     <  2:265858‑M1/44‑1 (MQ=255)
                              aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                                                   <  1:265840‑M1/81‑31 (MQ=255)
                              aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                                                   >  1:691413‑M1/1‑51 (MQ=255)
                                    tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcc                                             <  2:265841‑M1/81‑37 (MQ=255)
                                     gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                                             <  1:265865‑M1/52‑1 (MQ=255)
                                       cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttc                                          <  1:265842‑M1/81‑40 (MQ=255)
                                       cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                          <  2:265843/81‑1 (MQ=255)
                                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttca                                         <  2:265844‑M1/81‑41 (MQ=255)
                                         cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                                        <  2:265845‑M1/81‑42 (MQ=255)
                                           gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                                       >  2:266173‑M1/1‑38 (MQ=255)
                                            tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                                      <  1:265846‑M1/80‑44 (MQ=255)
                                            atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                      >  1:266205‑M1/22‑80 (MQ=255)
                                            atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                      >  2:266218‑M1/22‑80 (MQ=255)
                                              cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                   <  2:265847/81‑1 (MQ=255)
                                                 aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccg                                <  1:265849‑M1/81‑50 (MQ=255)
                                                 aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                <  2:265848/81‑1 (MQ=255)
                                                   cattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                              <  2:265866‑M1/67‑1 (MQ=255)
                                                    attttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                              >  2:1014509‑M1/14‑80 (MQ=255)
                                                       ctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa                            >  2:266171‑M1/1‑26 (MQ=255)
                                                             aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt                    <  1:265850/81‑1 (MQ=255)
                                                               ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct                  <  2:266168‑M1/81‑64 (MQ=255)
                                                               ccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                  >  2:265861/1‑81 (MQ=255)
                                                                cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct                  >  2:266167‑M1/1‑17 (MQ=255)
                                                                   GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt               >  1:266208‑M1/1‑80 (MQ=255)
                                                                    gTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaat             <  2:266166‑M1/81‑69 (MQ=255)
                                                                    gTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaat             >  1:266165‑M1/1‑13 (MQ=255)
                                                                     TTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg            <  2:265869‑M1/81‑1 (MQ=255)
                                                                      TTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt           <  2:265871‑M1/81‑1 (MQ=255)
                                                                       tGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc           <  1:266198‑M1/80‑71 (MQ=255)
                                                                        gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa         >  2:266162‑M1/1‑9 (MQ=255)
                                                                        gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa         >  2:266163‑M1/1‑9 (MQ=255)
                                                                        gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa         >  1:266164‑M1/1‑9 (MQ=255)
                                                                          GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa       >  1:266203‑M1/1‑81 (MQ=255)
                                                                            GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt     >  1:265851‑M1/1‑81 (MQ=255)
                                                                            GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt     >  2:265852‑M1/1‑81 (MQ=255)
                                                                               ggctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgcc  <  2:266161‑M1/81‑80 (MQ=255)
                                                                               GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg  >  2:266211‑M1/1‑81 (MQ=255)
                                                                                GAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg  <  2:265875‑M1/80‑1 (MQ=255)
                                                                                |                                                                               
AAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAG  >  NC_000913/1207725‑1207884

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