New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 1207790 =10 (0.180)39 (0.790) 37/252 0.4 82.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_tpiA 1209619 = 8 (0.160)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_pae_tpiA/1209471‑1209764
                                                                                                                                                    |                                                                                                                                                 
ccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                                   <  2:256552‑M2/1‑1 (MQ=255)
ccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgcgcaatacaaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaatcaac                                                                                                                                                   <  2:749849‑M2/1‑1 (MQ=255)
      gcagatgcttttgcggcattccaccccacctgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                             <  2:877187‑M2/7‑1 (MQ=255)
          atgcttttgagccattccacccccactgaacgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgacataatgtagccgttcctgcattgagttttgccctttctccccagccaaCGTTTAtgaca                                                                                                                                         >  2:336515‑M2/139‑147 (MQ=255)
                ttgcggcattccaccctcactgaatgataactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcttcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                   <  2:623643‑M2/17‑1 (MQ=255)
                 tgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgcccctcctccccacgcaaCGTTTATGAAAATaaaga                                                                                                                                  >  1:427412‑M2/132‑149 (MQ=255)
                              ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                     <  1:113653‑M2/31‑1 (MQ=255)
                                gtcttcttcACAAGAAAAAGCGGCACTATTGCATTTCCGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                                   <  1:332977/146‑1 (MQ=255)
                                   gtcttcAGCAGAAAAAGCGGCACCAGTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACAAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                                <  2:136131/148‑1 (MQ=255)
                                    tgaatgagtactgattgctcaataccacgaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                               <  2:135146‑M2/37‑1 (MQ=255)
                                         agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAGGAAGAAATAACAAGCAAATGACATCATTc                                                                                                          >  2:287798/1‑149 (MQ=255)
                                             gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                      <  1:679707/149‑1 (MQ=255)
                                                gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctctttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                                   <  1:547405‑M2/49‑1 (MQ=255)
                                                      tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             <  2:286512‑M2/55‑1 (MQ=255)
                                                      tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             >  1:843206‑M2/95‑149 (MQ=255)
                                                      gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAACAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             <  1:68686/149‑1 (MQ=255)
                                                       cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                            <  1:739670/149‑1 (MQ=255)
                                                        aataccaccaatcataaccggtatagctatgaaaccggtgcgtgctaatctagccgtagctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                           >  1:233801‑M2/93‑149 (MQ=255)
                                                         ataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgtttctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                          <  1:932122‑M2/58‑1 (MQ=255)
                                                            ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                                       <  2:200060‑M2/61‑1 (MQ=255)
                                                                caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                   >  1:134583‑M2/85‑149 (MQ=255)
                                                                     ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTCCTTTTACCGGGGGGATTTAACAtga                                                                              >  2:821947‑M2/80‑148 (MQ=255)
                                                                     ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgaattgagttttgctccttctcccaaacccaCGTTTATGAAAATGAAGAAATAACAAGCACATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACCTGc                                                                              >  2:330742‑M2/80‑149 (MQ=255)
                                                                        accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                                           >  1:522442‑M2/77‑149 (MQ=255)
                                                                         ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                          >  1:107586‑M2/76‑149 (MQ=255)
                                                                             tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                      >  1:504031/1‑149 (MQ=255)
                                                                                  ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                 <  2:843206‑M2/83‑1 (MQ=255)
                                                                                     tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttactccttctcacaaaccaaCGTTTATAAAAATGAAGAAATAACAAGCAAATGGCATAATTCCTGCTTTTACCAGGGGGATTCAACATGATTATTGGCTATGTAcg                                                              >  1:831030‑M2/64‑149 (MQ=255)
                                                                                       aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTCCCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                            >  2:586789‑M2/62‑149 (MQ=255)
                                                                                        gggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt                                                           <  2:399908/148‑1 (MQ=255)
                                                                                                  tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGtatcaacacaa                                                 >  2:492954‑M2/51‑146 (MQ=255)
                                                                                                   cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAACATTGCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACCGTTATCAACAAatc                                                >  2:36587‑M2/50‑148 (MQ=255)
                                                                                                    ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                               >  2:758224‑M2/49‑149 (MQ=255)
                                                                                                     tcatgtggctgttgcggcattgaggtttgctgcttctccgacaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                              <  2:678016‑M2/102‑1 (MQ=255)
                                                                                                     taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                              >  1:421424‑M2/48‑149 (MQ=255)
                                                                                                      aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                             <  2:522442‑M2/103‑1 (MQ=255)
                                                                                                        tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAg                                           >  1:903435‑M2/45‑149 (MQ=255)
                                                                                                        tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAg                                           <  1:859625‑M2/105‑1 (MQ=255)
                                                                                                                   ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                >  1:954833‑M2/34‑149 (MQ=255)
                                                                                                                     gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                       >  1:714904‑M2/32‑108 (MQ=255)
                                                                                                                     gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                       <  2:714904‑M2/77‑1 (MQ=255)
                                                                                                                       attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                              >  2:423135‑M2/30‑83 (MQ=255)
                                                                                                                       attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                              <  1:423135‑M2/54‑1 (MQ=255)
                                                                                                                        ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGCTCTACCAc                           >  1:330216‑M2/29‑149 (MQ=255)
                                                                                                                          gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                         >  2:342060‑M2/27‑149 (MQ=255)
                                                                                                                           agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATCACAAGCAAATGGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTATTGGCTATGTACGCGTACCAACAAATGACCAGAACACAGAtctacacagaa                        >  1:828288‑M2/26‑144 (MQ=255)
                                                                                                                                   tcctgctcccaaaccaaCGTTTATGTAAATTCCGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAGGAGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATATACACCGTAATGCGCTg                >  2:751976‑M2/18‑149 (MQ=255)
                                                                                                                                    ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa               >  1:573855‑M2/17‑149 (MQ=255)
                                                                                                                                      ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGGTGAAc             >  1:311363‑M2/15‑149 (MQ=255)
                                                                                                                                      ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc             >  1:656219‑M2/15‑149 (MQ=255)
                                                                                                                                      ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc             >  1:546956‑M2/15‑149 (MQ=255)
                                                                                                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACCAGAAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc        >  2:460799‑M2/10‑149 (MQ=255)
                                                                                                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc        >  1:565597‑M2/10‑149 (MQ=255)
                                                                                                                                            cAAACCAACGTTTATGAAAATGGAGAAATAACAACCAAATGGCAGCATTCATGCTTTTAGCAAGGAGGTATAACATGATAATTGGCTATGTACGCGTATCAACAAATGACCAGACCACAGTCCGACAAAGTAATGAGCTGTCCTGTGCa       >  2:4486/1‑149 (MQ=255)
                                                                                                                                             aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg      >  2:940620‑M2/8‑149 (MQ=255)
                                                                                                                                                 caaCGTTTATCAAAATGAAGAAATAACAAGCAAATGTCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTACTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATg  <  1:600659‑M2/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_pae_tpiA/1209471‑1209764

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.