New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =12 (0.450)12 (0.500) 12/252 0.9 48.4% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 15 (0.630)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGT  >  NC_000913/1209477‑1209756
                                                                                                                                              |                                                                                                                                         
ggcgcaTAAGCCTTCGCGAGTGCCGGGAGTCGCTCAGCAGAAAAAGCGGCACCATTGCATTTCCGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                                                                                           <  1:201186/138‑1 (MQ=255)
            tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                                               <  2:126631/143‑1 (MQ=255)
                  aggtCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAGTGCAGAGCCTTCTCCCAAGCCCACGTCTCTGAAAATgaagac                                                                                                                         >  2:264516/4‑142 (MQ=255)
                  aGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATggagaa                                                                                                                         >  2:304555/1‑143 (MQ=255)
                     cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttctgctacttatcacaaaccgaCATTTAGGAAAATGAAGAAATa                                                                                                                      >  1:87952‑M2/122‑143 (MQ=255)
                       cccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                    <  1:316082‑M2/24‑1 (MQ=255)
                        gggTATTCTTCAGCAGAAAAAACGGAAACATTGCATTAAAGCCAACCTGTAGGCAGAGTGGCTGAAGACAACGAAACCGTGACAACAACAGGTAAAGCAGAGCAGTCTACAAAATCAGCGAAGATGAAAAAGAAGAAAGAAcc                                                                                                                   >  2:148221/1‑142 (MQ=255)
                        ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagtttcgctccttctcccaccccaaCGTCTATGAAAACGAAGAAATAACa                                                                                                                   >  2:151650‑M2/119‑143 (MQ=255)
                         cccactgaatgattactgattgctcaatacccccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagcagttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                  <  1:459881‑M2/26‑1 (MQ=255)
                          ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                                 >  1:485624‑M2/117‑143 (MQ=255)
                             ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                                                                                              >  2:35131/1‑143 (MQ=255)
                                ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGGAAATg                                                                                                           >  2:87587/1‑143 (MQ=255)
                                       gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAACATGAAGAAATAACAAGCACATGGcatcag                                                                                                    >  1:261288/1‑142 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                           >  1:19488‑M2/95‑143 (MQ=255)
                                                     aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                      <  1:371753/143‑1 (MQ=255)
                                                        accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAACGGCATCATTCCTGCTCTTACCAggg                                                                                   >  2:250917‑M2/87‑143 (MQ=255)
                                                          aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                 <  1:95980/143‑1 (MQ=255)
                                                            atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                               >  2:198929‑M2/83‑143 (MQ=255)
                                                                            gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAACATTCCTGCTTTTACCAGGAGAATTTAACAAGCTCATAgg                                                               >  2:384337/1‑143 (MQ=255)
                                                                                      tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCCGAGCCTTCTCCCAAACCAACGTTAATGAAAATGACGAAATAACAAGCACATGGCATCACTCGAGCATTTACCAGGGGACTTTCACATGATTATATGCTATGTAcgc                                                     >  2:143624/1‑143 (MQ=255)
                                                                                        gtcctcctaatgtagcctttgctgctttgagttttgctccttctccccaaccaaCGTTTATGAATATGTAGAAATAACTATCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt                                                   <  2:139154‑M2/89‑1 (MQ=255)
                                                                                            tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATCCCAACTTCTACCAGGGCGCTCCAACATGCCTACTCGCTATGTGCGCGTATca                                               >  2:98646‑M2/51‑143 (MQ=255)
                                                                                              aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCATCGTTTATGAAAATGAAGAAATTACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                             <  2:480122/143‑1 (MQ=255)
                                                                                                 acgtagccgttgctgcattgagttttgctgcttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                          <  2:87952‑M2/98‑1 (MQ=255)
                                                                                                     cGGGACACCAACAGGTAATGCAGAGCCTTCGCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCCAATGGCATCATACCTGATTATGCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                      <  2:215363/143‑1 (MQ=255)
                                                                                                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                          <  1:357309‑M2/114‑1 (MQ=255)
                                                                                                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCCACAAATGACCAGAACACAGATCt                        >  1:452233/1‑143 (MQ=255)
                                                                                                                              ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg             >  1:293937‑M2/17‑143 (MQ=255)
                                                                                                                                  ctcccaacccaaCGTTTATGAACATGAAGAAACAACAAGCAAGTGGCATCATCCCTGCTTTTACCAGGGGGATTTAACATTCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt         <  1:230008‑M2/131‑1 (MQ=255)
                                                                                                                                    cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACATAATGCGCTGa       >  1:369512‑M2/11‑143 (MQ=255)
                                                                                                                                     ccgggccaaCGTTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa      <  2:418625‑M2/134‑1 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACACGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATATAACATGCTTATTGGCTATGTACGCCTATAAACAAATGACCAGCACCCAGATCTACAAAGTACTGCGCTgac      >  2:52108‑M2/10‑142 (MQ=255)
                                                                                                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt  <  2:351939‑M2/138‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                         
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGT  >  NC_000913/1209477‑1209756

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