New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =33 (0.630)7 (0.150) 7/250 5.1 20.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 26 (0.560)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

TTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913/1209474‑1209767
                                                                                                                                                 |                                                                                                                                                    
tttGACGGATAAGCCTTTGCCAGCTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCAACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGtt                                                                                                                                                   <  1:261096/149‑1 (MQ=255)
 ttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                                  <  1:736949/149‑1 (MQ=255)
   gTCGGATAAGCCTTTGCCAGTTCCGGGTTTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                                <  1:565385/149‑1 (MQ=255)
   gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                                >  1:451712/1‑149 (MQ=255)
     cGGATAAGCCTTTGCCAGTTCCGGGTATTGTGCAGCAGAAAAAGGGGCACCATTGCATTTCAGCCAGCCTGTTGGTGGAGTGGCTGATGGCCACGGAATCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCTACCTTTATGa                                                                                                                                              <  2:799533/149‑1 (MQ=255)
         tAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCATGAAGACGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaag                                                                                                                                          >  1:360526/1‑148 (MQ=255)
              cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                                                                                                     <  1:243944/149‑1 (MQ=255)
                                      agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                             <  2:360526/149‑1 (MQ=255)
                                                    cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTGCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                               <  1:588143/149‑1 (MQ=255)
                                                      ccATTGCATTTCCCACAGGCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                             <  2:604702/149‑1 (MQ=255)
                                                      ccATCGCATTTCAGCCAGGCTGTTGGGGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                             <  2:506113/149‑1 (MQ=255)
                                                        aTTGCATTTCAGCCAGCCTGTTGGCGGAGTTGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCACATGGCACCATTCCTGCTTTTACCCggcgcg                                                                                           >  2:914274/1‑145 (MQ=255)
                                                             aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                      <  2:816074/149‑1 (MQ=255)
                                                                  agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACGAACAGGTAATGCAGAGCCTACTCCCAAACCAACGTTTGTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTACTTTTACCAGGGGGATTTGACATGc                                                                                 >  1:828648/1‑149 (MQ=255)
                                                                    ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGGtt                                                                               >  2:721266/1‑149 (MQ=255)
                                                                      agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                             >  1:423620/1‑149 (MQ=255)
                                                                               gcggcgTGGTGGAAGGCCAAGGAACCGGGAAACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                    <  2:423620/149‑1 (MQ=255)
                                                                               gcggagcggaTGAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCACAACCCAACGTTTATGAAACGGAAGAAATAACAAGCAAATGGCATCATTCGTGCTGGTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                    <  2:451712/149‑1 (MQ=255)
                                                                                       gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                            >  1:662442/1‑149 (MQ=255)
                                                                                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                     >  2:532878/1‑149 (MQ=255)
                                                                                                  cGGAAGCAGGACACCAAAAGGTAATGCAGAGCCTTCTCCCAAACCATCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                                 <  2:764997/149‑1 (MQ=255)
                                                                                                                cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCGACAAATGACCAGAACACAGATCt                                   >  1:14092/1‑149 (MQ=255)
                                                                                                                cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                   >  1:665581/1‑149 (MQ=255)
                                                                                                                   cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGACATAACAAGCAAATGTCAGCATTCCTGCTTTTACCAGGCGGATTTAACATGCTTAATGGCTATGTACGCGTCTCAACAAATGACCATAACACAGATCTACa                                <  2:828648/149‑1 (MQ=255)
                                                                                                                     ggTAATGCAGAGCCTTCTCCCAAACCAACGTATATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                              <  1:359497/149‑1 (MQ=255)
                                                                                                                           tttgctccTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg                        <  2:842456/143‑1 (MQ=255)
                                                                                                                               agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                    >  1:145820/1‑149 (MQ=255)
                                                                                                                                         cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa          >  2:137789/1‑149 (MQ=255)
                                                                                                                                              cAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATg     <  1:375393/149‑1 (MQ=255)
                                                                                                                                                aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGCAATGCGCTGAACTGTGCAGAATGCg   >  1:586631/1‑149 (MQ=255)
                                                                                                                                                 cGTTTATGCAATTGAAGAAATAACAAGCAAATGCCATCATTCCTACTTGTACCAGTGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  2:639215/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913/1209474‑1209767

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