New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3363001 =41 (2.120)33 (1.720) 28/274 0.0 61.9% coding (195/2382 nt) yhcD putative fimbrial usher protein YhcD
?NC_000913 3766087 = 0 (0.000)coding (3905/4134 nt) rhsA rhs element protein RhsA

TTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAAT  >  NC_000913/3765948‑3766222
                                                                                                                                           |                                                                                                                                       
caatatatagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactt                                                                                                                                         <  1:68627‑M2/1‑1 (MQ=255)
    cAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGgg                                                                                                                                     >  1:146879‑M2/1‑138 (MQ=255)
    atatagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGat                                                                                                                                     <  2:196141‑M2/5‑1 (MQ=255)
      atagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGatat                                                                                                                                   <  1:148293‑M2/7‑1 (MQ=255)
      atagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGatat                                                                                                                                   <  1:178471‑M2/7‑1 (MQ=255)
      atagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGatat                                                                                                                                   <  1:1981‑M2/7‑1 (MQ=255)
         gtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAAc                                                                                                                                <  2:3938‑M2/10‑1 (MQ=255)
          tcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACa                                                                                                                               >  2:138067‑M2/130‑140 (MQ=255)
             taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  2:285986‑M2/127‑140 (MQ=255)
             taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  1:309654‑M2/127‑140 (MQ=255)
             taacagcaatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgcgatttactTAGATATAACAATg                                                                                                                            >  1:301593‑M2/127‑140 (MQ=255)
             taacagcaatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  1:136829‑M2/127‑140 (MQ=255)
             ggCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGgga                                                                                                                                    <  2:326553‑M2/132‑4 (MQ=255)
             ggCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGgga                                                                                                                                    >  1:326553‑M2/1‑129 (MQ=255)
               acagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGAt                                                                                                                          >  2:23250‑M2/125‑140 (MQ=255)
                  aaaaGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaag                                                                                                                       >  2:261800‑M2/1‑124 (MQ=255)
                   ccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCt                                                                                                                      <  1:124194‑M2/20‑1 (MQ=255)
                     aGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    >  1:138476‑M2/1‑121 (MQ=255)
                     aGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    >  1:119829‑M2/1‑121 (MQ=255)
                            tcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGa                                                                                                             <  1:285986‑M2/29‑1 (MQ=255)
                             cttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                            <  1:280385‑M2/30‑1 (MQ=255)
                                 cagcaAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    >  2:8946‑M2/1‑109 (MQ=255)
                                 cagcaAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    <  1:8946‑M2/128‑20 (MQ=255)
                                  gaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGtt                                                                                                       <  1:138067‑M2/35‑1 (MQ=255)
                                     ggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCa                                                                                                    <  1:323858‑M2/38‑1 (MQ=255)
                                          ttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGacac                                                                                               <  1:19743‑M2/43‑1 (MQ=255)
                                            ttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAg                                                                                             <  2:276250‑M2/45‑1 (MQ=255)
                                                 aTGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgg                                                                                        >  2:110983‑M2/1‑93 (MQ=255)
                                                      gcAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcca                                                                                   >  1:240794‑M2/1‑88 (MQ=255)
                                                                tGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacaggg                                                                         >  1:309775‑M2/1‑78 (MQ=255)
                                                                  tgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGa                                                                       <  1:23250‑M2/67‑1 (MQ=255)
                                                                      tgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACata                                                                   >  1:183516‑M2/70‑140 (MQ=255)
                                                                       gatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATAACataa                                                                  >  2:17083‑M2/69‑140 (MQ=255)
                                                                         ggtggAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                                                                >  1:255903‑M2/1‑69 (MQ=255)
                                                                         ggtggAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                                                                >  2:194999‑M2/1‑69 (MQ=255)
                                                                          caagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATAACAtaataa                                                               <  1:17083‑M2/75‑1 (MQ=255)
                                                                               aatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGAc                                                          <  2:86117‑M2/80‑1 (MQ=255)
                                                                                aaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACa                                                         <  2:183516‑M2/81‑1 (MQ=255)
                                                                                       gccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTg                                                  <  1:223151‑M2/88‑1 (MQ=255)
                                                                                          tCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacctt                                               <  2:92495‑M2/140‑89 (MQ=255)
                                                                                             caAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagacatagaacagggttcatcatgagtca                                                           <  1:84638‑M2/125‑77 (MQ=255)
                                                                                             caAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagacatagaacagggttcatcatgagtca                                                           >  2:84638‑M2/1‑49 (MQ=255)
                                                                                                                 tgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAAccc                               >  2:148306‑M2/27‑133 (MQ=255)
                                                                                                                 tgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAAccc                               <  1:148306‑M2/107‑1 (MQ=255)
                                                                                                                     gccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGa                    <  2:309654‑M2/118‑1 (MQ=255)
                                                                                                                      ccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGAc                   >  2:222830‑M2/22‑140 (MQ=255)
                                                                                                                        ttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGt                 >  1:305608‑M2/20‑140 (MQ=255)
                                                                                                                        ttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGACCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGt                 <  1:43096‑M2/121‑1 (MQ=255)
                                                                                                                        aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:223432‑M2/140‑119 (MQ=255)
                                                                                                                                cgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAgg         <  2:265001‑M2/129‑1 (MQ=255)
                                                                                                                                       gCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaaga  <  2:223818‑M2/140‑134 (MQ=255)
                                                                                                                                       gCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaaga  <  1:82792‑M2/140‑134 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAAT  >  NC_000913/3765948‑3766222

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.