New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780542 (0.640)18 (0.330) 17/166 3.7 34.6% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 33 (0.600)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.

ACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCA  >  NC_000913/1207707‑1207901
                                                                                                  |                                                                                                
aCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggc                                                                                                 <  2:1294778/100‑2 (MQ=255)
 cTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                >  2:1306488/1‑100 (MQ=255)
 cTTCGCGCGCTGCTGTTAAGCCGCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                <  1:156703/100‑1 (MQ=255)
   tCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCa                                                                                              >  2:1479274/1‑100 (MQ=255)
     gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                                                                            >  2:169747/1‑100 (MQ=255)
      cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgca                                                                                           <  1:518040/100‑8 (MQ=255)
         gctggtGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                                                        <  1:1118863/100‑1 (MQ=255)
              tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                   <  1:1479259/100‑1 (MQ=255)
               gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                  <  1:1138993/100‑1 (MQ=255)
                 tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                                <  1:169352/100‑1 (MQ=255)
                  agggcGCAAGAATCCGGTTGGCGAGATCAAATGGGATGGCACTGTGCAAACGGCTCTCGAAAACCGTGGTTTGGGAGAAGGAGCAAAACTCCATgcag                                                                                 <  1:243588/94‑1 (MQ=18)
                   aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                                                              <  1:1245408/100‑1 (MQ=255)
                      cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                           >  2:856141/1‑100 (MQ=255)
                          aGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                                                       <  1:1547285/100‑1 (MQ=255)
                          aGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                                                       <  1:847803/100‑1 (MQ=255)
                             aTCCGTTTGGTGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                                    >  1:255095/1‑100 (MQ=255)
                                cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                                                 <  1:774832/100‑1 (MQ=255)
                                  tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                                               <  1:794499/100‑1 (MQ=255)
                                     ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                            >  2:366870/1‑100 (MQ=255)
                                       cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGGAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                                          >  1:1021489/1‑100 (MQ=255)
                                         atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                                        >  2:829871/1‑100 (MQ=255)
                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                     >  2:377612/1‑100 (MQ=255)
                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                     <  1:936704/100‑1 (MQ=255)
                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                     >  2:377472/1‑100 (MQ=255)
                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                     <  2:750208/100‑1 (MQ=255)
                                                             gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                                    >  2:558108/1‑100 (MQ=255)
                                                                cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                 >  2:1128560/1‑100 (MQ=255)
                                                                  aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                               >  1:536504/1‑100 (MQ=255)
                                                                  aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                               <  1:1029981/100‑1 (MQ=255)
                                                                   aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTTGGACGCACCGGTTTCATAGCTATACCGGTTATGATtg                              <  1:856134/100‑1 (MQ=255)
                                                                      gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                           >  1:300438/1‑100 (MQ=255)
                                                                           cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCTACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                      >  2:1080741/1‑100 (MQ=255)
                                                                           cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                      >  1:251454/1‑100 (MQ=255)
                                                                             taaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                    >  2:1177596/2‑100 (MQ=255)
                                                                               aaCCTTGGTTTGGGAGAAGGAGCAAACCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGc                  <  2:317451/100‑1 (MQ=255)
                                                                                      gTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt           >  2:310931/1‑100 (MQ=255)
                                                                                      gTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt           >  2:564674/1‑100 (MQ=255)
                                                                                       tttGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTa          >  1:1527909/1‑100 (MQ=255)
                                                                                        ttGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTaa         <  2:945855/100‑1 (MQ=255)
                                                                                          gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc       >  2:715110/1‑100 (MQ=255)
                                                                                          gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc       >  1:1308372/1‑100 (MQ=255)
                                                                                            gagaAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAt     >  1:599003/1‑100 (MQ=255)
                                                                                                aGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCa  >  1:427280/1‑99 (MQ=255)
                                                                                                 ggAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCa  >  2:431865/1‑98 (MQ=255)
                                                                                                  |                                                                                                
ACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCA  >  NC_000913/1207707‑1207901

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