New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3363001 =40 (1.000)34 (0.850) 25/274 0.2 62.5% coding (195/2382 nt) yhcD putative fimbrial usher protein YhcD
?NC_000913 3766087 = 1 (0.030)coding (3905/4134 nt) rhsA rhs element protein RhsA

TTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTC  >  NC_000913/3765948‑3766212
                                                                                                                                           |                                                                                                                             
ttGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCAtt                                                                                                                               >  2:629649/1‑140 (MQ=255)
    cAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGgg                                                                                                                           >  1:206873‑M2/1‑138 (MQ=255)
      atagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGatat                                                                                                                         <  1:350774‑M2/7‑1 (MQ=255)
       ttCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaag                                                                                                                        <  1:664501‑M2/140‑6 (MQ=255)
        agtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATaa                                                                                                                       <  1:96262‑M2/9‑1 (MQ=255)
        agtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATaa                                                                                                                       >  1:105788‑M2/132‑140 (MQ=255)
         gtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAAc                                                                                                                      >  1:567846‑M2/131‑140 (MQ=255)
          tcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGatat                                                                                                                         >  1:418855‑M2/130‑136 (MQ=255)
          tcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGatat                                                                                                                         <  2:418855‑M2/7‑1 (MQ=255)
           aTGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgc                                                                                                                    >  1:77569‑M2/1‑131 (MQ=255)
            tGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcg                                                                                                                   <  1:221707‑M2/140‑11 (MQ=255)
            tGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcg                                                                                                                   <  2:477826‑M2/140‑11 (MQ=255)
             taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                  >  2:125313‑M2/127‑140 (MQ=255)
             taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                  >  1:284136‑M2/127‑140 (MQ=255)
                   ccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCt                                                                                                            >  2:546976‑M2/121‑140 (MQ=255)
                     aGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                          >  2:664501‑M2/1‑121 (MQ=255)
                      gCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccc                                                                                                         >  1:200830‑M2/1‑120 (MQ=255)
                           gCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgata                                                                                                    >  2:221707‑M2/1‑115 (MQ=255)
                             cttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                  >  1:529541‑M2/111‑140 (MQ=255)
                             cttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                  >  1:464157‑M2/111‑140 (MQ=255)
                             cttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                  <  1:257011‑M2/30‑1 (MQ=255)
                              cagcagCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatga                                                                                                 >  1:413675‑M2/1‑112 (MQ=255)
                                   aaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTg                                                                                            >  1:663101‑M2/105‑140 (MQ=255)
                                     ggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCa                                                                                          <  1:546976‑M2/38‑1 (MQ=255)
                                        caGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgttt                                                                                       >  2:177159‑M2/1‑102 (MQ=255)
                                          ttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGacac                                                                                     <  2:375691‑M2/43‑1 (MQ=255)
                                          ttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGacac                                                                                     <  2:567846‑M2/43‑1 (MQ=255)
                                            taaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgcctttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAg                                                                                   >  1:668889‑M2/96‑140 (MQ=255)
                                             taaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                  >  2:490916‑M2/95‑124 (MQ=255)
                                             taaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                  <  1:490916‑M2/30‑1 (MQ=255)
                                             taaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAgg                                                                                  <  2:105788‑M2/46‑1 (MQ=255)
                                                 aTGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgg                                                                              >  1:629649‑M2/1‑93 (MQ=255)
                                                 aTGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgg                                                                              <  1:66255‑M2/140‑48 (MQ=255)
                                                     tgcAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcc                                                                          <  1:121751‑M2/140‑52 (MQ=255)
                                                      ttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGa                                                                         <  2:533299‑M2/55‑1 (MQ=255)
                                                         aaaGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatag                                                                      <  1:498762‑M2/140‑56 (MQ=255)
                                                          tgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCgt                                                                     <  1:393819‑M2/59‑1 (MQ=255)
                                                          tgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCgt                                                                     <  2:590526‑M2/59‑1 (MQ=255)
                                                             cATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                                                                  >  1:258270‑M2/1‑81 (MQ=255)
                                                                 gACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                                                              <  2:660702‑M2/140‑64 (MQ=255)
                                                                  aCAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                                                             >  1:392524‑M2/1‑76 (MQ=255)
                                                                    cctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGAt                                                                        >  1:242506‑M2/72‑127 (MQ=255)
                                                                    cctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGAt                                                                        <  2:242506‑M2/56‑1 (MQ=255)
                                                                    cctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACa                                                           >  2:166630‑M2/72‑140 (MQ=255)
                                                                         tcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACataata                                                      >  1:230689‑M2/67‑140 (MQ=255)
                                                                         tcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACa                                                           <  2:456482‑M2/69‑1 (MQ=255)
                                                                         tcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACa                                                           >  1:456482‑M2/67‑135 (MQ=255)
                                                                           tggAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgag                                                    >  1:660702‑M2/1‑67 (MQ=255)
                                                                             gAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                                                  >  2:181680‑M2/1‑65 (MQ=255)
                                                                               caaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGAc                                                <  1:125313‑M2/80‑1 (MQ=255)
                                                                                  atgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGa                                             >  2:288881‑M2/58‑140 (MQ=255)
                                                                                  aTACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcat                                                       >  1:660166‑M2/1‑60 (MQ=255)
                                                                                  aTACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcat                                                       <  2:660166‑M2/130‑71 (MQ=255)
                                                                                           cacaAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc                                    <  1:600458‑M2/140‑90 (MQ=255)
                                                                                             caAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgc                                  <  1:36972‑M2/140‑92 (MQ=255)
                                                                                             caAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                                                  >  2:189285‑M2/1‑49 (MQ=255)
                                                                                             caAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                                                  <  1:189285‑M2/124‑76 (MQ=255)
                                                                                              tcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGAt                                 <  1:166630‑M2/95‑1 (MQ=255)
                                                                                                    tatTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagt                           >  2:184073‑M2/1‑42 (MQ=255)
                                                                                                    tatTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagt                           >  1:181456‑M2/1‑42 (MQ=255)
                                                                                                    tatTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagt                           >  1:386711‑M2/1‑42 (MQ=255)
                                                                                                       cgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGaaaaaa                        <  2:83804‑M2/104‑1 (MQ=255)
                                                                                                         gACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaatt                      >  2:121751‑M2/1‑37 (MQ=255)
                                                                                                            agcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGa                   <  2:663101‑M2/109‑1 (MQ=255)
                                                                                                                    atTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagc           >  1:477826‑M2/1‑26 (MQ=255)
                                                                                                                        aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagtgccg       <  2:482363‑M2/140‑119 (MQ=255)
                                                                                                                        aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagc           <  1:340691‑M2/136‑115 (MQ=255)
                                                                                                                        aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagc           >  2:340691‑M2/1‑22 (MQ=255)
                                                                                                                             tCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcaga  <  1:286930‑M2/140‑124 (MQ=255)
                                                                                                                                           |                                                                                                                             
TTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTC  >  NC_000913/3765948‑3766212

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