New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =16 (0.360)12 (0.300) 8/244 NT 42.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 18 (0.450)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGT  >  NC_000913/1209481‑1209756
                                                                                                                                          |                                                                                                                                         
atgcttttgcggcattccacccccactgaatgattattgattgctcaataccaccaatcaaaaccggaatagctgagcaaacggtgcgtcctaatgcagacgttgctgcattgcgttttgctcctgctcccaaaacaaCGTTTAGaaaa                                                                                                                                >  2:47551‑M2/139‑149 (MQ=255)
atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTCTGgaa                                                                                                                                 >  2:689766‑M2/139‑149 (MQ=255)
atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcagtgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                 >  1:703661‑M2/139‑149 (MQ=255)
atgcttttgcggcattccacccccactgaatgattactgatggctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                 <  1:175294‑M2/11‑1 (MQ=255)
 aTAAGCCTTTGCCAGTTCCGGGTATGCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAGGaaaa                                                                                                                                >  2:606472/1‑149 (MQ=255)
         ttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                                                                                        >  1:122829/1‑149 (MQ=255)
            ccAGTTCCGGGTATTCTTGAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                                     <  2:152518/149‑1 (MQ=255)
            cattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                                                                     >  2:516770‑M2/127‑149 (MQ=255)
               tcgccccccacggaattattactgattgctcaataccaccaatcataaccggcatagctatgaaaccggtgcgtcctaatatagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                  <  2:537597‑M2/26‑1 (MQ=255)
                            ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAACACATGGcaaca                                                                                                     >  2:535224/1‑149 (MQ=255)
                                    ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATGCCTGCt                                                                                             >  1:506996‑M2/103‑149 (MQ=255)
                                        ttgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAAAAGCAAATGGCATCATTCCTGCTTTTa                                                                                         >  1:759412‑M2/99‑149 (MQ=255)
                                          gctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTGTGAAAATGAAGAAATACCAAGCAAATGGCAGCATTCATGCTTTTAcc                                                                                       >  2:445420‑M2/97‑149 (MQ=255)
                                                  cctccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAAGGGGAt                                                                               >  2:99772‑M2/89‑149 (MQ=255)
                                                      aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                           >  2:747948/1‑149 (MQ=255)
                                                       tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                          <  2:535344/149‑1 (MQ=255)
                                                        atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                         <  2:49202‑M2/67‑1 (MQ=255)
                                                        atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                         <  1:189377‑M2/67‑1 (MQ=255)
                                                                   tggtGGCGGAGTGGCCGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                              <  1:290873/146‑1 (MQ=255)
                                                                   tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACCGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                              <  1:64539/149‑1 (MQ=255)
                                                                    gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                             >  1:126203/1‑149 (MQ=255)
                                                                         cGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTGCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTaggac                                                        >  2:512775/1‑145 (MQ=255)
                                                                                     aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca                                            >  1:62997/1‑149 (MQ=255)
                                                                                            ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                     <  1:718336/149‑1 (MQ=255)
                                                                                                   ggACACCAACAGGTAATGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAACCTTGAATGACCAGAacac                              >  1:217259/1‑149 (MQ=255)
                                                                                                             aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                    <  1:714586/149‑1 (MQ=255)
                                                                                                                tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                 >  1:3177/1‑149 (MQ=255)
                                                                                                                tAATGCAGAGCCTTCTCCCAAACCAACGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGGACACAGATCTACAACGt                 >  1:3729/1‑149 (MQ=255)
                                                                                                                  aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa               >  1:611738/1‑149 (MQ=255)
                                                                                                                  aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAAACAGATCTACAACGTaa               >  2:6358/1‑149 (MQ=255)
                                                                                                                        ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGGTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCCACAAATGGCTAGAACACAGATCTACAACGTAATGCGCt         >  1:252647‑M2/19‑149 (MQ=255)
                                                                                                                        agCCTTCTCCCAAACCAACGTTTATGAAAATCAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACACGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt         <  1:535224/149‑1 (MQ=255)
                                                                                                                              ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtg   >  1:312446‑M2/13‑149 (MQ=255)
                                                                                                                               tcccaaaccagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt  <  2:551470‑M2/138‑1 (MQ=255)
                                                                                                                               tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt  <  1:747948‑M2/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGT  >  NC_000913/1209481‑1209756

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.