New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_017635 4805124 =39 (0.960)31 (0.940) 21/80 0.3 56.0% intergenic (+48/‑432) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NC_017635 4805438 = 17 (0.520)intergenic (+362/‑118) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AAAAGTTTGTGCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGT  >  NC_017635/4805391‑4805487
                                               |                                                 
gagaaaaaatacgtaacttatttatgatatggacagtttggccccaaa                                                   <  2:799290‑M2/1‑1 (MQ=255)
aaaaGTTTGTGCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                 >  2:1190611/1‑50 (MQ=255)
 agaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGt                                                >  2:978465‑M2/47‑50 (MQ=255)
 aaaGTTTGTGCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                <  1:728282/50‑1 (MQ=255)
  aaGTTTGTGCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGtt                                               <  1:419976/50‑1 (MQ=255)
    aaaaatacgtaacttatttatgatatggacagtttggccccaaATGttt                                              >  1:1836526‑M2/44‑49 (MQ=255)
         tgCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtt                                        <  2:109601/50‑1 (MQ=255)
         tgCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCt                                         <  2:1530869/49‑1 (MQ=255)
           cgtaacttatttatgatatgtacagtttggccccaaATGTTTCATCtttt                                      >  2:1720966‑M2/37‑50 (MQ=255)
           cgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                      <  2:1115801‑M2/14‑1 (MQ=255)
           cgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                      <  2:34531‑M2/14‑1 (MQ=255)
            gtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTg                                     <  2:1310000‑M2/15‑1 (MQ=255)
            gtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTg                                     >  2:892223‑M2/36‑50 (MQ=255)
            gtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTg                                     <  2:1356648‑M2/15‑1 (MQ=255)
            gtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTg                                     >  1:822708‑M2/36‑50 (MQ=255)
             taacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTgg                                    <  2:391163‑M2/16‑1 (MQ=255)
              gATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                   <  1:1618989/50‑1 (MQ=255)
               acttatttatgatatggacagtttggccccatATGTTTCATCTTTTgggg                                  <  2:1646573‑M2/18‑1 (MQ=255)
               acttatttatgatatggacagtttggccccaaATGTTTCATCTTTTgggg                                  <  2:1321160‑M2/18‑1 (MQ=255)
                cttatttatgatatggacagtttggccccaaATGTTTCATCTTTTggggg                                 <  2:2046650‑M2/19‑1 (MQ=255)
                   atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGaa                               >  1:1094091‑M2/29‑49 (MQ=255)
                   atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGaa                               >  2:743373‑M2/29‑49 (MQ=255)
                   atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGa                                <  1:1268006‑M2/20‑1 (MQ=255)
                   atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGa                                <  1:853841‑M2/20‑1 (MQ=255)
                   atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGa                                <  1:1582471‑M2/20‑1 (MQ=255)
                     tCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                              <  1:1944050/48‑1 (MQ=255)
                     tCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                              <  1:985564/48‑1 (MQ=255)
                      tatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAAc                            <  2:1443841‑M2/24‑1 (MQ=255)
                       atgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACtg                          >  1:1115597‑M2/25‑50 (MQ=255)
                       atgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACtg                          >  1:1984974‑M2/25‑50 (MQ=255)
                        tgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACtgt                         <  1:599029‑M2/27‑1 (MQ=255)
                          atatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGc                       <  2:1402708‑M2/29‑1 (MQ=255)
                            tGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                       <  1:936972/48‑1 (MQ=255)
                             tggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                    <  1:1961338‑M2/32‑1 (MQ=255)
                             gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                    >  1:1870430/1‑50 (MQ=255)
                             gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                    >  2:1676804/1‑50 (MQ=255)
                               gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                  <  2:1903460/50‑1 (MQ=255)
                                tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                 <  1:1417837/50‑1 (MQ=255)
                                acagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                 >  1:805070‑M2/16‑50 (MQ=255)
                                acagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                 <  2:493005‑M2/35‑1 (MQ=255)
                                 cagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                <  1:548087‑M2/36‑1 (MQ=255)
                                 cagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                <  2:801970‑M2/36‑1 (MQ=255)
                                 cAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                >  2:1704412/1‑50 (MQ=255)
                                  agtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg               <  2:715692‑M2/37‑1 (MQ=255)
                                  agtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg               <  1:985272‑M2/37‑1 (MQ=255)
                                  aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGGGTTgg               >  1:1387239/1‑50 (MQ=255)
                                   aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGc              <  2:593103/50‑1 (MQ=255)
                                     aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg            >  2:1463887/1‑50 (MQ=255)
                                      tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt           <  2:860877‑M2/41‑1 (MQ=255)
                                      tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt           >  2:1790036‑M2/10‑50 (MQ=255)
                                      tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt           <  1:867202‑M2/41‑1 (MQ=255)
                                         ccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaa        <  2:2007600‑M2/44‑1 (MQ=255)
                                         ccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa         >  2:1827113‑M2/7‑49 (MQ=255)
                                          cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaaa       <  2:1480313‑M2/45‑1 (MQ=255)
                                           ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc      <  2:53302‑M2/46‑1 (MQ=255)
                                            caaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt     <  1:629622‑M2/47‑1 (MQ=255)
                                            caaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt     >  1:1026176‑M2/4‑50 (MQ=255)
                                             aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTc    <  1:1528338‑M2/48‑1 (MQ=255)
                                              aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCg   >  1:724435‑M2/2‑50 (MQ=255)
                                               aTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt  >  1:1150566/1‑50 (MQ=255)
                                               aTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt  <  1:1549844/50‑1 (MQ=255)
                                               |                                                 
AAAAGTTTGTGCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGT  >  NC_017635/4805391‑4805487

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