New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =7 (0.240)17 (0.670) 16/252 0.6 75.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 5 (0.200)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGG  >  NC_000913/1209483‑1209761
                                                                                                                                        |                                                                                                                                              
tAAGCCGTGGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                          <  1:32300/143‑1 (MQ=255)
             tccatccccactgcatgattactgattgctcaccaccaccattcatacccggtatagctatgaaaccgttgcgtcctaatgtagccgttgctgcattgagtttagctcctgctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                             <  2:520523‑M2/20‑1 (MQ=255)
                  ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgtcgctgcatggagttttgctccttcccccaaaccaaCGTTTATCAAAATGAAGAAATACCa                                                                                                                        >  2:473945‑M2/119‑143 (MQ=255)
                   gttATTCCTCAGCATAAAAGGCGGGGCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                       <  2:317617/141‑1 (MQ=255)
                           cggacttctgtttgcgcaatagcatcaatcataaccggtatagctatgaaaccggtgcgtcctaatgtatccgttgctgcattgagttgtgctccttctcccaaaccaaCTTTTATGAAAATGAAGAAAGAACAAGCCAATgg                                                                                                               <  2:333970‑M2/34‑1 (MQ=255)
                             gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                             >  1:77982‑M2/108‑143 (MQ=255)
                                       tgctcaataccaccaatcagaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                   <  1:8146‑M2/46‑1 (MQ=255)
                                         ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagctctgctccttctcccaaaccaaCGTTTATGAAAATGAAGGAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                                 >  2:383865‑M2/96‑143 (MQ=255)
                                          tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGATTATGAAAATGAAGAAATAACAAGCAAATGGAATCATTACTgctat                                                                                                >  2:284822‑M2/95‑141 (MQ=255)
                                                   ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                       >  1:464489‑M2/86‑143 (MQ=255)
                                                   ccaatcataaccggtatagctatgaaaccggtgcgtactaatgtagccgttgctgaattgagttttgctccttctcccaaaccaaCGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTCCTTTTACCAgggg                                                                                       >  1:339448‑M2/86‑143 (MQ=255)
                                                         ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAAACCTTCTTTTACAAGGGGGATTTa                                                                                 >  2:339264‑M2/80‑143 (MQ=255)
                                                           aaccggtatagctatgaaaccggagcgtcctaatgtagccgttgctgcattgagttttgctccttctaacaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAGTCCTGCTGTTACCAGGGGGGTTtaag                                                                               >  2:426844‑M2/78‑142 (MQ=255)
                                                                            aaccggtgcgtcctaatgtagccgttgctgcattaaattttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTACTATTACCAGGGGGATTTAACATACTTATTAGCtatat                                                              >  2:204931‑M2/61‑141 (MQ=255)
                                                                              gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGcgacgaac                                                            >  2:244676/1‑136 (MQ=255)
                                                                                ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGAGAATTAACAAACTTATTGGAAATGTAagc                                                          >  1:917‑M2/57‑143 (MQ=255)
                                                                                     gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTGCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                     >  1:365831‑M2/52‑143 (MQ=255)
                                                                                                      tctgcattgagttttgctcctactcccaaaccaaCGTCTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTACCAGGGGGAGTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                    <  2:161505‑M2/109‑1 (MQ=255)
                                                                                                        tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGATTTAACATGCTTATTGGCTATGTACGCGTAGCAACAAATGACCAGAAcaca                                  >  2:520401‑M2/33‑143 (MQ=255)
                                                                                                        aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTCATTGGCTACGTACGCGTATCAACAAATGACCAGAAcaca                                  >  2:270386/1‑143 (MQ=255)
                                                                                                        aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATATAACATGCTTATTGGCTATGTACACGTATCAGCAAATGACCAGAAcaca                                  >  2:429722/1‑143 (MQ=255)
                                                                                                            ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTCACGAGACCGAGGGAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACACGACCACAGATc                              >  1:119218‑M2/29‑143 (MQ=255)
                                                                                                                               cccaaccaaCCGTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAACGACCAGCACACAGATCTACAACGTAATGCGCTGaa           <  2:104584‑M2/134‑1 (MQ=255)
                                                                                                                               ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa           >  2:92394‑M2/10‑143 (MQ=255)
                                                                                                                                        cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAgg  <  1:244676/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                              
TAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGG  >  NC_000913/1209483‑1209761

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