New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =7 (0.150)35 (0.850) 32/240 NT 78.5% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 13 (0.320)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_bme_pgi/1209476‑1209765
                                                                                                                                               |                                                                                                                                                  
cgTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGTCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                                                                                                                               <  1:307880/148‑1 (MQ=255)
       tAAGCCTTTGCCAGTTTCGGGTTTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                                                        <  1:690605/149‑1 (MQ=255)
         ttttgcgggattccacccccactgaatgattactgattgctcaatatcaccaatcataaccggtatagctatgaaacaggtgcgtcctaacgcagccggtgcttcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                      <  1:397660‑M2/15‑1 (MQ=255)
            cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGAAGTGGCTGAAGGCCACGGAACAGAGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                                                                                                   >  1:70493/1‑149 (MQ=255)
              cggcattcgaccacctttgaatgattaatgatttctcaataccattaatcataaccggtatagctatgaaaacgttgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccacaccaaCGTTTATGAAAATGAAGaaa                                                                                                                                 <  1:42750‑M2/20‑1 (MQ=255)
                   ttccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatcgctatgaaaccggtgcgtcctaatgtagccgttgcttcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                                                                                                            >  2:217223‑M2/125‑149 (MQ=255)
                   ttccaccaccactgaatgattactgattgcgcaaaaccaccaatcataaacggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                                                                                                            >  2:305486‑M2/125‑149 (MQ=255)
                    gatgcgtgTATTATACAGCAGAAATAGCGGCACCATTGCATTTCAGCCAGCCTTTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                           <  2:349915/142‑1 (MQ=255)
                     ccagccccactcaatggttactgtttgcttaataccaccaatcataaccggtatagctatgaaatcggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                                          <  2:139364‑M2/27‑1 (MQ=255)
                      cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAATTAACAaga                                                                                                                         >  2:128063‑M2/122‑148 (MQ=255)
                       acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                                        >  2:629278‑M2/121‑149 (MQ=255)
                            cactgaatgattacggattggtcaatagcaccaatcataaccggtatagctatgaaaccggcgcggcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAACTgg                                                                                                                   <  2:459072‑M2/34‑1 (MQ=255)
                               tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtggggcctaatgtagccgttcctgcattgagttttgctacttctccgaaaccaaCGTTTATGAAAATGAAAAAATAACAAACAAATGGcat                                                                                                                >  2:327104‑M2/113‑149 (MQ=255)
                               tcttcAGCAGAAAAATCGGCTCCCTTGCATTTGAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGACCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                                <  1:585793/149‑1 (MQ=255)
                                      ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtaaccgttgctgcattgagttttgctccttctcccaaatcaaCGTTTATAAAAATGAAGAAATAACAACCAAATGGCATCATacca                                                                                                         >  2:672763‑M2/106‑148 (MQ=255)
                                          tgattgctcaataccaccaatcataaccggtatagctatgaaacaggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctaccaaaccaaCATTTATGAAAATGAAAAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                                     >  1:190631‑M2/102‑149 (MQ=255)
                                            attgctcaataccaccaatcataaccggtatagctatgaaaccggagcgtcctaatgtagccgttgctgcattgagttgtgctccttctcccaaagcaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTgattgt                                                                                                   >  2:813002‑M2/100‑147 (MQ=255)
                                                ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                               <  2:703213‑M2/54‑1 (MQ=255)
                                                  tgataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                             <  2:218927‑M2/56‑1 (MQ=255)
                                                  caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGATGAAATAACAAGCAAATGGCATCATTACTGCTTTTACCAgg                                                                                             >  1:633871‑M2/94‑149 (MQ=255)
                                                    ataccaccaaacataaccggtatagctatgaaaccggtgcgtcctaatgtagacgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                           <  1:844018‑M2/58‑1 (MQ=255)
                                                     cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACAACAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                          >  1:695667/1‑149 (MQ=255)
                                                        caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                                    >  2:345300‑M2/88‑88 (MQ=255)
                                                        caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                                    <  1:345300‑M2/1‑1 (MQ=255)
                                                        caccaatcataaccggtatagctatgaaaccggtacgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAGGAAGAAATAACGAGCAAATCGCATCATTCCTGCGTTTACCAGGGggact                                                                                       >  2:10859‑M2/88‑147 (MQ=255)
                                                          ccaatcataaccggtatagatatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATAACAATCAAATGGCATCATTCCTGCTTTTACCACGGGGatttt                                                                                     >  2:689620‑M2/86‑148 (MQ=255)
                                                           caataataaccggtagagctatgaaacgggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                    <  2:310178‑M2/65‑1 (MQ=255)
                                                             ttCAGCCAGTCTGTTGGCGGAGTGGCTTAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                  <  1:266621/149‑1 (MQ=255)
                                                             atcataaccggtatagttaagcaaccggtgcgtcctcatgtagccgttggtgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAATTAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                  <  2:542422‑M2/67‑1 (MQ=255)
                                                             atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                  <  1:503376‑M2/67‑1 (MQ=255)
                                                             atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                  <  1:449976‑M2/67‑1 (MQ=255)
                                                              tcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                                 <  1:751078‑M2/68‑1 (MQ=255)
                                                                aaccacccTGTTGTCGGAGTGGCTGAAGACCACGGAACCGGTACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                                               <  1:422832/147‑1 (MQ=255)
                                                                 taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                                                              >  1:453446‑M2/79‑149 (MQ=255)
                                                                   accgggatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                                            <  2:751456‑M2/73‑1 (MQ=255)
                                                                        tatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                       <  2:856891‑M2/78‑1 (MQ=255)
                                                                              tactaaaccggtgcgtcctaatgtcgccgtggctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                                 <  2:788560‑M2/84‑1 (MQ=255)
                                                                                tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAACATGAAGAAATAACAAGCAAATGGCATAATTCCTGCCTTTACCAGGGGGATTTAAACTCCTTATTACCTATGTAcg                                                               >  2:747034‑M2/64‑149 (MQ=255)
                                                                                  aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccaaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAtgcg                                                             >  2:704305‑M2/62‑149 (MQ=255)
                                                                                   tGGCTGAAGGCCACGGAACCGGGACAACAACAGGTAATGCAGAGCATTCTCCCAAATCAATGTTTATGAAAATGAAGAAATAACAAGCAAGTGGCATAATTCCTGCTTTTCCCAGGGGGATTTAACATGATTGTTGGGTGTGTGTGCGt                                                            >  2:245972/1‑149 (MQ=255)
                                                                                      cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTACTGCTTTTACAAGGAGGATTTAACATACTTTTTGGCTATGTCCGGATATc                                                         >  1:413472‑M2/58‑149 (MQ=255)
                                                                                        gAAGGCCACGGAATCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAGTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                       <  1:32503/149‑1 (MQ=255)
                                                                                                 aaagtagccgttgctgcattgaggtttactccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                              <  2:703678‑M2/103‑1 (MQ=255)
                                                                                                     tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                          <  1:502372‑M2/107‑1 (MQ=255)
                                                                                                     ccGGGACAACAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                          <  2:695667/149‑1 (MQ=255)
                                                                                                            ggctgcattcagttttgtcacttctaccaaaccaaCGTTTATTAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                   <  2:349034‑M2/114‑1 (MQ=255)
                                                                                                                  attgagttttgctccttctcccaaaccaaCGTTTATGCAGATGAAGAAATTACAAGCAAATTGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                             <  1:244851‑M2/120‑1 (MQ=255)
                                                                                                                   ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                            <  1:22107/149‑1 (MQ=255)
                                                                                                                        tGCAGAGCGTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGAACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                       >  2:795424/1‑149 (MQ=255)
                                                                                                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg                      >  2:225783‑M2/23‑149 (MQ=255)
                                                                                                                          ttgctctttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgc                     <  2:781804‑M2/128‑1 (MQ=255)
                                                                                                                             ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                  >  2:225586‑M2/19‑149 (MQ=255)
                                                                                                                                cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGTGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa               >  1:332784‑M2/16‑149 (MQ=255)
                                                                                                                                cttctcccaaacacaCGTTTATTAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa               <  2:119893‑M2/134‑1 (MQ=255)
                                                                                                                                 ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGTATCATTCCTGCTTTAACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAACGACCAGAACACAGATCTACAACGTAATGCGCTGAAc              >  1:791009‑M2/15‑149 (MQ=255)
                                                                                                                                    tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAATTAACATGCTTATTGGCTATGTACGCGTATGAACAAATGACAAAAACACATATCTACAAGGTAATACCCTGAActgg           >  2:285463‑M2/12‑148 (MQ=255)
                                                                                                                                      ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATGTACCACACAATGCGCTGAAcggtga         >  2:588328‑M2/10‑148 (MQ=255)
                                                                                                                                           ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGAACTAACATGCTTATTGGCTATGTCCGCGTATCAACAAATGACCAAACCACAGATCTACAACGCAATGCGCTGAACTGTGCAGGAt    >  2:398449‑M2/5‑149 (MQ=255)
                                                                                                                                             aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAAATGCTTATTGGCTATGTACGCGTATCAACAAATGACCACAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGCTGc  >  2:398101‑M2/3‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_bme_pgi/1209476‑1209765

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.