New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 454269144 (1.030)8 (0.200) 8/266 3.9 17.2% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913_3_bme_pgi = 4542987 36 (0.900)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913_3_bme_pgi/4542555‑4542838
                                                                                                                                        |                                                                                                                                                   
ccAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCatat                                                                                                                                         >  1:614998/1‑149 (MQ=255)
      aTGCTGCTCCTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACTAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                                                                                                                   <  1:317052/149‑1 (MQ=255)
       tgctgcTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaat                                                                                                                                  >  1:313614/1‑149 (MQ=255)
        gctgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                                                                                                                 >  2:418308/1‑149 (MQ=255)
         ctgctCGTTTTGACGGATTATGGTAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAATTTATTGATAAGAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                                <  1:707931/149‑1 (MQ=255)
            ctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGtt                                                                                                                             <  1:418308/149‑1 (MQ=255)
                     ccgtagtATGGGAAAGAAAAAACCACATAAACGAAAAAGTAAAAAGAGAAGAGGTTTGATGTAAATTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAtttt                                                                                                                    <  2:433949/143‑1 (MQ=255)
                      cGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAttttt                                                                                                                   >  2:47413/1‑149 (MQ=255)
                      cGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGCCCAAACTGTCCATATCATAAATAAGTTACGTAttttt                                                                                                                   >  1:255202/1‑149 (MQ=255)
                         aTTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtct                                                                                                                <  1:47413/149‑1 (MQ=255)
                            ataggAAAGAAATAATCGCCAAAACGATAAATTAAAAAGAGAAGAGGTTTGATTTATCTTATTTATAATAAAAATAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                                             <  2:425409/146‑1 (MQ=255)
                              gggAAAGAAATAATTTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                                                                                           <  2:239567/149‑1 (MQ=255)
                              gggAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                                                                                           <  2:614998/149‑1 (MQ=255)
                                gaaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCAt                                                                                                         >  2:487682/1‑149 (MQ=255)
                                gaaagaaaTAATCTCATAAACGAAAAACAAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCAt                                                                                                         >  2:645290/1‑149 (MQ=255)
                                     aaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAg                                                                                                                               >  2:193441/1‑122 (MQ=255)
                                     aaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAg                                                                                                                               <  1:193441/122‑1 (MQ=255)
                                         aatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                                                                                                >  2:568233/1‑149 (MQ=255)
                                            ctcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaa                                                                                             <  1:487682/149‑1 (MQ=255)
                                               aTAAACGAAAAATTAAAAAGAGAAAAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaaa                                                                                          >  1:370750/1‑149 (MQ=255)
                                                           ttAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCa                                                                              >  1:96582/1‑149 (MQ=255)
                                                            tAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt                                                                             <  1:710252/149‑1 (MQ=255)
                                                              aaaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                                                                           <  2:313614/149‑1 (MQ=255)
                                                                 agagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAACTACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAc                                                                        >  1:703417/1‑149 (MQ=255)
                                                                        agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGa                                                                 <  1:586614/149‑1 (MQ=255)
                                                                          gTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCAAAAATCTCAATATCATAAATAAGTTAAGTATTATTTCTCAATCATACAAATATTAAAAAACGACAAAAAGCTTATAATCGTTTgctt                                                               >  2:478174/1‑146 (MQ=255)
                                                                               aTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTaa                                                          >  1:646875/1‑149 (MQ=255)
                                                                                       tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTc                                                  >  1:685938/1‑149 (MQ=255)
                                                                                         ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTa                                                <  1:645290/149‑1 (MQ=255)
                                                                                          tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAt                                               <  2:685938/149‑1 (MQ=255)
                                                                                                     tcaaaaaaataaaaaaatATAAGACAATAGGGGACAAACTGTTCATATCATAAATAAGTTACGTATTTTGTCTCAAGCATAACAATATTAAAAAACAACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaattaat                                    <  2:467073/147‑1 (MQ=255)
                                                                                                       aaaaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAAttt                                  <  1:682147/149‑1 (MQ=255)
                                                                                                       aaaaaaaCAAATAAATACAAGAAAATTAGGGCCAAACTGTCCATATCAGAAATAAGTTACGTATTTTTTCTCACGCATAAGAATATTAAAAAACGACAAAAAGTATCTAACTATTTGATATGTAAATTAACTATATTGTAAATTAAgtt                                  >  1:120412/1‑149 (MQ=255)
                                                                                                           aaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcaca                              <  2:646875/149‑1 (MQ=255)
                                                                                                           aaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcaca                              <  2:667389/149‑1 (MQ=255)
                                                                                                               aaataaatACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAc                          >  1:70124/1‑149 (MQ=255)
                                                                                                                         gaGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtata                <  2:712656/148‑1 (MQ=255)
                                                                                                                         aaGACAATTGGGGCCAAACTGTCCATATCATAAATACGTTTCGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtata                <  1:206431/149‑1 (MQ=255)
                                                                                                                          aaaCATTTGAGGCCAAACGGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTatat               <  2:490510/147‑1 (MQ=255)
                                                                                                                          aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTatat               <  2:96582/149‑1 (MQ=255)
                                                                                                                               gttGGGGACAAACTGTCAATATCAAAAATAAGTTACGTATTTTTACTCAAGTATAAAAATATTAAAAAACAAAAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTaaa          <  2:255202/148‑1 (MQ=255)
                                                                                                                                 tGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTTAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCAGCACTAGGTaaatc        <  2:196606/149‑3 (MQ=255)
                                                                                                                                     gCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAAc    >  1:481422/1‑149 (MQ=255)
                                                                                                                                       cAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAACTCCGCTATTTGTAAAGCTAACGt  >  1:286812/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
CCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913_3_bme_pgi/4542555‑4542838

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.